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ABSTRACT 


The  operation  of  an  aircraft  antiskid  wheel  brake 
control  system  has  the  potential  for  producing  adverse 
aircraft  dynamic  behavior  and  structural  damage.  Antiskid 
operation  is  also  a  major  influence  upon  stopping  perfor¬ 
mance.  Unless  the  characteristics  and  effects  of  antiskid 
operation  can  be  defined,  an  aircraft's  capability  for 
safe,  reliable  and  economical  accomplishment  of  its  inten¬ 
ded  usage  cannot  be  assured.  This  report  presents  an 
analysis  procedure  for  predicting  antiskid  operational 
character  i.stics  and  the  inter-related  effects  upon  the 
aircraft  and  its  performance.  The  analytical  procedure 
is  the  development  of  mathematical  equations  for  a  com¬ 
prehensive  description  of  the  antiskid  system  components, 
the  significantly  influencing  aircraft  systems  and  the 
characteristics  of  the  surface  upon  which  the  aircraft  is 
operating.  The  mathematical  description  includes  such  con¬ 
siderations  as  landing  gear  dynamics,  tire  elasticity, 
brake  torque  response  characteristics,  antiskid  electronic 
circuitry,  brake  hydraulic  control  system  dynamics,  runway 
surface  profile  and  tire-to-runway  friction  characteristics. 
Both  on-off  and  "modulated**  antiskid  systems  are  analyzed. 
Procedures  for  quantitative  evaluation  of  the  influencing 
parameters  and  examples  of  their  usage  ai.e  also  presented. 
The  implementation  of  the  analytical  prediction  procedure 
by  simultaneous  solution  of  all  the  mathematical  equations 
on  an  electronic  computer  is  described. 
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SECTION  I 


INTRODUCTION 

An  antiskid  system  Is  provided  as  a  part  of  the  landing 
gear  wheel  brake  control  system  of  most  large  aircraft, 
particularly  those  having  full  power  brake  actuation. 
Aircraft  operational  experience  has  shown  that  an  anti¬ 
skid  system  is  required  because  there  are  many  occasions 
where  the  maximum  available  friction  force  between  the 
tires  and  runway  surface  is  insufficient  to  react  the 
applied  brake  torque.  For  c^ses  where  excessive  brake 
torque  is  applied  the  antiskid  system  functions  to  control 
tire  motion  so  that  ^kids  are  prevented  and  so  that  the 
associated  problems  and  hazardous  circumstances  which  are 
detrimental  to  safe,  predictable  and  economical  aircraft 
operation  are  avoided.  The  antiskid  function  is  accom¬ 
plished  by  a  group  of  ancillary  components  which  provide 
an  automatic  means  for  detecting  a^'d  alleviating  an 
incipient  tire  skid  condition  by  controlling  brake  torque. 

An  incipient  skid  is  alleviated  by  temporarily  reducing 
brake  torque  to  a  value  less  than  the  torque  being  produced 
by  the  friction  force  at  the  tire-runway  interface.  Brake 
torque  reduction  is  sustained  for  a  time  interval  of  suf¬ 
ficient  duration  to  allow  the  wheel  to  regain  speed.  After 
Che  wheel  has  regained  speed,  brake  torque  is  reapplied. 

The  reduction  and  subsequent  reapplication  of  brake 
torque  results  in  an  oscillatory  braking  force  being  applied 
to  the  airplane.  This  oscillatory  force  has  the  potential 
for  causing  adverse  dynamic  loading  of  the  airplane  struc¬ 
ture,  for  causing  directional  control  difficulty  and  for 
degrading  the  aircraft's  stopping  performance.  Therefore, 
the  antiskid  system  must  control  tire  motion  in  a  way  such 
that  objectionable  or  unsafe  conditions  other  than  those 
related  to  tire  skidding  are  not  incurred.  The  need  for 
evaluating  the  potentially  deleterious  effects  of  an  oscil¬ 
latory  braking  force  is  now  recognized  because  there  have 
been  a  number  of  instances  where  failure  to  do  so  has 
resulted  in  severe  operational  difficulty  and  in  some  cases 
catastrophic  landing  gear  failure. 

The  objective  of  this  study  is  to  develop  analytical 
procedures  and  techniques  for  predicting  aircraft  antiskid 
operational  behavior  and  its  effects.  These  analysis 
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techniques  are  intended  to  help  overcome  some  of  the  pre¬ 
viously  experienced  problems  or  uncertainties  and  to 
provide  a  foundation  for  a  comprehensive  evaluation  of 
aircraft  antiskid  performance  and  total  system  compati¬ 
bility.  It  is  also  intended  that  these  procedures  be 
capable  of  application  during  the  conceptual  design  phase 
of  new  airplanes.  In  the  initial  design  of  a  new  airplane 
the  capabilities  of  various  candidate  equipment  which  might 
be  used  for  stopping  during  the  landing  sequence  or 
rejected  takeoff  should  be  evaltiated  with  respect  to  the 
airplane's  mission  requirements.  Factors  such  as  stopping 
performance,  weight,  cost  and  reliability  should  be  consi¬ 
dered  when  the  influence  of  the  braking  equipment  is  being 
examined  to  establish  the  overall  effect  upon  the  aircraft's 
configuration.  In  such  an  evaluation,  the  performance  of 
the  wheel  braking  system,  including  any  applicable  antiskid 
equipment,  is  a  major  consideration.  Use  of  an  analysis 
procedure  whereby  the  effects  of  antiskid  operation  can  be 
accurately  predicted  provides  the  means  for  minimi'^.ing  the 
technical  and  financial  risks  of  both  the  aircraft  manu¬ 
facturer  and  the  aircraft  user.  Inaccurately  predicting 
the  wheel  braking  system's  performance  can  result  in  an 
airplane  design  unsuited  for  its  intended  usage,  a  costly 
redesign  program,  or  both. 

This  study  mathematically  describes  the  physical  oper¬ 
ation  of  antiskid  equipment  in  conjunction  with  the  airplane 
and  its  other  applicable  components.  The  basis  of  the 
mathematical  relationships  is  the  description  of  actual 
(or  conceivable)  hardware  behavior  rather  than  a  compilation 
of  equations  relating  various  parameters  in  a  desirable  or 
compatible  manner  without  regard  to  detail  design  features. 
This  approach  is  taken  to  assure  all  influencing  parameters 
are  accounted  for  and  to  provide  criteria  for  equipment 
detail  design  and  test.  Also,  by  examining  the  individual 
component  behavior,  the  evaluation  can  include  such  consid¬ 
erations  as  cost  and  weight  along  with  performance  charac¬ 
teristics. 

The  essence  of  antiskid  operation  is  the  cumulative 
effect  of  a  number  of  successive  events,  where  the  inter¬ 
vening  occurrences  and  outcome  of  each  is  influenced  by  and 
dependent  upon  the  conditions  resulting  from  preceding 
events.  Since  these  events  occur  quite  rapidly  and  involve 
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the  behavior  of  the  aircraft  and  many  of  its  components, 
the  instantaneous  condition  of  a  very  large  number  of 
variables  must  be  continually  maintained  with  high  accuracy 
so  that  they  are  available  when  needed.  Consequently,  one 
of  the  major  problems  associated  with  analyzing  antiskid 
operation  is  the  magnitude  of  the  computation  task.  It 
will  be  noted  that  the  study  has  analytical  components 
encompassing  several  engineering  and  scientific  disciplines 
such  as  electronics,  aerodynamics,  mechanics  and  hydraulics. 
Each  of  the  individual  analytical  components  is  often 
deserving  of  considerable  more  elaborate  and  complete 
treatment.  However,  to  provide  an  economically  feasible 
and  comprehensible  composite  solution,  the  scope  of  the 
individual  analytical  components  has  been  limited  to  account 
for  only  those  effects  or  influencing  factors  which  are  of 
traditional  interest  and  which  are  required  to  achieve 
reasonable  agreement  between  observed  operational  behavior 
and  analytical  results. 
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SECTION  II 


ANALYTICAL  APPROACH 

The  analytical  approach  of  this  study  is  directed 
toward  predicting  the  existence  of  adverse  circumstances 
which  have  caused  various  problems  in  the  past  and  t-oward 
providing  information  which  is  typically  needed  to  estab¬ 
lish  detail  design  criteria  and  to  define  aircraft  operating 
procedures.  Specific  consideration  is  given  to  providing 
the  means  for: 

(a)  Establishing  the  magnitude  and  frequency  of  dynamic 
loading  applied  to  the  landing  gear. 

(b)  Establishing  the  value  of  the  braking  force  which  can 
be  predictably  and  dependably  achieved  for  various 
runway  surface  and  aircraft  operating  conditions. 

(c)  Determining  individual  component  a.id  system  operational 
characteristics  which  are  required  so  that  overall 
aircraft  performance  objectives  are  achieved. 

(d)  Establishing  th«i  effects  of  varying  performance 
characterisrics  of  individual  components  within  the 
brake  control  system  to  assure  no  incompatibilities 
exist. 

1.  PROBLEM  DEFINITION 

Figure  1  is  a  block  diagram  showing  the  typical 
arrangement  of  an  antiskid  system  and  its  relationship 
within  the  total  aircraft  oystem.  This  arrangement  is 
representative  of  most  antiskid  systems  in  current  use 
and  the  various  types  of  airplanes  on  which  they  are 
installed.  The  major  components,  the  significant  forces 
and  their  controlling  elements  are  shown  for  a  single 
wheel  main  gear  configuration  of  a  nose  wheel  type  air¬ 
plane  which  is  the  usual  case  for  fighter  type  aircraft. 

For  airplanes  having  multiple  wheeled  landing  gears  and/or 
multiple  landing  gears  the  same  basic  relationships  prevail 
with  the  addition  of  similar  type  components  as  appropriate. 
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Antiskid  systems  usually  operate  by  measuring  a 
wheel's  motion,  comparing  the  measurement  to  an  index 
of  acceptability  and  causing  broke  torque  to  be  decreased 
or  increased  in  accordance  with  some  function  of  the  dif¬ 
ference  between  the  measured  motion  and  comparison  index. 

A  detailed  description  of  the  operational  behavior  and 
influence  of  the  individual  elements  is  presented  in 
Section  111. 

Since  antiskid  operation  is  basically  the  control  cf 
tire  motion  and  since  the  motion  of  a  tire  is  determined 
by  the  forces  imposed  (the  same  as  for  any  other  object) 
the  study  of  antiskid  operation  resolves  itself  into  (1) 
defining  the  forces  on  the  tire  and  wheel  and  (2)  estab¬ 
lishing  the  resultant  effects  of  these  forces.  It  is 
easily  observed  that  the  forces  acting  upon  an  airplane 
tire  and  wheel  are  the  forces  between  the  tire  tread  and 
runway  surface  and  the  forces  from  the  airplane's  landing 
gear  and  brake.  The  values  of  these  forces  are  established 
by  tho  wheel's  relative  position  and  relative  motion  with 
respect  to  the  runway  surface  and  to  the  airplane.  The 
wheel's  relative  motion  and  position  is  determined  by  con¬ 
sidering  simultaneous  and  interrelated  actions  of  the 
aircraft  and  a  number  of  its  systems.  The  effects  of  the 
following  parameters  are  considered  in  this  study. 

(a)  Tire  circumferential  deformation  and  its  rate 

(b)  Tire  radial  deformation  and  its  rate 

(c)  Brake  torque  as  a  function  of  velocity,  the  brake's 
inertia,  and  actuation  pressure 

(d)  Brake  actuation  pressure  as  a  function  of  the 
actuation  media's  compressibility  and  inertia,  line 
restrictions  and  elasticity,  variable  flow  areas 
within  valves  and  the  actuation  media's  containment 
vessels'  (lines,  brake  housing,  valve  bodies)  volume 

(e)  Elastic  and  inertia  properties  of  the  landing  gear 

(f)  Aerodynamic  forces  upon  the  airplane 

(g)  Runway  surface  profile 

(h)  Tlre-to-runway  friction  coefficient  as  a  function  of 
relative  velocity  and  runway  surface  condition 
Including  hydroplaning  effects 

(1)  The  aircraft's  inertia  and  control  surface  position 
including  stability  augmentation  system  effects. 
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2 .  BACKGROUND 


During  the  initial  design  and  system  development  phase 
for  most  new  aircraft,  it  has  become  a  customary  practice 
to  analyze  antiskid  operation  to  define  its  effects  and 
thereby  assure  compliance  with  the  airplane's  stopping 
performance  objectives  and  assure  adverse  dynamic  loading 
conditions  or  directional  control  problems  will  not  be 
encountered.  These  analyses  have  usually  been  accomplished 
by  utilizing  a  set-up  composed  of  hardware  representative 
of  aircraft  components  interfaced  with  an  electronic  com¬ 
puter  (most  often  an  analog  computer).  The  computer  is 
used  to  solve  mathematical  equations  describing  the  motion 
of  the  aircraft  and  the  landing  gear,  forces  on  the  air¬ 
craft,  tire  and  wheel  motion  and  tire-to-runway  friction, 
etc.  The  actual  behavior  of  a  laboratory  set-up  including 
such  components  as  the  antiskid  control  circuit,  hydraulic 
brake  valves  and  interconnecting  lines  is  measured  by 
suitable  instrumentation  and  fed  into  the  computer  to  obtain 
a  composite  solution.  This  analysis  procedure  is  used 
because  a  complete  mathematical  computer  setup  requires 
greater  computer  capacity  than  is  usually  available  and 
because  an  accurate  mathematical  description  for  some  com¬ 
ponents  such  as  the  electronic  antiskid  control  circuit  is 
often  unavailable. 

Some  antiskid  analyses  have  been  performed  using  an 
"all  mathematical"  approach;  however,  these  have  usually 
been  associated  with  academic  endeavors  or  a  comparative 
evaluation  of  a  specific  device  and  did  not  account  for 
all  of  the  known  significant  influencing  parameters  and 
constraints  for  an  actual  aircraft  antiskid  sytem  instal¬ 
lation.  While  the  hybrid  hardware -computer  analyses  have 
often  satisfied  their  objectives,  several  factors  have  led 
to  a  number  of  uncertainties  for  which  the  bounds  are  not 
adequately  established,  either  because  of  great  difficulty 
and  expense  or  because  of  inadequate  knowledge.  These 
uncertainties  tend  to  obscure  the  analysis  results  and 
generally  detract  from  their  credibility.  The  most  signi¬ 
ficant  factor  causing  uncertainty  is  that  the  usual 
definition  for  the  friction  force  between  the  tire  and 
runway  surface  does  not  account  for  all  the  observed 
variations.  A  second  factor  is  the  analytical  limitations 


associated  with  the  u^se  of  actual  hardware.  The  use  of 
actual  hardware  dictates  that  the  analysis  be  performed 
"real  time"  and  complicates  or  prevents  examination  of 
some  parameter  variations.  Since  some  parameters  have 
a  very  high  rate  of  variation  with  res;'^ct  to  time,  the 
outputs  from  a  "real  time"  solution  can  ue  extremely 
difficult  to  observe  and  interpret.  Also,  the  instru¬ 
mentation  used  to  interface  the  hardware  with  the  computer 
introduces  additional  variables  to  an  otherwise  very 
complex  system.  This  study  is  intended  to  provide  the 
means  for  overcoming  these  problems  and  for  minimizing 
uncertainty. 

3.  ANALYTICAL  PROCEDURE  AND  RATIONALE 

The  evaluation  of  antiskid  operation  is  conducted 
using  a  modular  analysis  technique  whereby  the  problem  is 
divided  into  a  number  of  modules  or  component  parts,  each 
having  defined  inputs  and  outputs  such  that  the  outputs 
from  one  or  more  components  are  provided  as  inputs  to  other 
components.  By  combining  all  the  analytical  components,  a 
composite  simultaneous  .lolution  is  obtained.  The  analytical 
modules  are  formulated  so  as  to  correspond  to  various  air¬ 
craft  components  or  systems.  The  modules  can  be  arranged 
in  a  number  of  combinations  representative  of  a  variety  of 
aircraft  configurations.  In  addition,  the  modular  approach 
allows  maximum  computation  flexibility  in  that  changes  can 
be  made  within  individual  modules  without  affecting  the 
overall  analysis  program.  The  predominate  influencing 
factors  governing  the  choice  of  each  analytical  component's 
content  and  treatment  are  experience  and  judgment  as  to  the 
degree  of  detail  which  is  required  to  accurately  establish 
the  timing  or  relative  sequence  of  significant  events,  bach 
analytical  module  is  formulated  so  that  particular  effects 
or  circumstances  can  be  examined  and  so  that  its  outputs 
will  supply  the  information  needed  as  inputs  to  other  modules. 
It  will  be  noted  that  some  relatively  insignificant  para¬ 
meters  must  be  considered  to  achieve  mathematical  continuity. 
To  exemplify  the  analysis  procedure  antiskid  operation  for 
a  fighter  type  aircraft  having  a  single  wheel  main  landing 
gear  arrangement  is  evaluated.  All  of  the  analytical  com¬ 
ponents,  except  for  the  antiskid  control  circuit,  are 
expressed  in  general  terms  and  could  be  applied  to  almost 
any  airplane.  The  antiskid  control  circuits  considered  are 
those  specifically  utilized  on  the  F104  and  the  F-111. 
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For  the  case  of  the  F-'104  on-o££  antiskid  control  circuit, 
the  wheel  speed  input  signal  is  arbitrarily  adjusted  to 
account  for  the  difference  between  the  F-104  and  F-111 
tire  sizes.  All  parameter  values  used  to  prove  the  validity 
of  the  analysis  procedures  are  those  associated  with  the 
F-111  airplane  so  that  the  analytical  results  can  be  com¬ 
pared  to  available  records  of  actual  aircraft  operation. 

To  analyze  other  control  circuits  will  require  that  their 
mathematical  models  be  formulated  and  incorporated  in  the 
composite  solution.  The  detail  assumptions  and  procedures 
for  establishing  parameter  values  are  presented  in 
Section  IXI  within  the  description  of  each  analytical 
module. 


4 .  parameter  INv'ESTIGATIONS 

The  basic  intent  of  this  study  is  to  account  for  the 
influence  of  parameters  and  effects  which  have  been 
identified' as  responsible  for  previously  experienced 
operational  difficulties  or  which  are  otherwise  known 
to  significantly  affect  antiskid  performance.  Such 
items  as  tire  radial  and  circumferential  spring  rate, 
the  characteristics  of  brake  torque  variations  with 
velocity  and  actuation  pressure,  brake  chatter  and 
squeal,  hydraulic  system  response  as  affected  by  line- 
sizes,  component  flow  restrictions  and  metering  valve 
characteristics,  the  airplane's  response  to  aerodynamic 
forces  and  runway  roughness,  landing  gear  elastic 
characteristics  and  the  characteristic  of  the  tire-to- 
runway  friction  force  variations  are  given  particular 
attention.  The  treatment  of  most  parameters  is  that 
which  experience  has  proven  gives  satisfactory  results. 
However,  to  overcome  some  previous  antiskid  evaluation 
analytical  difficulties  associated  with  tire-to-runway 
friction  and  hydraulic  system  operation  and  to  examine 
the  effects  of  brake  chatter  and  squeal,  some  prelimi¬ 
nary  investigations  were  conducted. 

A.  Brake  Investigation 

Since  an  antiskid  system  controls  brake  torque 
implicitly  by  controlling  brake  application  pres¬ 
sure,  the  hysteresis  in  the  brake's  torque  response 
to  pressure  changes  must  be  accounted  for.  This 
hysteresis  results  from  inertia  of  the  brake  moving 
parts, friction  forces  on  the  actuating  pistons  due 
to  hydraulic  seals  and  piston  tide  loading,  and 
from  friction  in  the  splined  connections  between 
the  brake  discs  and  the  wheel  and  between  the 
discs  and  the  torque  tube.  To  evaluate  a  typical 
brake's  torque  response  to  rapidly  changing  actua¬ 
tion  pressure  and  to  briefly  investigate  brake 
chatter  and  squeal  effects,  a  relatively  complex 
six-degree  of  freedom  brake  mathematical  model  was 
initially  formulated.  In  this  model  six  discs  were 
treated  as  separate  masses  with  individual  axial 
position,  velocity  and  acceleration  computation, 
non-linear  keyway  and  piston  friction  as  a  function 
of  axial  velocity,  non-linear  brake  lining  friction 
as  a  function  of  rotational  velocity,  and  variable 
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elasticity  to  simulate  the  effects  of  disc  warpage. 

The  model  was  set  up  on  an  analog  computer  and  sub¬ 
jected  to  step  input  pressures  and  to  sinusoidal 
pressure  oscillations  of  various  amplitudes  and  mean 
values  at  frequencies  from  10  cps  to  1000  cps.  The 
computer  setup  also  included  rotational  and  longitu¬ 
dinal  elastic  deformations  within  the  tire  and  brake 
supporting  structure.  The  set  up  was  operated  at 
l/lOO  real  time  and  at  a  number  of  aircraft  veloci¬ 
ties.  By  suitable  choice  of  elastic,  damping  and 
friction  characteristics,  both  chatter  and  squeal 
were  produced  at  low  aircraft  speed.  Using  a  key¬ 
way  friction  coefficient  varying  from  0.15  at  zero 
velocity  to  0.10  at  high  velocity,  it  was  found  that 
the  brake  torque  oscillated  in  response  to  oscillat¬ 
ing  pressure  at  all  frequencies  up  to  1000  cps.  At 
low  brake  rotational  velocities  (20-40  rad/sec)  with 
low  frequency  pressure  oscillation  where  the  minimum 
pressure  was  the  value  for  full  brake  release,  the 
brake  torque  oscillation  had  considerable  deviation 
from  a  sinusoidal  variation.  The  phase  lag  between 
instants  of  maximum  torque  and  maxlumm  pressure  varied 
from  15-20  degrees  at  10  cps  to  40-50  degrees  at  100 
cps  to  110-150  degrees  at  1000  cps.  The  oscillatory 
component  of  the  brake  torque  exhibited  appreciable 
attenuation  at  high  frequency  such  that  the  amplitude 
at  1000  cps  was  about  20  percent  of  the  10  cps  ampli¬ 
tude  with  constant  pressure  amplitude.  Even  though 
there  was  noticeable  phase  lag  in  the  pressure-torque 
characteristic,  it  was  found  that  throughout  the  10- 
1000  cps  frequency  range  there  was  no  appreciable 
phase  difference  between  the  displacement,  velocity 
or  acceleration  of  the  individual  discs.  Therefore, 
a  simplified  model  was  formulated  where  all  the  discs 
were  treated  as  a  single  mass.  The  simple  model  was 
set  up  and  tested  on  the  analog  computer  where  its 
torque  response  to  varying  pressure  was  confirmed  to 
be  identical  to  the  more  complex  model.  The  more 
simple  brake  mathematical  model  is  used  in  this  study 
and  is  described  in  Section  III.  A  significant  and 
somewhat  unexpected  finding  of  this  investigation  is 
that  a  typical  airplane  brake  can  be  expected  to  have 
appreciable  torque  response  when  subjected  to  pres¬ 
sure  oscillations  in  the  100-200  cps  frequency  range 
as  might  be  associated  with  a  hydraulic  line  resonance. 


B*  Tire-to-Runway  Friction  Investigation 


The  usual  and  relatively  arbitrary  function  relat¬ 
ing  coefficient  of  friction  to  tire  or  wheel  slip 
ratio  has  been  used  in  most  prior  antiskid  analyses 
to  establish  the  tire-to-runway  friction  force. 

While  there  are  many  circumstances  where  the  slip 
ratio  approach  is  adequate  for  examining  most  of 
the  aspects  of  antiskid  operation,  a  number  of 
difficulties  and  undesirable  effects  are  associated 
with  its  use,  A  major  analytical  problem  is  that 
examination  of  antiskid  operation  at  low  aircraft 
speed  is  prohibited  because  the  slip  ratio  compu¬ 
tation  would  require  division  by  zero.  In  addition, 
the  large  differences  in  the  friction  coefficient - 
slip  ratio  characteristic  variation  which  have  been 
observed  for  changes  such  as  aircraft  speed,  r  ..iway 
surface  condition  and  tire  properties  lead  to  a 
number  of  uncertainties,  particularly  with  respect 
to  stopping  performance  predictions. 


To  satisfy  the  objectives  of  this  study,  it  was 
considered  necessary  that  a  mathematical  description 
of  the  tire-to-runway  friction  coefficient  be  used 
which  would  not  have  the  above  undesirable  qualities. 
To  develop  such  a  description,  several  hypotheses  were 
formulated  considering  the  tire's  elastic  deformation 
and  its  response  to  ground  friction  forces.  Effects 
such  as  tread  stretch,  tread  circumferential  displace¬ 
ment  and  variation  of  relative  velocity  between  tire 
tread  particles  and  the  runway  surface  throughout  the 
footprint  were  examined  mathematically.  Because  of 
the  extremely  complex  nature  of  a  tire's  elastic 
behavior,  these  examinations  quickly  lead  to  an 
analytical  task  at  least  equal  to  the  scope  of  the 
entire  antiskid  study.  Even  though  this  subject 
deserves  further  investigation,  a  more  simple  hypo¬ 
thesis  accounting  for  most  known  variations  and 
effects  was  adopted  to  comply  with  this  program's 
objectives.  For  the  purpose  of  this  analysis,  it 
is  assumed  that:  (1)  the  tire  tread  is  a  perfectly 
flexible  inelastic  belt  with  radial  and  torsional 
elastic  attachment  to  the  wheel,  (2)  All  tread 
particles  within  the  footprint  have  the  same  rela¬ 
tive  velocity  with  respect  to  the  runway  surface 
and  the  coefficient  of  friction  between  the  tire 
tread  and  runway  surface  Is  a  function  of  relative 
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velocity.  (3)  The  function  defining  the  friction 
coefficient  variation  with  relative  velocity  is  that 
established  by  testing  a  tire  in  a  full  skid. 

A  description  of  the  tire  and  wheel  mathematical 
model  utilizing  these  assumptions  is  contained  in 
Section  III.  The  equations  listed  show  that  the 
relative  velocity  between  the  tire  footprint  and 
runway  surface  is  determined  by  computing  the  tread 
belt's  C.  G.  (center  of  gravity)  translational 
velocity  component  parallel  to  the  runway  surface 
and  the  angular  velocity  of  a  point  on  the  tread 
belt  about  the  C.  G.  The  footprint  horizontal  velo¬ 
city  component  relative  t #  the  C.  G.  is  computed  from 
the  angular  velocity  and  an  apparent  rolling  radius. 
The  apparent  rolling  radius  is  the  unbraked  rolling 
radius  plus  a  fraction  of  the  tread  belt's  C.  Go 
horizontal  displacement  with  respect  to  the  wheel's 
rotational  axis.  The  net  footprint  velocity  rela¬ 
tive  to  the  runway  surface  is  then  the  sum  of  the 
tread  belt  C.  G.  translational  velocity  and  the 
velocit>  of  the  footprint  relative  to  the  tread  belt 
C.  G.  The  mathematical  expression  for  friction 
coefficient  as  a  function  of  relative  velocity  is  of 
exponential  form  vith  coefficients  chosen  to  fit  test 
data. 

This  model  was  set  up  on  an  analog  computer  and 
examined  statically  and  dynamically.  Statically,  the 
friction  coefficient  versus  slip  ratio  (with  respect 
to  the  wheel)  characteristic  varies  with  axle  velocity 
in  accordance  with  observations.  This  observed 
variation  is  that  the  slip  ratio  value  associated  with 
maximum  friction  coefficient  is  greater  at  low  axle 
velocity  than  at  high  axle  velocity,  and  the  value 
of  friction  coefficient  at  maximum  slip  ratio 
decreases  as  axle  velocity  increases.  Figure  2A 
shows  friction  coefficient  versus  slip  ratio  (with 
respect  to  the  wheel)  recorded  dynamically  during  an 
analog  computer  run  with  an  ON-OFF  antiskid  system. 
Figure  2B  ic  a  similar  curve  recorded  dynamically 
during  wheel  spinup  from  a  full  skid.  For  both  cases 
shown  on  Figure  2  axle  velocity  is  constant. 


C ,  Hydraulic  System  Investigation 

From  experience  gained  in  conjunction  with  practi- 
cally  all  antiskid  development  programs,  it  is 
generally  accepted  that  one  o£  the  more  predominate 
influences  upon  antiskid  operation  and  aircraft 
stopping  performance  is  the  time  lag  between  the 
antiskid  control  device's  command  for  a  brake  tor¬ 
que  change  and  the  actual  brake  torque  response. 
Hydraulic  flow  restrictions  and  the  response  charac¬ 
teristics  of  che  antiskid  control  valve  and  other 
hydraulic  system  elements  are  responsible  for  most 
of  this  time  lag.  In  an  attempt  to  minimize  the 
effects  of  the  time  lag  many  antiskid  control  de¬ 
vices  actually  issue  commands  in  anticipation  of  a 
predicted  circumstance.  Confident  prediction  of 
antiskid  overall  operational  effects  including  the 
resultant  airplane  stopping  performance  requires 
that  the  hydraulic  time  lag  be  accurately  accounted 
for.  Therefore,  to  comply  with  the  objectives  of 
this  study,  a  preliminary  exploration  was  conducted 
to  establish  a  suitable  mathematical  model  permitting 
evaluation  of  antiskid  control  valve  and  pilot's 
metering  valve  response  characteristics  and  such 
effects  as  hydraulic  line  resonant  oscillation. 

During  these  explorations  the  operation  of  the 
pilot's  metering  valve,  antiskid  control  valve  and 
the  hydraulic  line  connecting  the  control  valve  to 
the  brake  were  examined.  In  each  case  several 
different  mathematical  descriptions  were  formulated 
and  investigated  on  an  analog  computer. 

For  both  the  pilot's  metering  valve  and  antiskid 
control  valve  mathematical  descriptions  accounting 
for  all  component  characteristics  of  an  actual  phy¬ 
sical  device  and  simpler  descriptions  eliminating 
spool  mass  considerations  were  examined.  While  by 
suitable  choice  of  parameter  values  either  mathe¬ 
matical  model  can  produce  an  accurate  description, 
the  second  order  equations  resulting  from  consi¬ 
deration  o£  spool  mass  cause  analytical  difficulty 
because  the  inertia  is  very  small  in  comparison 
with  hydraulic  pressure  and  spring  forces.  These 
very  high  gain  second  order  systems  necessitate  very 
rapid  integration;  therefore,  using  the  "massless" 
first  order  equations  is  highly  desirable  to  achieve 
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computation  economy-  In  Section  III  the  pilot's 
metering  valve  description  (a  part  of  the  hydraulic 
system)  is  the  simpler  first  order  system  while  the 
control  valve  equations  account  for  spool  mass- 
This  approach  is  taken  to  permit  easy  recognition 
of  the  relationship  between  the  control  valve's 
physical  construction  and  its  performance  charac¬ 
teristics-  While  having  the  same  facility  for  the 
metering  valve  is  desirable,  it  was  considered  ana¬ 
lytically  too  extravagant.  A  metering  valve  having 
satisfactory  performance, by  whatever  physical  means 
it  is  achieved,  will  exhibit  behavior  in  accordance 
with  the  "massless"  equation- 

To  explore  hydraulic  line  resonant  oscillation  and 
"water  hammer"  effects,  a  ten  element  hydraulic 
line  model  (ten  degree  of  freedom)  was  initially 
formulated  and  examined  on  an  analog  computer  with 
On-Off  antiskid  operation  at  one  hundredth  real 
time-  This  model  produced  very  excellent  results; 
however,  the  low  intensity  of  the  higher  frequency 
harmonics  (above  100  cps)  showed  that  a  more  sim¬ 
plified  model  would  probably  be  satisfactory- 
Accordingly,  a  single  degree  of  freedom  model  was 
formulated  and  tested  in  the  same  manner  as  the 
ten  element  model.  For  the  purpose  of  antiskid 
evaluation,  the  single  degree  of  freedom  model  gave 
satisfactory  results  and  is  described  in  Section  III, 


SECTION  III 


DEVELOPMENT  OF  MATHEMATICAL  MODELS 


This  section  Is  devoted  to  the  exposition  of 
caathematlcal  models  for  each  of  tLe  following  total 
system  components: 


1.  Brake  System 

2.  Hydraulic  System 

3.  Airplane  System 

4.  Wheel  and  Tire  System 

5.  Wheel  Speed  Sensor 

6.  Antiskid  Control  Circuit 

7.  Antiskid  Control  Valve 

8.  Horizontal  Tail  Control 

9.  Runway  System 

For  some  of  the  system  components  alternate  models  are  pro¬ 
vided,  These  alternate  models  ar^  listed  alphabetically 
within  each  section.  For  example,  3a  describes  an  airplane 
system  modeled  as  a  laboratory  flywheel,  3b  describes  an 
airplane  which  has  three  degrees  of  freedom,  and  3c  describes 
an  airplane  with  six  degrees  of  freedom.  Each  component 
model  is  discussed  as  a  self-contained  unit  without  any 
particular  reference  to  the  total  system  and  each  model,  in 
general,  contains  its  complete  mathematical  description  such 
that  it  is  essentially  immune  to  changes  within  other  models 
of  the  total  system. 

Format  and  Convention  Useage 

The  presentation  of  the  various  sytems  follows  a  common 
format.  Each  system  discussion  begins  with  an  introductory 
explanation  of  its  function  or  its  characteristics  relevant 
to  antiskid  operation.  Following  this  introduction  is  the 
main  bedy  of  the  discussion  under  the  heading,  "A,  Mathemat¬ 
ical  Description,"  containing  the  derivation  of  the  equations 
that  describe  the  system  dynamically.  This  section  is  con¬ 
cluded  with  an  equation  flow  diagram  showing  the  relationship 
among  the  various  system  equations.  A  final  discussion 
follows  under  the  heading,  "B,  Parameter  Evaluation,"  which 
sets  forth  methods  of  determining  the  values  of  the  constants 
appearing  in  the  system  equations.  The  system  presentation 
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closes  with  a  "Table  of  Par 'iieters'*  which  lists  all  of  the 
system  variables  and  constaiits. 

The  flow  diagram  which  appears  at  the  end  of  Section  A  is 
provided  principally  as  an  aid  in  the  preparation  of  the 
digital  computer  program  which  solves  the  system  equations. 
This  flow  diagram  could  also  be  used  for  an  analog  solution 
although  other  flow  diagram  arrangements  would  be  more 
efficient  for  that  purpose.  The  following  conventions  apply 
as  to  the  usage  of  the  flow  diagrams:  The  triangles  outside 
the  enclosing  phantom  line  denote  variables  which  are  u^ed 
as  inputs  and  outputs  to  other  systems.  The  numbered 
rectangles  refer  to  equations  within  the  system.  As  an 
example,  in  Figure  5  the  rectangle  numbered  9  indicates 
that  Ter  is  a  function  otMB  and  F'a  and  that  the  equation 
that  gives  the  exact  relationship  is  equation  1.9.  No 
constants  are  shown  in  these  diagrams.  The  triangles 
denoting  integrators  do  not  always  contain  an  equation 
number.  If  the  input  to  an  integrator  is  Xp  and  its  output 
is  xp  ,  then  the  equation  is  implied.  Thus,  as  in  Figure  63, 
if  the  input  to  an  Integrator  is  Ra-  and  the  output  is>^R4, 
then  the  equation  yU/i4.  ^  fRA.  Jt  ,  or  equivantly,  ^^4.  ^  Ra-j 
is  implied.  Because  of  the  size  of  the  six  degree  airplane 
system,  the  flow  diagram  in  Figure  32  is  slightly  different. 
Its  use  is  strictly  limited  to  the  digital  program  generation. 
It  says  that  all  equations  within  one  block  must  be  written 
before  proce  'ling  to  the  next  block.  Thus,  the  first  vari¬ 
ables  to  be  solved  for  are  ^  2’<vv,  yotw  ,•••,  Smu  ,  After 
this  Fw  ,  Fln  ,  *  • ' ,  are  solved  for. 

After  this  ,  Xaau  Fni^  etc. 

The  "Table  of  Parameters"  is  a  listing  of  all  variables  and 
constants  found  in  the  equations  of  that  system.  Each 
variable  is  identified  by  its  symbol,  description,  units, 
and  "Type."  The  "Type"  is  listed  as  v,  v(i),  and  v(o) 
depending  on  whether  the  variable  is  only  used  within  the 
system,  is  received  as  an  input  from  another  system,  or  is 
an  output  to  another  system.  Each  constant  is  identified 
by  its  symbol,  units,  description,  "type,"  and  value.  The 
"type"  for  each  constant  is  always  "c"  and  its  value  is 
that  used  with  the  F-111  antiskid  system. 

Table  1  lists  the  mathematical  conventions  utilized 
throughout  this  study. 
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Table  1  Explanation  of  Mathematical  Convention 


Convention 

Description 

A  dot  over  a  variable  denotes 
differentiation  with  respect  to  time 

Computer 

Notation 

All  variables  are  expressed  in  a 
form  to  harmonize  with  Fortran 
chax'acter  utilization.  Thus  a 
variable  would  appear  as  WTE 

Also,  in  general,  the  following 
practice  is  adhered  to.  If  Xtt 
is  a  variable,  then  XTT  is  its  For¬ 
tran  form.  The  symbol  for  ^tt  is 
XTTD.  The  symbol  for  Xrr  is  XTTDD. 
The  initial  condition  is  denoted  by 
adding  0  (zero) .  Thus  Xtt  at  time  - 
0  is  denoted  by  XTTDO  . 

ZSD<^> 

The  brackets  "<^  "  are  used  exclu¬ 
sively  to  denote  the  position  of  a 
function  argument.  The  script  x 
is  used  CO  denote  an  arbitrary 
variable.  The  parentheses  "(  )" 

are  normally  used  to  denote  multi¬ 
plication. 

Parameter 

Type 

Within  each  table  of  parameters  is  a 
column  which  lists  the  parameter 
"type." 

V  a  variable 

C  a  constant 

vCo)  a  variabl€:  used  as  output  to 
another  system. 

vCt)  a  variable  received  as  an  input 
from  another  system. 
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Table  1  Explanation  of  Mathematical  Convention 


Convention 


Description 

For  symbols  appearing  in  equations 
the  following  conventions  are  used. 


I  »  Capital  "i" 

1  -  One 

&  -  Capital  "Oh" 

0  Zero 
^  “  Capital  "zee" 

2  “  Two 

0  ■  Greek  letter  "Theta"  but  is 

treated  in  Fortran  as  capital©^. 


o 


Oit 


Placing  a  parameter  symbol  between  two 
vertical  bars  denotes  the  absolute  value 
of  the  paramet'*.r.  The  absolute  value  of 
a  signed  number  N  is  defined  as  N  when  N 
is  positive  and  as  -N  when  N  is  negative 
For  example:  J3|  =  3  and  | -3|  ■  3. 

X  The  braces  preceded  by  "MIN"  or  "MAX" 

^  denote  the  value  of  the  least  (or 
largest)  of  the  constant  or  the  para- 
meters  enclosed  within  the  braces. 


1. 


BRAKE  SYSTEM 


The  convenClonal  airplane  brake  consists  of  a  series  of 
discs  vhlch  are  alternately  stators  and  rotors.  The  stators 
are  restrained  from  rotating  about  the  axle  by  splines  or 
keyways.  The  rotors  are  similarly  connected  to  the  wheel 
and  hence  rotate  with  the  wheel  and  tire.  The  brake  torque 
Is  produced  by  axially  compressing  Che  disc  stack;  usually 
by  hydraulically  actuated  pistons.  Many  brakes  use  return 
springs  to  release  the  brake  stack  against  the  return  pres¬ 
sure  of  the  hydraulic  system. 

A.  Mathematical  Description 

In  this  analysis  Xpwlll  denote  the  brake  piston  linear 
displacement.  The  pistons,  rotors,  and  stators  are  treated 
as  a  single  mass  system  In  the  axial  mode  ( Xp  direction). 
The  forces  acting  on  Che  brake  mass  In  the  axial  mode  are: 

a.  Brake  actuation  force:  equals (brake  pressure)  x (piston 

area) 

b.  Force  due  Co  axial  restraint 

c.  Keyway  friction  force 

d.  Brake  piston  seal  friction  force 

e.  Brake  return  spring  force 

f.  Brake  piston  bottoming  force 

Figure  3  shows  the  brake  system  and  the  forces  acting  In 
the  axial  mode.  Each  of  the  axial  forces  Is  established  as 
follows: 

a.  Brake  Actuation  Force 


The  brake  actuation  pressure  Pe  Is  received  as  an  Input 
from  the  hydraulic  system.  The  brake  actuation  force  Is 
given  by  Ps  .  where  is  the  total  brake  piston  area. 

b.  Force  due  to  Axial  Restraint 


The  axial  restraining  force  reflects  the  elasticity  In  the 
brake  discs,  the  back  plate,  and  the  piston  housing  and  Is 
a  function  of  their  cumulative  displacements.  A  way  to 
derive  this  characteristic  is  from  a  curve  of  brake  volumetric 
displacement  vs.  brake  pressure.  This  characteristic  does 
not  Include  friction  or  return  spring  effects. 
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Let  Fp  denote  the  force  due  to  axial  restraint.  And  be 


defined  by 

(1,1)  Fb  ~  Fb2 

(1,2)  =  f  Cfi/ 

1 1  ^  p  ^  S31 

1  o 

1  ^  X  p  <  Sg, 

(1,3)  ffie  - “  Sgj)  -DBe-^p 

Xp  ^  S02 

\o 

I'i  Xp  ^ 

c.  Keyway  Friction  Force 


Let  the  keyway  friction  characteristic  be  defined  by  a  function, 
Gp  ,  where : 
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Figure  4  shows  Gf  as  a  function  of  Xp 


Figure  4  Keyway  Friction  Characteristic 


The  brake  torque,  Tsr  , is  tranferred  to  the  wheel  and  tire 
through  the  rotor  keyways.  Torque,  Tar  ,  is  also 
transmitted  to  the  axle.  The  major  portion  is  transmitted 
through  the  stator  keyways.  The  remaining  portion  of  the 
torque  is  transmitted  as  piston  side  loading  which  results 
from  friction  between  the  pistons  and  the  pressure  plate. 
Let  100  He/  denote  the  percentage  of  brake  torque  transfer¬ 
red  through  the  stator  keyways  and  let  100  Her  denote  the 
percentage  of  torque  transferred  through  the  pistons. 
Naturally,  He.  +  Hei-/.  The  normal  force  on  the 

stator  keys  is  thus  HailTeriyRet  ^  while  the  normal  force 
on  the  rotor  keys  is  (Ter'i  /  Rea .  The  total  keyway  friction 
force  is  then  given  by 


(1.5)  Fkp  iTerj  CJf (Hsty/Rsi  +  l/Rea) 
d.  Brake  Piston  Seal  Force 

Let  Pod  denote  the  seal  frictiai  force.  Then 

(1.6) 

For  =  Gf  (  Hofc  Hgfp  Pe  +  fXstjj^Kp  Hez 
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The  piston  return  force  Fa«  is  given  by 
(1.7)  heft  =  Far©  Csr  Xp 


Brake  Piston  Bottoming  Force 


In  the  brake  released  condition,  an  axial  force  is 
developed  between  the  pistons  and  housing  to  balance 
return  spring  preload.  This  piston  bottoming  force 
is  defined  as: 

(1.8)  f'C6B(Xe-SBe)~De8  X?  FoR 

L  O  FOR  XP^SBd 

This  concludes  the  discussion  of  the  axial  brake  forces. 


Let  Rnr  be  the  number  of  rotors.  Let  VV(3  be  the 
relative  angular  velocity  between  the  rotors  and  stators 
as  received  from  the  wheel  and  tire  system.  The  brake 
torque  Tbt  is  then  given  by 

(1.9)  Tsr  =  Fj9  Rbt^^ 

Where  is: 


(F  Xe'^o 

-G 

/F  Ve  <  O 


Where  Vfe  is: 

(1.11)  Var  Rar  We 


Summing  the  forces  in  the-  axial  direction  yields: 

( 1 . 12)  V/aa  XV  =  p0  ~  Fe  -  h<p-  -  /"©«  -  fe/z  ■*-  Fed 

In  Equation  (1.12)  Wsr  is  the  brake  mass  which  experiences 
axial  motion.  Generally,  is  the  brake  heat  sink  mass. 
Figure  5  shows  the  relationship  of  the  brake  system 
equations.  Table  2  lists  the  system  parameters. 
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Brake  System  Equation  Flow  Diagram 


B«  Parameter  Evaluation 


Figure  6  shows  a  plot  of  brake  piston  displacement  as  a 
function  of  brake  application  pressure  for  a  new  brake. 
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■■■  If  itmiMnaMiiiiiMiMliiMaii 
aiiMimiaMBiuniiiiBiinniianiaiiiBiMimii 

■■iniiiiiiMNiiimiiiiiiiiBaiiiiiiBiminK 

•■iiniMlianillllitiiiliBIliiliiiiiil 


{*  $  yo  \z  )4>  le  zo 


Brake  Application  Pressure  (p«'0 
Figure  6  Brake  Pressure  Volume  Characteristic 

Assuming  that  no  frictional  effects  are 

present,  Cor  and  Cbi  can  be  derived  as  follows  :since  the  initial 
slope  is  due  to  spring  return  force  only,  then 

(1.13)  Cbr  =^^jAap  =^eoj(l3.3)^  =  8850  Ib/in 


From  the  other  slope  on  the  curve, 
(1.14)  Cb, 


-  9950  =  lo.2C^\0  lU/in 


For  a  new  brake  C  *  0  • 
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Assuming  that  the  discs  all  move  together,  since  the  heat 
sink  weight  is  138  LBM,  then  Wee  *  138/386  -  .358  LBF 
SEC^/IN.  The  natural  frequency  is  then  UJ^  «  f 

Otl  tUn'  =  (5/5  Rf^o/S^c. 

Assuming  that  ^  “  .01  (see  page  117), 

(1.15)  081*=  -  ¥',7/ JbF  sec //f? 

U/ri  C/J/S') 

It  is  assuried  that  Xp  *  0  when  the  brake  pressure  is  100 
psi.  Thus 

(1.16)  FtP9  =  flepPs=  (/3.3)00<^=  mo /bF 

3 

Since  the  brake  piston  displacement  is  1.55  IN  before  the 
brake  discs  come  into  contact,  then  Set  ■  1.55/13.3  ■  .1165  in. 

Since  the  F-111  brake  has  8  stators  with  14  rubbing  s*^r£aces, 
Hei  cannot  be  greater  than  1/14.  A  conservatively  high 
value  of  Hdi  05  has  been  assumed  and  it  follows  that 
Hez  =  .95. 

The  brake  piston  seals  are  equivalent  to  MS28775-219.  The 
seal  friction  force  is  established  using  the  procedures 
described  in  Reference  4.  The  seal  sliding  friction  force 
is  a  function  of  rubber  compound  hardness,  amount  of  in¬ 
stalled  compression,  length  of  rubbing  surface,  seal  groove 
projected  area  and  applied  hydraulic  pressure.  For  the 
MS28775-219  size  seal  having  10  percent  installed  compression 
and  70  degree  Shore  A  hardness  the  sliding  friction  force  is 
2.88  Ibf  plus  0.02  Ibf  per  psi  applied  pressure  per  seal. 

There  are  10  pistons  in  the  brake  housing;  therefore, 

(1.17)  Hofc  =  (io)LZ.83)=  ze.e  IbF 

(1.18)  HoFP  =  0o)(.0.O2)=r  0.20  Ibf/ps/ 

Conservatively  high  values  for  the  friction  coefficients  .4^ 
are  estimated  as  Mk  ■  J-5  a'a.d/T.Ke  ■  .10.  is 

estimated  to  be  1.50. 

Values  for  the  following  brake  dimensional  characteristics 
are  then  from  the  appropriate  brake  component  drawings: 

“  4.40  IN,  Rg^  “  6.25  IN,  and  Rgp  *  8.25  IN. 


Observations  of  braking  stops  indicate  that  for 
average  F-111  brake  lining, 
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2.  HYDRAULIC  SYSTEM 


The  hydraulic  system  is  the  brake  actuation  power  source 
and  is  made  up  of  the  four  components  as  shown  in 
Figure  7  :  the  pilot's  metering  valve,  the  antiskid  con¬ 

trol  valve,  the  control  line,  and  the  brake  piston  housing. 
The  pilot's  metering  valve  is  a  pressure  regulator,  usually 
having  a  mechanical  input,  which  has  a  steady  state  output 
pressure  (Pmv)  at  a  level  commanded  by  the  pilot  (Pcom) . 

The  antiskid  valve  is  a  pressure  regulator  which  has  a 
steady  state  output  as  dictated  by  the  antiskid  control 
device.  For  a  modulated  antiskid  system,  the  control  valve 
is  a  variable  pressure  servo  type  regulator  and  for  an 
ON-OFF  antiskid  system  the  control  valve  is  an  ON-OFF  valve. 
The  control  line  is  simply  the  fluid  transmission  line  or 
containment  vessel  connecting  the  control  valve  to  the 
brake  housing.  The  brake  ht-asing  is  a  collection  of 
cylinders  and  pistons  which  act  to  compress  the  brake  discs. 
From  a  hydraulic  system  aspect,  the  control  valve  is  a 
variable  area  orifice,  where  the  orifice  area  is  a  function 
of  spool  position.  The  control  valve  spool  position  is 
received  as  an  input  from  cOTiputations  described  in  a 
section  devoted  to  Che  operation  of  the  control  valve. 

In  the  description  of  the  brake  actuation  system,  there  are 
two  principal  effects  which  should  be  accounted  for.  The 
first  is  the  time  lag  which  exists  between  the  control  valve 
output  pressure  (Pev)  and  the  actual  brake  pressure  (Pb). 
This  lag  is  caused  by  the  fluid's  resistance  to  flow  due  to 
inertia  and  friction  and  by  the  brake  pressure's  dependence 
upon  fluid  volume  within  the  pressure  cavity.  The  second 
effect  is  the  instantaneous  brake  pressure  intensity  as 
influenced  by  fluid  inertia  and  the  combined  elasticity  of 
the  fluid  and  the  pressure  cavity.  Rapid  valve  operation 
can  cause  pressure  overshoot  and  oscillation  du?  to 
"water  hammer"  effects.  This  overshoot  can  cause  excessive 
brake  torque  and  may  interfere  with  proper  control  valve 
operation.  The  pilot's  metering  valve  pressure  drop  and 
response  characteristics  are  included  in  the  actuating 
system  description  so  that  these  effects  upon  antiskid 
operation  can  be  examined.  To  allow  for  a  variety  of  brake 
actuation  systems  which  might  be  encountered,  provision  is 
made  to  accommodate  both  hydraulic  and  pneumatic  actuation 
media.  The  line  connecting  the  control  valve  and  the  brake 
can  be  treated  as  a  separate  fluid  cavity  or  the  effects  of 
its  volume  may  be  lumped  with  the  brake  as  would  be  appro¬ 
priate  for  a  short  line. 
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A.  Mathematical  Description 


Figure  8  is  a  schematic  of  the  brake  hydraulic  system. 

The  analytical  procedures  of  References  5  and  6  are 
utilized  to  mathematically  describe  the  system. 

Let  PcoM  denote  the  brake  pressure  which  is  commanded  by 
the  pilot  and  define  PcoM  increases  from  a 

minimum  value,  Pr  ,  (reservoir  pressure)  to  the  desired  steady 
state  value  P^p,  as  a  linear  function  of  time  over  an 
interval,  T^^p,  as  follows: 


(2.1) 

PcofA  =  T  C  PcP  -  Pfi)/ TeP 

Pcp 

+  IF  c>  -  Tep 

IF  Tcp  ^  T' 

The  metering  valve  attempts  to  maintain  Pj^  at  the  level 
of  P^oi^.  The  metering  valve  spool  displacement  is 

defined  by  equations  (2.2)  and  (2.3). 

(2.2) 

Vmv  = 

/Til/  f  Rccm  ~  Rri !/  } 

(2.3) 

%v  " 

J  yW/V  J"  Oj  VfAv} 

<  Vmv 

/  fYiAK  1  Vmv} 

1  f  S/v\Vu  ^  Xmi/ 

Jp  S/y\tfl,  ^  ^  rn\/U 

XmV  S rriVi. 

Let  be  a  function  defined  as  follows: 


(a)  For  hydraulic  fluid 

(2.4)  <  X  -  5-/6W  (X-  Y)  ij  K-  Y( 


(b) 


For  compressible  pneumatic  fluids 


(2.5) 


IF  X>Y  X  ^  Y/Rcrit  Rc/nr  - 

=  k[i-  Vr  i 

IF  anci  x^r/R.  cRvr 


x[''  ( 


IF  Y^X  3i^d  Y-  ^/RcR\T 

f<.x^Y}  =  -<^<Y^xy 
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Let  A[vjv('ic>be  defined  by: 

(2.6) 'A|viV^'X>  -  f  Am  Vt'  ^  VO 

'  JY12tx{  AmvL,  ^  A  (Villi'/  S^iy ,5  ^  ti  X<-C>MV0 
Let  and  /V,yji  be  defined  by: 

(2.7)  Am  Vo  “  Aivjv  <  X  MV  ^ 

(2.8)  Amvi2.  "  Amv^“ 

Then 

(2.9)  (ji''s  =  A,v,vs  5^  ^  ^  >  Ff^iv^ 

(2.10) <+i?_  -  Aivi  vE  ^  A  j  ^ 

Let  Vmvv  be  the  fluid  volume  from  the  output  of  the  meter¬ 
ing  valve  up  to  the  input  of  the  control  valve. 

Then 

(2.11)  P|viv  "  (  BMv/VMVvX(?i  ~Qt^)y  ■bQcv'l  ) 

Let  be  defined  by; 

(2.12)  *  /  Acvo  t-f  ^^‘-'cvc 

'  max  it  Acvu,  i.  Acv'o/Scvoi  <-Smvo 

Let  Acvs  ^nd  Acv^  be  defined  by 

(2.13)  Acvs  =  Acv<>^cv“Sci.> 

(2.14)  Acvt^  -  Ac  V  ^  "  2cl  "  ^cv  ^ 

Then 

(2.15)  G?niv  ~  Acvi  ^  ^  FWiv  y  Pcv  y 

(2.16)  C?c VR  “  Acve  ^  ^  Pcv  >  P:v  12.  ^ 

(2.17)  ^cvR  --  (B  C  VR  /  Vcve  )(Qcvi2  QrC  "^QcvJ^ 

(2.18)  Cec  "  A^c^c  <Pcv»^,  PR> 
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Tae  volume  of  the  cavity  occupied  by  the  brake  actuation 
media  is  established  by  equation  (2,19)  as  follows: 

(2.19)  Ve  ^  Vflo  + Xp 

Three  options  for  the  control  line  mathematical  description 
are  provided  to  cover  a  variety  of  circumstances  which  may 
be  encountered.  The  third  option  is  representative  of  a 
typical  aircraft  installation  and  is  used  in  analyzing  the 
F-111  system. 

The  first  option  is  for  a  control  line  with  hydraulic  fluid 
considering  volume  effects  only.  This  option  will  not  pre- 


diet  'water  hammer'  but  is  satisfactory  for  many  cases, 
particularly  for  the  case  of  a  short  control  line  50  inches 
or  less  in  length.  The  following  equations  describe  the  . 
first  option: 

(2.20a) 

Q,CV  ~  QmV  -  +  Cjcvt 

(2.21a) 

Pev  *  C  6©/Va)  (  Qcv  -  //bps  ^  p) 

(2.22a) 

Pei  =  Pev 

(2.23a) 

P0  r  Pbx 

(2.2^a) 

Qg  =■  Qcv 

The  following  equations  are  applicable  to  the  second 
option  for  the  control  line  using  compressible  pneumatic 
fluid. 

(2.20b) 

(2.21b) 

Pev  =  -  Pev  /idPsKp/Se) 

(2.22b) 

Par  3  Pev 

(2.23b) 

Pe  ~  Pei 

(2.24b) 

C^e  =  Qcv 
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The  third  option  is  for  a  control  line  with  hydraulic  fluid 
where  both  volume  and  inertial  effects  are  considered  and 
is  described  by  the  following  equations: 

(2,20c)  C^cv  -  (AsL/RHoSBL)^R:t>'-fhi-DR.8LC^c\y‘’D'^6LQc\/\(^c\^ 

(2.21c)  Pci^  ~ 

(2.22c)  4r  -  /f^^X Qa/ 

(2.23c)  Qs  =  /!bo  P^<Pe2^  Ps} 

(2.24c)  Xa  ^  (Ss/vXiQB-/?»^skp) 

In  this  study  the  brake  system  hydraulic  supply  pressure, 
Pg,  is  treated  as  a  constant.  If  Pg  varies  significantly 
due  to  operation  of  other  aircraft  hydraulic  system  equip¬ 
ment,  this  variable  pressure  defined  as  a  function  of  time 
may  be  used. 

B.  Parameter  Evaluation 

For  this  study  the  third  optional  control  line  description 
as  applied  to  the  F-111  is  of  primary  interest.  For  this 
case  MIL-H-5606  hydraulic  fluid  is  used.  The  hydraulic 
fluid  properties  for  a  mean  temperature  of  100®F  and  1500 
psi  are: 

(1)  Adiabatic  bulk  modulus:  8-  248,000  psi 

(2)  Density:  /?wo^.781  x  10“^  LBF  SEC^/lN^ 

(3)  Kinematic  viscosity:  ,0267  IN^/SEC 

The  system  supply  pressure  is  3000  psi  and  the  return  pres¬ 
sure  is  100  psi.  Initially,  all  flows  are  zero  and  all 
pressures  except  the  supply  pressure  are  at  100  psi.  The 
pilot* s  input  command  pressure  Pq()[^  is  also  100  psi.  The 
pilot's  input  PcOM  go  from  100  to  1500  psi  in  0.2 

seconds.  Thus  T^^p  =  0.2  sec  and  P^p  =  1500  psi. 

Metering  Valve 

When  the  metering  valve  spool  is  centered,  the  flow  area 
is  essentially  zero  for  both  the  return  and  supply  lines. 

In  this  spool  position  Xmv  =  0.0.  From  equation  (2.3)  the 
spool  is  constrained  to  stay  between  Smvl  and  Srvwu. 
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For  the  meterlM  valve,  ■  -  .06  In  and^m^i;-  ,06  In. 

However,  when  Am^  la  at  +  .05,  the  valve  area  has  reached 
Its  maximum  for  the  flow  .  When  Xiki/  -  -  .05,  the  area  Is 
maximum  for  the  return  £1ow(2r*  Thus  ■  .05.  By  actual 
mea^iirement,  with  the  valve  full  open  (area  *  )  at  100^ 

F,  th'j  flow  Is  9.23  In^/sec.  at  200  pslAP.  Thus  from  (2.9) 

-  (2.10), 

(2.25)  =  a/'flf =  =,  6i3  ,;^y^ec-)L  lU) 

In  the  F-lIl  system,  the  metering  valve  is  situated  next  to 
the  control  valve  so  that  the  volume  V^vv  is  quite  small. 

was  calculated  from  the  valve  drawing  as  being  about 
I.O  in^.  Also,  the  valve  body  is  considered  to  be  much 
stlffer  than  the  hydraulic  fluid  so  that  the  effective  bulk 
modulus  is  the  fluid  modulus.  Thus,  Bfr-  »  248,000  psi. 
was  estimated  from  analog  studies  to  be  about  .05. 

Control  Valve 

For  the  control  valve,  Xcv  -  0.0  when  the  spool  is  centered. 
At  this  point  the  flow  area  is  zero  so  that/)c/<.-  0.0.  The 
flow  area  remains  zero  for  -  .005  s  Xcv's  ,00'.  Thus  the 
valve  has  an  overlap  of  .005  in.  and  Stx  “  .005.  An  addi¬ 
tional  movement  of  .030  in.  produces  full  area  so  -  .030. 
By  actual  measurement  at  this  position  at  100°F,  the  flow 
is  7.7  in^/sec.  at  50  psi  A  P.  Thus 

(2.26)  AiVo^Sl/'fiF^7.?/>fFb-  hO‘10 

The  following  values  are  estimates  of  the  return  character¬ 
istics  of  the  control  valve:  ■  2,0  in^,  -  248,000 

psi,  /?RC-  l.C  inV(sec)(lbf)l/2. 

Control  Line 

The  control  line  is  1/4  inch  outside  diameter  steel  tubing 
having  0.14  inch  wall  thickness  and  internal  cross  sectional 
area, equal  to  .0386  in^.  Because  of  the  thin  wall, 
the  tube  elasticity  greatly  reduces  the  bulk  modulus.  The 
equivalent  bulk  modulus, may  be  calculated  from 


Where  B  -  Fluid  bulk  modulus 

E  ■  Young's  modulus  of  tube  material 
D  «  Mean  tube  diameter 


t  »  Tube  wall  thickness 


Thus 

(2.28)  Bbl 


2i-QOOO _ 


2/7,  700  2s/ 


The  control/line  length,  is  191  inches  with  various 

types  of  flow  restrictors  according  to  the  following  table. 


Table  3  Control  Line  Restrictions 


Description 

"K"  Value* 

Number  n 

nk 

An815-4J  Union 

.54 

1 

.54 

AN832-4J  Union 

.54 

1 

.54 

AN821-4J  Elbow  (90°) 

1.23 

4 

4.92 

AN837-4J  Elbow  (45°) 

.89 

1 

.89 

90°  Tube  Bend 

.01 

12 

.12 

90*^  Hose  Fitting 

1.25 

1 

1.25 

Total 

8.26 

2 

«  KV  /2g  Where  V  is  the  velocity  in  the  line. 


The  "K”  values  in  Table  3  were  derived  from  information 
contained  in  Reference 

Equation  (2.20c)  is  the  result  of  summing  forces  on  the 
mass  of  fluid  in  the  control  line.  The  friction  losses 
are  depicted  by  a  turbulent  flow  loss  Dtbl  and  a  lami¬ 

nar  flow  loss  Drbl  Qcv  It  is  assumed  that  all  the  turbu¬ 
lent  flow  losses  come  from  elbows,  etc.,  which  are  listed 
in  Table  3.  The  loss  due  to  the  line  itself  is  considered 
to  be  always  laminar.  This  assumption  of  laminar  flow  for 
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the  line  is  justified  for  two  reasons:  (1)  the  loss  in 
the  line  is  small  compared  to  other  losses  in  the  system; 
(2)  the  flow  is  normally  laminar  anyway  (Re5molds  humber 
is  less  than  6000  for  the  F-111  system). 

For  the  turbulent  losses 

(2.29)  AP  “ 

-  Kp  1/^-2 

=  q!- 

Thus 

(2.30)  Ptbl  ^  K  P 

1 

Z.  i.  OSH)^ 

=  ^Z/io  /hi- 

For  laminar  losses,  at  temperatures  normally  encountered, 
the  ''oscillatory"  friction  is  higher  than  the  steady  state 
friction.  See  Reference  9.  The  pressure  loss  can  be 
written  as 

(2.31)  AP  =  Ru  ^ 

For  the  steady  state  case  as  shown  in  Reference  6, 

(2.32)  R^~8rffzy- 

In  Figure  ]0  values  for  this  theoretical  steady  state 
are  compared  over  a  range  of  temperatures  to  values  from 
Reference  9  which  were  experimentally  established  for 
oscillatory  flow.  Since  the  hydraulic  flow  in  the  brake 
control  line  associated  with  antiskid  operation  is  transi¬ 
tory,  the  laminar  flow  resistance  base  on  experimental 
measurements  for  oscillatory  flow  is  used. 
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Figure  10  Hydraulic  Fluid  Damping  Characteristic 

From  Figure  10  at  100°F  ^Lfor  the  experimental  oscillatory 
case  is  1.5  X  10*4.  LBF  SEC/IN^ 

Therefore: 

(2.33)  PffBi  °  (RLjjSed,  =  (/Sx/o'^Kifl) 

(^61-f  (.03Q^Y 

=  .19.  !6fsecy%^ 

When  a  "lumped  parameter"  type  analysis  as  described  by 
equations  (2.20c),  (2.21c)  and  (2.22c)  is  used  for  the  con¬ 
trol  line  the  resulting  natural  frequency  is  somewhat  lower 
than  the  actual  line,  if  the  actual  line  volume,  Vbl  ,  is 
used.  The  value  of  is  adjusted  as  follows  to  achieve 
the  correct  natural  frequency  for  the  "lumped  parameter" 
description. 
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Consider  hydraulic  fluid  flowing  through  a  line  with  cross 
sectional  area,  and  divided  into  segments  having  equal 
length,  >5*^  as  shown  below. 


If  each  segment  is  treated  as  a  separate  pressure  vessel 
having  volume,  with  a  flow  in  and  a  flow  out,  and  if 
equations  of  the  form  of  (2.20c)  (2.21c)  and  (2.22c)  are 
written  for  these  pressure  vessels,  neglecting  friction, 
the  following  expressions  are  obtained: 

(2.34)  Qz  -(/!//=  S)iP,-Pz) 

(2.35)  P,  =  CB/v)LQ^-aCj 

(2.36) 

By  substituting  equations  (2.35)  and  (2.36)  into  equation 
(2,34)  differentiat^^d  once  with  respect  to  time  the  fol¬ 
lowing  differential  equation  is  formed: 

(2.37)  Qz  - 
or 

(2.38)  2  {ft  a/f>S\/)  Q,.  --  i  fie/fSV)(Q,*Cli) 

Equation  (2.38)  establishes  that  the  natural  frequency  of 
each  line  segment  is: 

(2.39)  = 

zrr  \ 

However,  vibration  theory  considering  distributed  mass  and 
elasticity  establishes  the  speed  of  sound,  C,  in  the  line 
as : 

(2.40)  C  - 


42 


For  fundamental  mode  oscillation  in  a  closed  end  tube  having 
length,  5,  the  natural  period ,Tc,  is: 


(2.41)  Tc  =  2yc 

Therefore,  the  natural  frequency,  ,  of  an  actual  tube 
segment  is: 

(2.42)  =  l/Tc.  ^  Cl/2S)  'f~B/p  CPS 

By  equating  the  two  expressions  for  natural  frequency, 
equations  (2.39)  and  (2.42),  the  volume  of  the  line  segment 
which  will  have  the  same  natural  frequency  as  the  actual 
is  established  as; 

(2.43)  y=  2fiS'/7r^ 

Thus, 


(2.44)  Vbl--  A^  J^blSbl.  = 


Brake  Housini 


o 

The  brake  housing  has  ten  pistons  of  1.33  in^  area  each. 
Since  the  number  of  pistons  serviced  by  one  control  line 
is  five,  then  f\^?s  ^  5(1.33)  =*  6.65  in^. 

The  fluid  volume  in  the  brake  housing  with  the  pistons 
bottomed  (^p-d)  is  8.00  in^.  Thus  Viso  «  4.00  in^  or  one 
half  the  total  volume.  The  orifice  coefficient  was 
estimated  to  be  about  2.0 


Operational  Systems 


The  option  1  system  neglects  the  line  inertial  effects.  The 
parameters  have  the  same  value  as  the  corresponding  parameters 
for  the  option  3  system,  escept  that  \Jeo  should  include  any 
line  volume.  Thus,  for  the  F-111  system,  with  the  option  1 
system,  Vio  =  4',oo  +  .osec,  L^9f)  =  //.S4> 


The  option  2  description  is  used  for  systems  with  compressible 
pneumatic  fluid.  The  appropriate  parameters  will  be  evaluated 
for  nitrogen  at  100°F  as  the  fluid  media  and  isothermal 
processes  are  assumed  except  for  orifice  flow  calculations. 
While  the  heat  transfer  characteristics  of  the  brake 
system  components  have  not  been  rigorously  evaluated,  the 
usual  component  installation  is  such  that  assuming  isothermal 
processes  is  valid.  The  mathematical  description  of  the  brake 
actuation  control  system  using  compressible  pneumatic  fluid 
is  written  using  equations  of  the  same  general  form  as  for 
those  describing  the  hydraulic  system,  thereby  minimizing  the 


the  number  of  equations  and  enhancing  computation  flexibility. 
Utilizing  the  hydraulic  equations  when  pneumatic  fluid  is  used 
requires  that  the  appropriate  parameters  be  expressed  in  suit¬ 
able  mathematically  equivalent  terms.  Consider  the  character¬ 
istic  equation  of  state  for  a  perfect  gas: 


(2.45)  P  mRT 

V 


And  the  definition: 


(2.46) 


dp  _  ^  dm 

dt  dt 


jv  IP  ^ 
dV  If  dV  If 


For  the  assumed  isothermal  process,  substitution  of 
equation  (2.45)  into  equation  (2.46)  gives: 

(2.47)  p  ^ 

For  those  cases,  such  as  for  the  metering  valve  and  control 
valve  pressure  cavities,  where  the  volume  is  not  changing, 

V  is  zero  and  equation  (2.47)  reduces  to: 

(2.48) 

For  hydraulic  fluid,  P  is  described  by  equations  having  the 
form  of  equation  (2,49)  below.  (See  equation  (2,11)  for 
instance.) 

(2.49)  h  (i)Q 

Noting  the  similarity  between  equation  (2.48)  and  equation 

(2.49)  it  is  obvious  that  if  PT  is  used  in  place  of  8  and 
if  m  is  used  in  place  of  CR  ,  the  "Hydraulic"  equations  can 
be  used  for  computing  performance  of  a  system  using  pneu¬ 
matic  fluid.  Thus,  Be  »  =  Brr^v  *  R7",  ^ 

For  nitrogen  R  *  662.4  in  and  at  100  F 

RT  -  (662.4)  (460  +  100)  =  371  x  10^  ,»  \hf /ibm . 


Since  P/RT  =  M/V,  equation  (2.47)  can  be  written  as 
(2.50, 

Equation  (2,21b)  is  obtained  by  substituting  8b  for  RT  i 
f\BPs  X?  for  V'  ,  and  ^5  for  in  equation  (2.50),  thereby 
accounting  for  the  change  in  brake  volume  caused  by  piston 
movement. 


Equation  (2.51)  below,  from  Reference  6,  describes  the  mass 
flow  rate  of  a  gas  from  a  container  having  high  pressure, 

Ph  ,  through  an  orifice  of  area,/?o,  to  a  container  having 


44 


low  pressure,  R.  . 


(2.51) 


Equation  (2.52)  below,  from  Reference  6,  describes  the 
volximetric  flow  rate  of  hydraulic  fluid  through  an  orifice 
under  similar  circumstances. 

(2.52)  <5= 

Both  equations  (2.51)  and  (2.52)  can  be  written  in  the  form 
Gl  =  Pri,  Pu)  where  is  a  flow  function  as 

defined  by  equations  (2.4)  and  (2.5)  for  the  appropriate 
circumstances  and  where /9f is  a  flow  coefficient  accounting 
for  orifice  and  fluid  properties.  For  the  case  of  hydraulic 
fluids  a  value  of  C  o  ~ /OSS  has  been 

established  by  experience  as  being  representative  of  an 
average  orifice  (i.e.,  Cx>^  0.65).  The  metering  valve  flow 
coefficient,  Wmvo  ,  previously  computed  is  0.653  IN^/sec 
therefore,  the  apparent  actual  orifice  area,  /?©  ,  for  the 
metering  valve  is  -O.  6S3//oSf  = 

For  the  case  of  the  pneumatic  system  with  nitrogen  at  lOO^F 
as  the  working  fluid  and  using  Cp  =  2300  in  Ibf/lbm?  F,  and 

k  -  662.4  in  Ibf/lbm  °F; 

(2.53)  =  Co/^o 

.-e  '  r  “ 

f — 3^5 

=  O.  4-3  '*-\o  ^  Ibrv,  in  lb  f  jec 
Using  the  same  procedure  establishes  that: 

/ic^o  =  O,  '7/^  ^\o'^  Ibfyy  Ibf  see 
bi  f^C  -  0.^S8>!ID^  Ibm  Ibf  Jec. 
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Calculate  from  ,  see  Equation  (2.5) 


Table  4  (Contd) 


Option  3  Brake  actuation  system  using  hydraulic  flu.d  where  the  cont 
line  description  considers  both  volume  and  inertia  effects. 


3a  AIRPLANE  SYSTEM  (FLYWHEEL) 


Figure  11  shows  the  model  for  the  airplane  system  as  It 
might  be  simulated  with  a  dynamometer  flywheel-  set-up.  The 
mass  W/i  Is  supported  by  the  tire  and  Is  determined  by  the 
percentage  of  the  airplane  weight  carried  on  one  main  gear. 
The  mass  represents  some  part  of  the  airplane  structure 
whlcn  could  vibrate  In  sympathy  with  certain  ground  discon¬ 
tinuities  such  as  wing  mounted  fuel  tanks  or  armament.  The 
foices  Flo  and  Fal  act  on  Wa  because  of  gravity  and 
aerodynamic  lift,  respectively. 

A.  Mathematical  Description 

The  shock  strut  stroke  Is  denoted  by  . 

This  stroke  Is  determed  by  Z  and  Zwm* 

(3a.  1)  Esm  = 

(3a. 2)  EsAf 

The  shock  strut  force  Fvai  Is  given  by  equation  (3a. 3) 

(3a. 3) 

Let  Zgo  and  denote  the  height  and  slope  of  the  ground 
(or  flywheel  surface).  Let  5a7  denote  the  tire  deflection. 
Then  S/r\  and  are  determined  by 

(3a. 4)  S/v)  =  0.0  j  CXf^'Zw/v*} 

(3a. 5)  SfA  =  ^  Vt  -  Z'zzm 

The  force  Fvm  acting  vertically  upward  on  the  tire  Is  then 
given  by 

(3a. 6)  Fa/m  =  S/^  (Cmt  ^ 
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Summing  forces  in  the  vertical  direction  on  the  unsprung 
mass  W^v  ,  there  follows: 

(3a,  7)  Wvv-v  Z  WM  “  ~  ^RV 

Where  F^rv  is  the  tire  unbalance  force. 

For  the  mass  ,  summing  forces  vertically  gives: 

(3a. 8)  WftR 

(3a. 9)  ~  (  2  ~  ^Ap  ^  ^  ~  ^Aje) 

The  aerodynamic  lift  and  drag  forces  and  are 
defined  as  follows: 

(3a. 10)  Fal  “  ^al 

(3a.  11)  Fad  = 

The  equation  which  determines  Z  is  given  as 
(3a.  12)  (  Wa  "  Wae)  2  "  Fvm 


The  equation  for  the  flywheel  velocity  is  given  by 
(3a.  13)  W^r  Vp  =  Frw  "  ^  ^7 

Where  Fvh  is  a  force  equivalent  to  engine  thrust  and  Wat 
is  the  airplane  mass.  The  aircraft's  longitudinal  dis¬ 
placement  is  established  by 

(3a.  14)  Xp  =  {Vp  dt  +  Xpo 


The  equation  flow  diagram  for  the  airplane  system 
(flywheel)  is  shown  on  Figure  12. 


B.  Parameter  Evaluation 


Shock  Strut  Characteristics 


Figures  13  and  14  show  the  main  gear  load  and  damping 
characteristics  for  one  gear. 


gsi| 


4  a  10  )4-  14  18  20  22  24  24  28  30 

Stroke  in) 

Figure  13  Hain  Gear  Damping  Curve 


4  4  8  10  12  14-  )4  )  8  20  22  24  24  28  30 
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Figure  14  Main  Gear  Air  Load  Curve 


Vertical  Tire  Characteristics 

In  equation  (3a. 6)  It  has  been  assumed  that  the  tire  loading 
characteristic  Is  given  by  an  equation  of  the  form 

(3a. 15)  F  -  'S'(c^OS) 

Let  the  following  terms  be  defined  for  a  tire: 

Fn  =  Rated  load 
pR  =  Rated  pressure 
=*  Rated  deflection 

If  P  is  the  actual  pressure,  then  obviously  the  tire  spring 
rate,  C  ,  is 

(3..16) 

From  reference  1  (Equation  132)  the  damping  force,  Fp  , 

Is  established  as: 

(3a.  17)  Fp  C  5 


It  Is  assumed  that  the  damping  force  Is  related  to  the 
undamped  natural  frequency  at  rated  conditions.  The  un- 
dflmp‘>H  ''qtur“T  frequency,  a;  ,  is  established  as: 


(3a. 18) 


a> 


c  fF  ^  Ci 
f  /yi 


J  Sfi  Fa. 


Where  G  -  386  IN/SEC^.  Also  from  Equationsl37  and  138  of 
Reference  1  I 

(3a.  19)  '>1  -  ^  / [ I 

Where 


The  main  landing  gear  shock  strut  linear  damping  coefficient, 
Dv/v\,  is  set  equal  to  zero  for  the  example  problem. 


The  unsprung  mass,VVwv/j  experiencing  vertical  motion  is 
6,44  Ibra.  Thus,  Wws/ (644)/386  «  1.667  Ibf  sec^/in. 


As  previously  assumed  in  Equation  (3a.  15),  ho  = -S' OS' 
Equating  the  two  expressions  for  Fq  at  rated  deflection 

(3a. 20)  07  C  ^  Sr.  ^ 

uu 

Or 

(3a. 21)  D^V£_=  lf\ofl£. 

uj  S/i  ^ 

For  the  47  x  18  -  18  26  ply  rating  F-111  main  tire, 

■p=P^=  150  psi 

Fk  =  38,100  lb.  and  Sk.  =  4.00  IN. 

Thus 

(3a. 22)  CmT  =  [-£-]( =  (yjra){Sg/oo)  _  IU/m 

(ISOM4-.00)  ~  ^ 

(3a. 23)  D„r  -  /_P_\/  TIFA 

(  )(  V  '  ^ 

^  0:££)  {.l)(38l00]  'j/TTo  ^  ^4.Z-f-  MseJw^ 

(/so)  (i-.QSf  !  sei>  ' 

Aircraft  Characteristics 

For  the  example  problem,  an  airplane  weight  of  57,000  lb, 
is  used.  The  static  vertical  load  on  one  main  gear  is 
25,200  lbs.  so  that 

(3a. 24)  =  i>S’.o  lb(-  sec^^iN. 

For  a  velocity  of  \Jf  ^24-00  IN/SEC  and  a  representative 
tire-to-runway  braking  coefficient  of  .45  at  the  main  wheel, 
the  tire  load  is  21,400  lbs.  Thus  Flo  21,400  lb. 

The  total  aircraft  mass  is  \klf^r  -  S7O0O^G  =  14-l.S 


The  mass  V/aa  is  used  to  simulate  some  airplane  resonant 
effect.  For  illustrative  purposes,  it  is  assumed  that 
War  =  1000  LBM  =2.59  LBF  SEC  ^/IN  and  has  a  natural 
frequency  of  12  cps.  Therefore,  since  oj  -  arrC  =  7 s,4-  rad/s<sc 
and  k.  - 

(3a.  25)  Cftij  -  VflR.  =  (75.+)*(2.59)  =  14,720  Ik/in 

Using  3  percent  critical  damping  gives 
(3a.  26)  Pa/2.  ~  (.03^  2  '/Crc.Wae.  — 

=  C.03)  2  '1  (  |4j720)C?r5^)  =  /  1,72  /b  sec/tn 

The  initial  conditions  are  calculated  for  equilibrium. 

At  time  =  0,  let  =  0  so  that 

since  >  is  Slways  0.  Let  =  1200  IN/SEC  and 

assume  that  =  Cad  =  0 

From  equation  (3a.6),j 

(3a. 27)  5m  =  /Cmt  =  2/,400/9S30  ^  2.2^5  im 
From  equation  (3a. 4), 

(3a, 28)  Z-wMo  "  "2.24S’  in. 

From  figure  14,  when  Fy^.,5  *  2l,40oiL.  5 
ZsM  *  23.98  in  and  from  equation  (3a, 1), 

(3a. 29)  Zo  =  ZwMo  -  =  -2. .245  -  23.980  =-26.225  IN. 

Also  Z-flE.0  =  2o  =  -Zlo.ZZS'm. 

For  the  example  problem  the  effects  of  aerodynamic 
forces  are  not  included  in  the  flywheel  simulation; 
therefore,  Cao  =  0.0  andC/»i  =  0.0. 

The  unsprung  mass  moving  vertically,  Wwv  ,  is  che 
same  as  Wgw  described  in  the  Section  4a  Wheel  and 
Tire  System  (Flywheel)  for  horizontal  motion.  There¬ 
fore,  =  1.60  Ibf  sec^/in. 

The  average  en^^ine  idle  thrust  is  1000  Ibf. 

Therefore,  Fth  =  1000  Ibf. 
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SYMBOL  TYPE  '  VALUE  1  UNITS  DESCRIPTION 
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3b 0  AIRPLANE  SYSTEM  (3  DEGREE) 


The  three  degree  airplane  system  is  built  around  a  rigid 
body  airplane  which  is  allowed  to  move  vertically,  hori¬ 
zontally  (parallel  to  the  runway  centerline),  and  rota- 
tionally  in  the  pitch  mode.  This  model  provides  for  the 
interaction  of  the  anti-skid  system  with  those  effects 
which  are  related  to  airplane  pitch.  This  includes  such 
pitch  effects  as  change  in  the  aerodynamic  lift,  drag,  and 
moment  due  to  change  in  wing  angle  of  attack,  change  in 
the  aerodynamic  lift,  drag,  and  moment  due  to  changes  in 
elevator  deflection  as  dictated  by  the  stability  augmenta¬ 
tion  system  (pitch  mode),  change  in  tire  loaning  due  to 
braking  pitch  moment,  and  the  effect  of  ground  slope  and 
roughness  as  reacted  through  both  the  main  and  nose  gears, 

A.  Mathematical  Description 

Figure  15  shows  the  three  coordinates  which  describe 
the  airplane  position  relative  to  reference  points  on  the 
earth* s  surface. 


Figure  15  Airplane  Coordinates 
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Figure  16  shows  the  gear  extended  dimensions  as  mea¬ 
sured  in  the  airplane's  water  line-fuselage  station  refer¬ 
ence  system. 


Figure  16  Airplane  Geometry 


Let  denote  the  runway  profile  height  and  letZsop^'^'^ 

denote  the  runway  profile  slope. 


Nose  Gear 


Let  and  denote  the  nose  strut  ^stroke  and  stroke 
velocity.  From  Figure  17  ,  2.^^  and  are  given  by 


(3b, 1)  'Svw  "  2.  ~  ShKjQ 

(3b, 2)  ~  Z  “  Q 

The  nose  gear  shock  strut  force  is  then  given  by 

(3b. 3)  DmW  I  Zsw  I 

the  normal  ground  force  at  the  nose  gear  is  given  by 

i 

(3b, 4)  "  Skj  T  ^*0  ) 


' '  tr  ^  ^  -7 -.'7.'  ■:* 


V.%'  o 


where  Sw  is  the  nose  tire  deflection.  Su  and  Skj  are  given 
by; 

(3b, 5)  S^j  “  {  0,0  j  ^  ^ 

,  •  ' 

(3b, 6)  Snj  —  ^  ww 

Summing  vertical  forces  on  the  nose  wheel, 

•  4' 

(3b. 7)  Wwu  Lwpo  ~  (Oh)  “  Pvw 
(3b. 8)  Fd>j  ~  ^WK» 

Main  Gear 

Let  Zsm  denote  the  main  gear  stroke  and  stroke 

velocity: 

(3b, 9)  ^  ^  ^riivi  Q 

(3b,  10)  ■  Z-VM  ”2  "h  SmmQ 

The  main  gear  shock  strut  force  is  given  by: 

'  *1*1 
(3b, 11)  '*'Dvm  ^sm 

Let  Sm  denote  the  main  gear  tire  deflection.  Then  the 
tire  normal  force  is  given  by: 

(3bol2)F^u^  -  ?>mCCait  '^  ^ 

(3b, 1:3)  =  V>ia^  i  0.0  ^  2-<bIiXXwM/  -  2  ^jvy  J 

*  ,  • 

(3b, 14)  S,V1  ”  X 

Summing  vertical  forces  on  the  main  wheel, 

(3b, 15)  Ww/M  ^ivrt^  ~  Fmiwi  “  ^vm  ^ 

Figure  18  shows  the  model  of  the  main  gear.  With  the 
assumption  that  the  gear  weight  is  much  less  than  the  air¬ 
plane  weight  (that  is,W,^<<WA)>  it  follows  that; 

(3b, 16)  W(^  Li  "  Fa  Ss u,  “  ( ^<3 k'F  ^CrL  ^  ”  ~^s 
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Figure  18  Main  Strut  Model 
where  i-®  determined  by: 

(3b,  17)  'Z<jL  “  ^SM 

Fp^  can  then  be  computed  from 
(3b,  18)  Fdu.  ^  CFe.Z,6L^’^s 

where 

(3b. 19)  Fa--  SsaCCaCQ  -06)  +  D,,^O-0^)) 

and  are  outputs  from  the  tire  and  wheel  system.  The 
horizontal  axle  reference  location  is  denoted  by  Xax 

is  given  by: 

(3b, 20)  Xax  =  X  -  Sh^v,  b  (  Ssa  "*■  Z,6u  )  Qfr 
(3b,  21)  Xa;.  =  X+  (S,,a+  Ec.L)®fr 
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Thrust 


Referring  to  Figures  16  and  17  ,  if  is  the  thrust, 
then 

(3b, 22)  TH  (‘=^TM  ^ 

(3b,  23)  ^TH  ~  5  th  F7H 

Aerodynamics 

The  dynamic  Air  Force  is  given  by: 

(3b. 24)  Qa  =  X^AeefI^ha  / 288.0 
The  aerodynamic  lift,  drag,  and  moment  are  then  given  by: 
(3b, 25)  “  ^AL^A 

(3b, 26)  Fad  ^  QdQa 
(3b, 27) 

If  denotes  the  wing  angle  of  attack  relative  to  the  air, 
then: 

(3b. 28)  .xw  =  ^0  +  Cl80/Tr)(G5  -£/X) 

Let  ShT  denote  the  horizontal  tail  deflection.  Then  the 
aerodynamic  coefficients  are  given  by: 

(3b, 29)  (^Au  “  5al.  *■  ■^£At5Ht 

(3b, 30)  Cad  -  Cad  ^  ^  Ead^mt 

(3b, 31) 

Dynamics 

Referring  to  Figure  17  , 

(3b,32)  W^'i  =  rA»_+FTHv''WAG 

(3b, 33)  Wa  X  =  Ftm  "Fad  +  ^  "  2  Fu  -  Fow 

(3b, 34)  >^'iqC!  =  Fk^Shkj-2  F^^ShM  2  SvMU. 

+  -  2  Fj^CS&u,  +SvMu.)”  Fo*g  (Z 
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Airplane  System  (3  Degree)  Equation  Flow  Diagram 


(3b, 35)  Xvvtj  ~  5hw  ^VKJ  Q 
»  *  * 

(3b. 36)  Xwki  =  X  +  SvkjQ 

Figure  19  shows  the  system  flow  diagram, 

B.  PARAMETER  EVALUATION 

Shock  Strut  Characteristics 

Figures  20  and  21  show  the  nose  gear  load  and 
damping  characteristics. 


Stroke  Zsk,('«) 

Figure  20  Nose  Gear  Damping  Curve 


Nose  Tire  Characteristics 

See  also  page  57  of  the  flywheel  system.  The  22  x  6,6-10 
16-ply  rating  nose  tire  has  a  rating  of  9150  lbs,  at  190 
psi.  The  deflection  is  1,50  inches.  The  operating  pres¬ 
sure  is  190  psi.  Since  these  are  two  nose  tires, 

(3b.37)C^jr-  (  I50\(2)(9IS'0_)  =  I2,200lt/tn 


i 


"iT 


15.3  K 


B.6>0  ir\ 


Stroke 

Figure  21  Nose  Gear  Air  Load  Curve 


Since  ^  «  0,1, 


(3b,38)  C^jT-  =  7;  \ 

-  (.i)(2K9i5o)  /'PO\  ,/ri^ 

(l.?o)^  Vl9o/\3^ 


=•  5*01(5  llo  Stftl 


The  nose  tire  rolling  resistance  coefficient  is  - 

.020  and  the  unsprung  nose  tire  mass  (mass  of  tires,  wheels, 
axlg,  and  lower  shock  strut)  is  Wwtv  “  175/386  =  .453  LBF 
SEC  /IN.  The  nose  tire  undeflected  radius, Rer/V,  is  10.8  in. 

Main  Tire  Characteristics 


The  main  tire  undeflected  radius,  Raria,  is  23.32  inches. 
The  other  main  tire  characteristics  are  computed  as  shown 
on  page  56, 


Main  Gear  Characteristics 


The  F-111  main  gear  spring  rate  parameters  were  computed 
from  load-deflection  data  recorded  during  structural 
testing  and  correlated  with  data  from  jig  drop  tests  and 
from  flight  tests. 

Figure  22  shows  the  model  which  has  the  same  form  as  that 
described  in  equations  (3b. 16)  through  (3b. 21)  and  in  the 
wheel  and  tire  system.  The  rotational  spring  rate  of  one 
main  gear  about  its  pivot  is  26  ^  ui 

The  remaining  values  are  calculated  (at  static  position)  as 


(3b, 39) 


'  -2-1 


~  .72  3 

=  VVwi4(  *  lh/*i  ■*  /  ^^7  /bf 
Cc^  -  2.00^000  \\>/ 


Thus  from  figure  22  ,  Cu.  given  by 

(3b. 40)  Cu_  :  C^,urT)/S<:^  ^  24- X  IfV  S-l""  ^  55,lXC- lb/ 

The  first  mode  natural  frequency  of  the  model  is  21,84  cps. 
Assuming  that  ^  is  .054  (about  3%  critical),  then  evalua¬ 
ting  the  damping  at  =  ^.nrrX  )  -  »  3  ?.  s' 

there  follows: 

(3b. 41)  :  ?iCo  -  (.  r<5  4)(  2.0  0, ore  )  -  ILsi't 

^l27-S)  I*') 

(3b.42)  Do  ^  =  (:0v'?K  55  ,ccy)  Zl.l  Ihjc;. 

CO  C  1 37  5  )  In 


Figure  22  Main  Gear  Strut  and  Wheel  Model 


Aerodynamic  Data 


For  finding  the  aerodynamic  data,  the  F-lllA  is  landing 
with  flaps  at  34  ,  wings  swept  to  26°,  and  spoilers  ap¬ 
plied.  ^An  equilibrium  airplane  condition  of  =  2°  and 
SuT  ="5  is  assumed.  For  these  conditions, 


(3b.43)C^^o.,s 
(3b.44)CD'  ,256 
(3b.45)CM^= 


«^i_  —  .  1 28  cLzj^ ' 


Xp  =  .000 


—  '",025 


“  *  O  2  2  cU 

^Sur 


~~.c^SZdc'^ 

T"'' 

n>bHT 


The  aerodynamic  reference  point  is  F.S.  526.8,  WL  197.2  . 
Assuming  the  airplane  C.G.  at  F.S.  519.0,  WL  180.0,  if 
Ax  andAV  are  given  by; 

(3b. 46)  AX  ^  F5A  -FSCG  ^  52(^.8-519.0  ^  7^  tnchej 

(3b. 47)  AY  =  WlA-WU:G=  =  i7.Z  jnchc-i 

Then  if  C  =  108.5  inches  is  the  length  of  the  M.A.C.,  th 
Crr^C  at  the  airplane  C.G.  is  given  by: 

(3b ,48)  C  -CfviftC  ■“Ci_AX'+^pAY 

-  Mim.s)  -(o.\3){:7.8)i-02s'ffX^^^) 

Also, 

( 3b , 4 9 )  c  ~  <) 0 c  AX  +  A 4 

i^'Xvv 

-  (,-.tJ2b)(i08,s)  -  Lae)(7.8)  f  Cc,o)U7.2)  ;  -..^7 


(3b.50) 

«^Srtr 


-  -^Ci_  Ax  t  ^Cp  Ay 
<^SrtT  c^SrtT  jSnT 

(-,O3S2)(i0g.s)  -(,022)(7.8)  -t,ti03c)(|  7,2)  r  -  3  7 


Thus  from  equations  (3b. 29),  (3b. 30),  and  (3b. 31), 


(3b. 51) 

-  -  o,,3 

Bai_ 

-  (^Ct_/fio<vV  )  “ 

<  |2S  cLo^  * 

.022  ct£^c^~' 

^AD 

-  Cp  ^  .258 

(3b.52)| 

Bad 

“  L  o  7^  ^  ” 

0.0 

^AD 

-  C  /.) - 

-.O0  2(o 

0 

'  Cam 

■=  C(v\C  -  3.424-  uTv 

(3h.53) 

j^Am 

~  (  ^Civic  / 

-  - .  3  71  in./ 

-C^C^c/^^Sht 

)  "  ~  3.  7S5^  »/f^/ ilcc. 

(3b. 54) 

^  .Gt3  -{.C2z)i~S.o)  =  -.0\(c 

(3b«55) 

^AD  =  “  E  ad  S/ir 

=  .25*?  -  CC\tO(2)  -  (-.OOBCOC'S'-^  ■“  .2.40 

(3b. 56) 

<*rtAr\  =  f^Aivi  “  ■  Eaa^Sht 

=  .3.424  -l-37l)(2.)  -{-,^7o^^X-^)  --  Z.2BC> 

initial  Conditions 

Assume  that  at  time  -  0*0  seconds  the  airplane  velocity  is 
2400  in/sec  =  Xo.  The  airplane  is  shown  in  Figure  23 
with  brakes  off. 


Figure, 23  Airplane  Initial  Equilibrium  Forces 
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Assume  that  *  3  and  “  -3^,  then  from  equations 
(3b. 29)  (3b, 30)  and  (3b, 31),  there  follows: 

(3b. 57)  Cau""  *-  (.afe)U)  1-t,C22)(-s)  :  O.llZ 

(3b.58)<^^M-  U,28(.)  <■  i-.37i')('5) +(-3.7Ss)C-S)  -  15,973 


Since  3 hk;  =  258.9,5km  =  32,6  inches,  Ith  =  20,000  in/lb. , 
and  if  the  estimated  value  for  H^t  is  97.2  inches,  then 


(3b. 59) 


(V^ 


f- 


Tru  TkM  )  (SkK)  HsT>-  Ki)iWAG  ~rAL)\  1 
(^HurSrtM)  +  Hir(p«->Jw)  / 


Now  from  equations  (3b, 24),  (3b. 25),  and  (3b. 26), 

(3b.60)(Jfl  -  t3.‘K:.o)^(52S)(,00Z36)/2fc8  =  2SOOO  |b 
(3b.61)F/yL=  (,22l)(2S006  )  -  5500  lt> 

(3b.62)T/\M^  (  19.973  )(25000  )  ^  455,300 

Thus , 

(3b. 63)  -  i  /(  .?!?  j  +  (2S7)(S7O00  -  S50o)\ 

2.  \  -t  t97.2)(0)  / 

So 

(3b. 64)  foM  22  58  8  lb 

and 


(3b,  65)  Fkj#j  -  ”  ^Av.  ~  ^ 

T  57000  -  5500  -  2  (i2  98e)  -  5524  lb 


Assume  that  when  time  -  0  that  Xwi^  =  O.'O  inches,  then 
Z&D<'^v^«>=  0,0,  Then  Xww  =  295.1  inches  so  that  Z&d<^u/w> 
=  (9,676-9,703)12  =  -.32  inches.  Refer  to  the  runway  sys¬ 
tem  for  values  of  Z(.o  .  From  equation  (3b, 12): 

(3b. 66)  S„  =(22999  )/(5S3o)  =  2.4)  lo 


Thus  from  equation  (3b. 13) 


(3b. 67)  IwMo  ^  23.52  -  ZA\  =  20.9  1  U 

From  Figure  14  in  the  flywheel  system,  if  F^ms  =  22,950 
lbs.  then  =  24.00  inches.  Now,  from  equation  (3b. 4) 

(3b. 68)  ^(5524-)/C12^ZCo)  -  .4G  In 

From  equation  (3b. 5),  there  follows; 

(3b. 69)  =  C~.32)-^-(  lO.feo)  ~  lO.OzLn 

Also,  from  Figure  21,  if  =  5,600  lbs.  then:  -  S  n 

Rearranging  equations  (3b. 1)  and  (3b. 9) 

(3b. 70)  2^  1-  Shkj  Qo  "  2w»Jo  * 

(3b. 71)  Z.0  -  S/-JM  Co  "  5  VM  ^  Z-LVMO  ”  2.SM 
Solving  these  two  equations, 

(3b, 72)  Qc  -  .0229 
(3b. 73)  2c  =  S2.3fc  tn 
Finally, 

(3b. 74)  Xo  “  Xvvfvio  “  3C),20  in 

(3b. 75)  ^  Qc  -  .0  329 

The  values  of  the  following  parameters  as  listed  in  Table 
6  are  established  by  the  airplane's  dimensional  and  mass 
characteristics:  /Idep  ^  S£:i-  ^  Sury^^  Sva/j 

Srt*  j  \Na  dhd  vOra . 

For  the  example  problem  the  density  of  air  at  standard 
conditions,  sea  level  and  59.6°F,  is  assumed.  Thus, 

Rh>^=  .00^30  Sj(/^s  /  J=t^ 

The  shock  strut  linear  damping  coefficients ,  D'v(4  for  the 
nose  gear  and  for  the  main  gear,  are  set  equal  to 
zero  for  th'^  example  problem. 
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3c. 


AIRPLANE  SYSTEM  (6  DEGREE) 


The  six-degree  airplane  system  is  built  around  a  rigid 
body  airplane  which  is  allowed  to  move  vertically  and 
horizontally  (both  parallel  and  perpendicular  to  the  run¬ 
way  centerline).  Also,  the  airplane's  yaw,  pitch,  and  roll 
effects  are  considered.  This  model  considers  all  the  ef¬ 
fects  found  in  the  three-degree  airplane  system.  The  pur¬ 
pose  of  the  six-degree  airplane  is  primarily  two-fold:  the 
first  is  to  evaluate  the  effects  of  the  anti-skid  system 
on  the  airplane's  directional  stability;  the  second  is  to 
evaluate  any  anti-skid  system  degradation  caused  by  air¬ 
plane  yaw  and  side  drift  movement. 

For  the  nose  gear,  the  model  considers  the  tire  and  strut 
characteristics  in  the  vertical  direction.  Also,  the  nose 
tire's  yawed  rooling  characteristics  are  included.  The 
steering  loop  is  closed  by  providing  a  "pilot"  function 
which  provides  an  input  to  the  nose  tire.  The  "pilot" 
function  depends  on  the  airplane's  yaw  angle.  The  two 
main  gears  are  treated  as  two  distinct  systems  except  for 
any  structural  coupling  which  may  exist  between  the  two. 
Provisions  are  made  for  side  wind  perturbation  and  for 
aerodynamic  effects  caused  by  airplane  yaw  and  roll, 

A,  Mathematical  Description 

Figure  ?.4  shows  the  six  coordinates  which  describe  the  air¬ 
plane  position  relative  to  reference  points  on  the  earth's 
surface. 


Figure  24  Airplane  Coordinates 
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Vv^y  is  a  crosswind.  The  runway  is  oriented  so  that  its 
centerline  coincides  with  the  x  axis  for  0  inch  runway 
heights  (Zgd  =  0).  This  analysis  assumes  that  the  pitch 
(Q)  and  roil  (P)  angles  are  small.  Let  Z*o<x,y>denote  the 
runway  profile  and  let  Z4op<x,y>denote  the  runway  slope 
(ZcDP<x,y>  =  ^  x,y>/^x).  Figure  25  shows  the  airplane 

as  measured  in  the  fuselage  station-water  line  reference 
system. 

Nose  Gear 

Let  2,^  and  Z^^  denote  the  nose  gear  stroke  and  stroke 
velocity.  Then  we  have  that: 

(3Cgl)  ^  ^VM  “  Z  “  Sm*jQ 

(3c. 2)  *  ZwM  -  i  -  5h..Q 

The  nose  gear  shock  strut  force  Fvm  is  then  given  by: 

(3c, 3)  ^VM  ^SkJ  ^Sn)  1  ^  SW  I 


Figures  26  ,  27  ,  and  28  show  the  forces  acting  on  the 
airplane  as  seen  in  the  different  planes.  Let  Fln  denote 
the  lateral  force  on  the  nose  wheel  at  the  axle.  Then: 

(3c, 4)  Wwrg  Vy  =  ^SN5  ~ 

Where  Fsms  is  the  lateral  component  of  the  sliding  or  cor¬ 
nering  force  of  the  nose  tire.  The  load  F^w  is  caused  by 
the  nose  wheel  trying  to  move  laterally  relative  to  the 
airplane.  If  this  lateral  displacement  is  denoted  by  YpLW  » 
then: 

(3c, 5)  Pusi  ~  ClKJ  VolKJ  XplW 

(3c.6)youM=  yN.-y  ^-(^.+SHK.Q)p  - 
(3c. 7)  (I:  +  S„k,Q)P  +  (i+SHM(3)P-SHMR 

Now  F,^^  is  given  by; 

(3c, 8)  “  Sw  (  CmT  ) 

where 

(3c,  9)  5m  =  rT)a.x  C.O  j  ~  ^  w»j  1 

(3c  ,  10  )Sm  =  ^  '^WKJ  ""  2  WN( 

Summing  vertical  forces  on  the  nose  gear  unsprung  weight: 
(3Cell}  —  FVj  )j  — 

Assume  that  the  pilot  positions  the  nose  wheel  with  a  rate 
proportional  to  the  airplane  yaw  angle.  Thus; 

(3c,  12)  ©w  =  I  mm{o,  -GpxlR  ©n  ^ 

j  “  R  *  i  1  ^  I  ©WN)AX  1 

1  -GpilR  &Ki  -  “©N.'-ORX 

©w  gives  the  yaw  angle  of  the  nose  wheel  with  respect  to 
the  airplane  ,  The  yaw  angle  of  the  tire  with  respect  to 
its  direction  of  motion  is  given  byGy^y/, 

(3c.l3)  -t  R  -  (  Ykj/ Xv-n) 


The  steering  characteristic  is  developed  from  Reference  1 
(p.  30).  Let  Ukjtf  the  coefficient  of  friction  between 
the  nose  tire  and  the  ground.  Then  the  maximum  force  nor¬ 
mal  to  the  tire  in  the  plane  of  the  ground  is  Fmtf  where: 

(3c, 14) 

Using  equation  (79)  and  80)  from  Reference  1  1 

(3c.l5)  IIrf  "  I  ^wc  ®VA  W if  >o 
(  o  i  f  F^J7-F  -  ^ 

(3c.l6)  Fmcfs"  Fs,rF  if  Uet  ^  )-S 

C  Urt  ”  UEr/27  )  i-f  I  Ugr  1  <  i.  s 

-  a  Llpr  -  -‘-S 

Thus,  F^gcFs  corresponds  to  in  Reference  1  and  Pwc 

is  the  cornering  power  given  by: 

(3c,  17)  Sm-Spj 

|Cp3“Cp4  5jj  if  Sw^Spi 

The  actual  normal  cornering  force  F^cr  is  not  F^cfs  j  but 
lags  Fncfs  because  of  the  tire  relaxation  length.  The  ex¬ 
pression  for  Fk,cf  is  given  by: 

(3c, 18)  *  C  Fncfs  “  ^cp)(X  ww  /£v«u) 

Having  obtained  F^,cF  >  then  from  Figure  29  , 

(3c, 19)  Fsvjs  "  ^CF  "  UEew 

(3c,  20)  ^ R ^  F  Uegfj  c/ru  R  / 


Main  Gear 


Let  Zsr^  and  Zsm  denote  the  stroke  and  stroke  velocity. 
The  additional  subscripts  L  and  R  refer  to  the  left  and 
right  side  of  the  airplane  (looking  forward), 

(3c.21)  ^SMI?  ~  ZwMR.  "  ^  ^VM  P 

(3c,22)  “  r.WN»tt  "  "2.  +■  ShmQ  ■*’  P 

(3c,23)  ZjMu  ”  “iv/ML  "■  2  5 ShmQ  “  P 

.  •  •  • 

(3c624)  ^5mu  '  2ttfML  "■  2  ShmQ  “ 


The  main  gear  shock  strut  forces  are  then  given  by: 

(3c. 25) 

•  •  • 

(3c,  26)  =  F 2 

4ml  y  ■•"Dvm  2sml'*'  2sml^  2-Smi_  1  I 


Let  Sm  denote  the  main  gear  ti'e  deflection  and  let  Fkim 
be  the  associated  load.  Thus,  in  the  vertical  direction, 
the  relation  between  the  load  and  tire  deflection  is  given 
as  follows: 

(3c. 27)  (Cmt 

(3c,28)FijML  '  ^ML.  (^MT  ^  ) 

(3c,  29)Sprf)/j_  =  Max  {  o,  2<;d^  ^wi^r.  ,'^mr^  ^  ~  2  wroR.  } 

(3c.30)Smr.  <x  WMHiVmr)  ^WMR  “  ZsvME. 

(3c,31)Sj^^  '  ^  -fw  ML  1 

(3c.32)Sml  =  ZsDP<X  WMLj^ML^  ^WML 

Summing  forces  in  the  vertical  direction  on  the  main  gear 
wheels , 

•  » 

(3l,  33)  WwM  ~  ^IMC.  ”  ^  ^R-VR 

(3c  o  34)  WwM  2  WML  "  ~  ^rMU  Fsrv/l 


Figure  30  shows  a  side  view  of  the  left  hand  main  gear. 
With  the  assumption  that  W^  ,  ©^r.  and  ©s^.  are  des¬ 

cribed  by: 
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I  - 1 . • . 


where 


(3c*35)  Wu  5feu  Qgb  "  Fur  ■*"  (  Ftr  "  )  -Ts;. 

(3Cb  36)Wn  Sgu^SL  ~  Sttu  FijL  Cs&u.  26LlX  f>E~  ''Tl“  ~Tsl. 

(3c,37)  ^GL&  ~  ^tL.“ 


(3c, 38)  “  ^fiL.  "■  2&MI- 

The  forces  Fr^  and  are  used  to  impart  the  correct  mo¬ 
ment  into  the  gear. 


The  overall  gear  system  model  is  shown  in  Figure  31  ,  In 
order  to  transmit  torque  properly,  the  forces  Fre^  and 
Ftl^  are  applied  equal  and  opposite  on  different  sides  of 
the  gear.  Thus, 

(3c ,  39)  Ftr©  "  F&r  ^  S&w  "  ^&s  )  /  2  S&s 

(3c,  40)  Ftl©  ~  PsL  (  ^6W  ”  )/  2  s 

If  it  is  assumed  that  100  percent  of  this  torque  is 
taken  directly  into  the  airplane,  then  100  =  100-100 

percent  is  transmitted  through  the  gear.  Thus, 

(3c, 41)  H&  Pre^ 

(3c. 42)  Ftl.  =  Hs  FrLs 
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Spring  System 


All  spring  loads 
positive  in  tension 


Let  Qr  and  be  the  difference  between  the  gear  rotation 
and  the  airplane  rotation.  That  is, 

(3c. 43)  Qr  =  Q  - 

(3c. 44)  Qu  =  Q  -  e^L 

(3c. 45)  Qr  =  Q  -  06R 

(3c. 46)  Ql.  =  Q  ”  0&L 

Then  we  can  find  constants  Cui  ,  ^uz  such 

that: 

(3c. 47)  Cut  Qc'^uzQu) +  (0U|Qe  ” 

(3c. 48)  ~  S,ju.  (Cai  Q  L  ” CuzQ  a)  CDuiQl  “ 

It  then  follows,  assuming  negligible  strut  moment  of  in¬ 
ertia  that: 

(3.49)  ^t>uz  ((F^r  ■*'  "  ^Tl)  ^SLR.  "^E  )/  S&u. 

(3c. 50)  ^DU-L  =  ((F.,.  ^Vl”  ^R  )  ^6LL  ~^Su  )/  ^<sa 

As  outputs  to  the  tire  and  wheel  systems  we  need  to  com¬ 
pute  Xax  and  Y^x  •  ^Ax  shown  in  Figure  30 .  Y^x  is  as¬ 
sumed  to  be  the  undeflected  tire  footprint  position  in  the 
y  direction. 

(3c.  51)  ^AXL  "  X  SvmuQ  '*■(^611.'^  ^SLL)^tL  ”  ShM 

(3c. 52)  )^Axc.  " X  "  ^sw ^  ■*"  ^vmuQ  ^SLa^Sea”’ Shm 

(3c.53)Xax  =  y  +  Sgw  Svmu(3‘*‘  (Sfia'^'ZgLL 
(3c. 54)  XAxa'^X-  Sau/R  +Sv/mu(5  +  (Scu."*"  ^&LB  )06R 
(3C.55)  Yax^^  y-S^;^-(gHM'SvMuQ-(S6L.-^Z6LLi©*u)R 

(SvMU  ■*"  ^Gu.  Z<i(.L  3  P 

(3c.  56)  YAjce."  Y'^S&n/  ~  (Shm  "SvmuQ  "  (Ssu"^  ^6ur3©cr.3R 
'^S/MU.'*'  ^Gtc  Zgle)  P 

(3c  .  57)  YaXL  ”  Y  ~(Shm“  SvMU-Gi  "(Ssu+  Z6Ll)©&l)  P 
■*'CSvmU-Q  ■*“  (  Z.sLl)06l)  P 

”  t  Scu.  ■*"  Z.g,i.L  )  P 

91 


(3c.58)yAXR 


y  ('Shm“SvmuQ  “  (Scu. +  2:61,12.) ©fie)  ^ 

(LU.+  2-6Le  )^iSB.  )  ^ 

"  (  ^VMm.  ^^Sll  +  2.&lk3  P 


Engine  Thrust 

Referring  to  Figures  26  and  27  ,  if  Fth  is  the  engine 
thrust,  then 

(3c. 59)  ^Hv  “  Pth  ( ®^TH  +  Q) 

(3c,60)  TVh  -  Pth 

(3c, 61)  Pths  ~  ^ 

Aerodynamics 

The  following  eight  equations  apply  as  in  the  three- 
degree  model. 

(3c, 62)  Qfi,  ==  X  ^ha/^ 28$ 

(3c. 63)  F>^  =  C^lQa 

(3c*64)  Fap  =CApOft 

(3c, 65)  -  CA^lOA 

(3c,  66)  0<w  =  ex'©  +•  Ci8o/tt)(  Q  -  fe/y  ) 

(3c, 67)  Cai,  “  ^Al.  Bal*^W  ■*“  EAL^Hr 

(3c, 68)  =  ^ad  Bado^w  +-  Eac?  S^t 

(3c. 69)  Cam  "  Cam  "P  +  ^am^ht 

Let  Vwv  denote  the  wind  gust  velocity  as  shown  in  Figure  24 
If  4^  and  f3  are  defined  by: 

(3c. 70)  =  (  Vv^v  +  y)/x 


(3c. 71)  /3  =  Oso/TT  )(4^-  R  ) 
Then  /3  is  the  angle  of  sideslip. 


Let  Qat  denote  the  dynamic  air  pressure  (including  side 
wind)  multiplied  by  the  reference  area.  Then; 

(3c. 72)  ((Vwv,^-y)^+-X^)AeeF!?HA/z88 

Then  the  aerodynamic  yaw  moment  is  given  by: 

(3c, 73)  Tav  “  Can /3  Qat 

and  the  aerodynamic  side  force  is  given  by: 

(3c. 74)  F^sy  ^Cav/SQat 

Finally,  an  aerodynamic  force  Fasu  to  a  combination  of 
lift  and  pitch  is; 

(3c. 75)  Fasl  =  F>l  P 

Refer  to  Figure  27  as  to  the  direction  of  these  forces. 
Dynamics 

Referring  to  Figures  26  and  27  ,  summing  forces  in  the  x, 
y,  and  z  direction, 

(3c, 76)  Wa^  ^  +fvME  +  + 

(3c,  77)  Wp,X  -  Fth  "  Fap  1-  +  Foul  “  "" 

(3c,  78)  W^  y  ^  Fist  “  Fasv 

Summing  moments  about  the  C,G.  we  have: 

(3c , 79)  V/cQ Q  "  R-nj  ^jiM  “  ^  ^ v'mu 

R>UL^  LL  "Rh  "Rm  “  Fib  ^  ^SU.  ^  VM  u.  ) 

~  (S(iLl+'  5y^/|  Lt  )  “  RfJ  (  Z  "  H.CJD  ^  ^  VVKJ  ^  3 

(3c, 80)  Wji^  S  =  FLn^hw  ^  C  Fut  “  Ril  ‘“oul"  F*Due3 

+  2  Ma  (  F'tr.®  -  Ftlo  ) 

(3c  ,81)  Y/jp  P  =  (Fvmu  “  Fyjvjg)  +  ("asy  Hap 

“  C  ^  <x,y>)C  FsTR.  F &TL  Fi_n  ) 

(3c,82)  -  X  ^  SvmQ 

(3c, 83)  =  X  Q 
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Airplane  System  (6  Degree)  Equation  Flow  Diagram 


B*  Parameter  Evaluation 


Nose  Gear  Characteristics 


Based  on  a  nose  gear  lateral  natural  frequency  of  12  cps, 
we  have  cOn  =  2Tr(12)  =  75*5  RAD/SEC.  Since  Www  =  .435,  thei 

(3c.84)Ci.n  =  Wwkj<^«  =  .435- (75.5-)^  =  2+80  Ib/in 

Using  =  ,054  as  in  the  calculation  of  Dq  in  the  three- 
degree  model. 

(3c. 85)  =  C.05-4)(2+80)/7S.5  =  1.78 

The  steering  or  cornering  characteristic  parameters  are 
obtained  from  Reference  1,  Based  on  Figure44(a)  in  Refer¬ 
ence  1,  the  value  for  U^rp  is: 

(3c. 86)  U-urp  =  f>,r,a(r>iax)/Fa  =  ZS-ooo/asZoo  -  ,5^3 

Using  equation  82  from  Reference  1,  if  K  =  CP+-.4^/^)w*  = 
(1.44)(190)(6.6)2  =  11,920  lb.  then: 

(3c.87)Ci’i  =  1.2.  C.K/J  •  ’/Mo  llo/HxA  ih. 

(3c.88)c«  asctK/d*-. aiss  ih/R&d ik 

(3c,91)Spt  =  .0875  d  =  (.C>875)f22)  =  1.925  in 

From  Figure  43  in  Reference  1  we  see  that  the  cornering 
force  lags  the  yaw  angle.  Equation  63  in  Reference  1 
shows  that  the  equation  which  describes  the  curves  in 
Figure  43  is  given  by; 

(3c,  92)  Pyj  r  =  (  I  “  f  )  /vtax  ^  ^ 

where  Ly  is  the  tire  yawed  rolling  relaxation  length. 
Differentiating  this  equation,  there  follows: 

(3c,  93)  d^jr-  =  G  Fy  r  rvtay  ^  ©^V(Sw) 


-Oc/^y 

Eliminating  e  results  in: 

(3c,94)  Fy  r  Ly  d. Fy.K'  =■  V  noax 

da 

Equation  (3c. 18)  is  obtained  by  using: 

(3c,  95)  d  Py,r  =  d  Fy^ir  /  doc  ^  ^^ci=~ 
dx  dt  /  dt 

where  it  is  assumed  for  large  airplane  velocities  that 
Xvvw  ”  dx/dt.  We  see  that  the  parameter  Sy^l  is  the  relax¬ 
ation  length.  From  Figure  39  of  Reference  1,  for  most 
conditions,  obtained  from: 

(3c.  96)  SyEt.=  .6uj(2.S  -.9  P/Pr) 

=  (.fc)(t.^)(2.g -.S)  -  7.3Z  In 


Main  Gear  Characteristics 


For  many  airplanes  which  have  a  conventional  strut  arrange¬ 
ment  (similar  to  a  B-58)  most  of  the  moment  about  the  shock 
strut  is  taken  out  through  the  shock  strut.  In  this 
case  equations  (3c. 41)  and  (3c. 42)  would  use  Hg^o.o. 

In  the  case  of  the  F-111  gear  the  opposite  result  occurs  so 
that  Hs  •  1.0  and  -  0,0.  The  following  values  apply  to 
the  F-111  gear: 


(3c. 97) 


1 


Wu.  =  .723  \\>  5€c^/ in 

WwM  ■  Wwv  -  lb  sec^/in 

Ss-L^  -  ei.oo  in 

S<sw  =  4>0.oo  Iv) 

^65  =  20.00  to 

Hg.  =  uo 
H  A  =  0.0 


If  loads  ■>  t)  =  ^Ct)  =  Fv  are  applied  as  shown  in 
figure  31  ,  then  because  of  symmetry,  the  result  will  be 

that  Qe  =  0l  .  But  then  equation  (3c. 47)  says  that 
Cu)~Cu2  -  / Sfiu  Qr  but  Co-  ^ / Sfru,  Q a 

as  shown  in  the  3  degree  model.  Thus 

(3c. 98)  Cui-C^i  =  Cm.  =  59000  It/ I'n 
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With  the  main  gear  at  static,  if  a  drag  load  of  18,000  lb, 
is  applied  to  the  left  gear  at  the  ground  and  -18,000  lb, 
is  applied  to  the  right  gear  at  the  ground  the  observed 
deflections  with  Q  ■  0  are  Ql  "  .0236  rad  and  Qr  “  -,0236 
v^aid,  (Assuming  a  lateral  beam  torsional  spring  rate  of 

43.0  ^  ]0‘^  In  lb/ rad  ). 

In  the  equations  which  describe  the  gear  loading  ~^stz  and 
can  be  chosen  as  0  if  and  are  the  dimensions 

to  the  ground  instead  of  the  axle.  Thus  Z&u  ■  •2.g.le  = 

2.2  +  12,9  =  21.4  in.  Equations  (3c. 35),  (3c, 36),  (3c. 39), 
(3c, 40),  (3c. 41)  and  (3c. 42)  can  then  be  combined  to  give 


(3c. 99)  Fu,2  -Fol 


=  -2 12^000  lb 


Subtracting  equation  (3c. 48)  from  (3.47)  results  in 
(3c. 100)  =  (C  ~  (Cui  +  Cuz) 

So  that 

(3c,101)Cu,  =  -2 >2000  =  214,000  Ib/in 

(2)(2I.O)(-.023^ 

Adding  and  subtracting  equations  (3c. 98)  and  (3c, 101) 
results  in 

(3c. 102)  Cw.1  -  5300^  1-^14  000  =  /3fcjS-00  'b/Ln 

2 

(3c.103)Cu2  =  214,000  -  55000  =  7?^  SOD  Ib/in 

2 

At  a  fore  and  aft  natu;*al  frequency  of  137,5 

the  damping  coefficients  D^i  and  are  given  as 

(3c. 104)  Du.  ^  /^Cui/uj  =  (■0^4)(l.3fa>tlQ^)  =  53.4  ib  sec/ln 

(137.5  ) 

=  C.054)(.77SxiO^^^  30.5  lb  se.c/ln 

(i37,5  ) 


(3C.105)  duz 


AerodynaiaxC  Characteristics 


The  coefficients  for  equations  (3c«62)  thru  (3c* 69)  have 
been  derived  in  the  3  degree  model.  For  the  F-lllA  in  the 
landing  configuration  and  wings  swept  to  26  degrees  as  des¬ 
cribed  in  the  3  degree  system,  Cwg  «  .0014  and  Cy^  -  -.021. 
Then  the  coefficient  is  calculated  from 

(3c.  106)  =  ~^Y/3  - 

Let  A  X  »  FSA  -  FSCG  as  in  the  3  degree  system  where 
FSA  =  526.8  and  FSCG  “  519.0.  Let  b  be  the  wing  span. 

If  b  =  756  in.,  then 

(3c.  107)  Cam  =  bCwa  -  AX  Cy/j 

(3c. 108)  Can  =  (75C.)(.OCl4)  -  (7. eo)(--,02 1  )  =  I.2Z2 

Airplane  Characteristics 

The  parameters  listed  in  Table  7  describing  the  airplane's 
dimensional  and  mass  characteristics  are  those  previously 
derived  in  the  3  degree  model  or  simply  a  listing  of  the 
appropriate  values  applicable  to  the  F-111  for  which  no 
derivation  or  computation  is  required. 
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4a„  WHEEL  AND  TIRE  SYSTEM  (FLYWHEEL) 

Figure  33  shows  the  components  of  the  wheel  and  tire  system. 


Figure  33  Components  of  the  Wheel  and  Tire  System 


In  the  vertical,  or  Z  direction,  the  axle,  brake,  wheel, 
tire,  and  lower  shock  strut  are  combined  and  operate  as  a 
single  mass  points  A  description  of  this  mode  is  found  in 
the  airplane  system.  The  airplane  system  furnishes  various 
inputs  to  the  tire  and  wheel:  Vf  the  airplane  (flywheel 
surface)  velocity;  fv/M  ,  the  vertical  load  between  the  tire 
and  pavement; 5h  >  the  tire  deflection.  The  brake  torque 
Tbt  is  an  input  from  the  brake  system. 


The  horizontal  displacement  of  two  mass  points  is  con¬ 
sidered,  One  mass  point  is  made  up  of  the  axle,  brake, 
wheel,  and  the  inner  part  of  the  tire  and  its  location  is 
designated  asXw.  The  other  mass  point  is  the  tire  tread 
and  its  location  is  designated  asXrr* 


In  rotation,  there  are  three  mass  points:  the  axle  and 
stationary  brake  elements  make  up  the  first;  the  brake 
rotors,  wheel,  and  inner  tire  make  up  the  second;  and  the 
tire  tread  makes  up  the  third.  The  angular  positions  of 
these  three  mass  points  are  denoted  respectively  as  Os, 
Ow^and^T*  be  the  horizontal  force  acting  on  the 
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axle  and  let  Pyr  be  the  net  horizontal  force  between  the 
wheel  and  tire  treads  Figure  34  shows  the  location  of 
these  forces e  Fsr  is  the  horizontal  force  between  the  tire 
and  the  flywheel  surface. 


Figure  34  Tire  Horizontal  Model 
A,  Mathematical  Description 

Equations  describing  the  tire  and  wheel  behavior  are  de¬ 
veloped  by  referring  to  Figure  34  .  Forces  Fs  and 

are  defined  by  equations  (4a, 1),  (4a, 2),  and  (4a. 3)  as 
follows: 

(4a, 1)  Fg  -  "Cgh  Xw  “DtfdXw 

(4a,  2)  "  CTr(XTr'Xv^/)  E-,t  (  X-fx  ”  ^ 

(4  a,  3)  Dtt  (  Xy  ”  )  =  ETr^XrT“Xy) 

Equations  (4a, 2)  and  (4a, 3)  describe  a  type  2  spring- 
damper  system  as  defined  by  Figure  38  and  discussed  in 
the  parameter  evaliiation. 
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Let  Wsw  denote  the  mass  of  the  axle,  wheel,  brake  and  inner 
part  of  the  tire.  Let  Wte  denote  the  appropriate  tire 
tread  mass.  Summing  forces  in  the  horizontal  direction 
giv2s: 

(4a. 4)  Wgvs/  r 

(4a. 5)  Wrg  Xtt  “  ”  F’tt  “  ’’’ 

Where  is  a  force  produced  by  tire  unbalance,  the  cor¬ 

responding  vertical  part  of  this  unbalance  force  is  de¬ 
noted  byFsBv  -  These  two  forces  are  given  in  equations 
(4a. 6)  and  (4a. 7). 

(4a. 6)  Vn/;- 

(4a.?) 

Wrander  are  the  rotational  speed  and  position  ot  the  tire 
tread. 

The  rotational  schematic  of  the  wheel  and  tire  system  is 
shown  in  Figure  35  . 


Let  Trt  and  Tj  be  defined  by  equations  (4a. S),  (4a. 9),  and 
(4a. 1C)  as  follows: 

(4a. 8)  Tp-j-  = 

(4a. 9)  DHr(0v~^T)  “ 


no 


(4a,  10)  7^  -  C ,50^65 

Let  Hr  be  the  height  of  the  axle  above  the  ground.  Let 
be  the  torque  on  the  tire  that  produces  rolling  resist 
ance.  These  two  quantities  are  given  by: 

(4a,  11)  Hr  = 

(4a, 12)  T'ee  =  j  Dyjj  Wj*)  if  NVr  ^  ® 

o  i£  Wr  =  o 

Sm  (~Dsr  +DvRVr)  if  Vr  <  o 


If  Tqt  is  the  brake  torque,  then  torques  can  be  summed  to 
obtain  the  following  three  equations: 

(4a, 13)  =  T3r"Ts 

(4a, 14) 

(4a.  15)  W^rBr  =  ^rf^r^T^r  Tze 

The  rolling  radius  of  the  tire  is  obtained  using  the 
methods  of  Reference  1.  Denoting  the  rolling  radius  as 
Ry ,  it  is  defined  as: 

(4a,  16)  Rr  "  Xrr  “  ^ 

Let  Vps  denote  the  velocity  of  the  tire  footprint  ,  Rela¬ 
tive  to  the  flywheel  surface  with  WV  =  0,  let  be  the 
relative  velocity  including Wy, 

(4a,17)  Ves  =  Vp’  +  ;<rr 

(4a,  18)  Vti  -  Ve5  “  RrV/r 

Here  Vp  is  the  velocity  of  the  flywheel  surface.  Adopting 
the  convention,  V/y  ;  W,i  =  ;  and  =  ©w>  the 

relative  angular  velocity  between  the  stators  and  rotors 
is  denoted  by  Wb  and  is  established  by: 

(4a,  *.9)  Wb  =  Ww  “Wj 


When  a  tire  is  moving  over  a  runway  with  any  appreciable 
amount  of  standing  water  or  slush,  a  hydrodynamic  "wedge" 


Figure  36  Wheel  and  Tire  System  (Flywheel)  Equation  Flow  Diagram 


of  V7ater  starts  separating  the  tread  and  runway  surface. 

It  is  assumed  that  the  length  of  this  "wedge"  is  propor¬ 
tional  to  Vel  and  at  hydroplaning  speed,  V^y  ,  the  tread 
is  completely  separated  from  the  runway.  In  equations 
(4a. 20)  and  (4a, 20)  the  coefficients  Chv  and  Dwy  are  used 
to  define  hydroplaning  effects  and  water  drag  on  the  wheel. 
For  dry  runway  conditions ,  Chv  and  Dhv  are  zero.  The  hori¬ 
zontal  force  between  the  tire  tread  footprint  and  the  run¬ 
way  surface  is  established  by  equations  (4a, 20),  (4a, 21), 
and  (4a.  22)  a.'^  follov?s: 

(4a.  20)  ^  “  ^HY  ^  ^es/ ) 

(4a.  21/  Fg-p  = 

(4a, 22)  U.T  -  U.TI  +  (  Ut2 "  ^tWs)  ^ 

o 

'  -  Uri  -  (  Ut2  ~  Et'^Rs)  ^ 


if  Vg  >0 
if  Ve  = 
if  Vg  ^  0 


Figure  36  is  an  equation  flow  diagram  showing  the  re¬ 
lation  between  equations  (4a, 1)  through  (4a, 22), 

B,  Parameter  Evaluation 

Gear  Characteristics 

The  mass  W^w  is  made  up  of  the  mass  of  half  the  shock  strut 
half  the  lateral  beam,  the  axle,  the  wheel,  tiie  brakes, 
and  all  but  one-third  of  the  tire  tread.  The  sum  of  the 
masses  of  these  components  totals  616  LBM.  Thus,V6w  = 
6IG/3F4.  -  1,4,0  lb  5ec^/  m  ■  The  fore  and  aft  natural  fre¬ 
quency  of  the  gear  (as  calculated  from  deflection  data)  is 
Zi  e^cps  =  137,5  rjid/scc  .  Using  the  gear  mass,  with  all  of 
the  tire  included  (644  LBM) ,  the  spring  rate  C&ncan  be 
calculated  as: 

(4a, 23)  Cgs=  moji  =  =  31,500  Ib/in 

\38<oJ 


A  typical  approach  to  estimate  the  damping  coefficient 
is  to  use  3%  critical.  Thus, 


(4a. 24)  (.03)zy^fi  ~  13.8  Ibs^c 


Tire  Tread  Characteristics 


The  principle  underlying  the  calculation  of  the  tire  fric¬ 
tion  coefficient  is  that  compared  to  the  rest  of  the  tire, 
the  tire  "footprint"  is  totally  inelastic,  (The  tire  "foot 
print"  is  that  portion  of  the  tire  tread  which  is  in  con¬ 
tact  with  the  ground).  Thus,  if  the  vi  oclty  of  the  foot¬ 
print  and  the  friction  vs,  velocity  curve  for  the  rubber- 
surface  interface  are  defined,  the  tire  friction  coeffi¬ 
cient  is  established.  In  order  to  predict  the  motion  of 
the  footprint,  the  tir^'.  tread  is  assumed  to  behave  like  an 
inelastic  ring  which  is  supported  on  the  wheel  as  shown  in 
Figure  37  . 


Torsional  springs 


Inelastic  ring 


Translational  springs 


The  horizontal  position  of  the  footprint  is  assumed  to  be 
the  same  as  the  h^'^rizontal  position  of  the  ring  center  of 
gravity.  A  frictional  force  applied  at  the  ground  can 
be  resolved  into  a  tranjlational  force  F  ■  Fh  acting  at 
the  ring  C.G. ,  plus  a  moment  M  *  Fh'H  acting  about  the  ring 
C.G. 

For  an  actual  tire  with  distributed  mass  and  elasticity  ap¬ 
proximately  one-third  of  the  tire  would  move  with  the  foot¬ 
print  in  response  to  force,  Fh.  Therefore,  it  is  assumed 
for  translation  the  mass  cf  the  rlng,WrE  »  one-third  'f 
the  tire  tread  mass,  which  is  84  LBM.  Thus,  Wte  ■=  (84)/^3x 
386)  =  0.0725  LBF  SEC^/IN. 

For  rotation  the  total  tire  tread  mass  is  assumed  to  move 
in  response  to  moment,  M.  Thus,  the  mo-^nt  of  Inertia  » 
about  its  center  of  gravity  is  Wn  =  MR^  *  (84/386)  (23  )  = 
115  LBF  SEC^/IN. 

References  1  (page  22)  and  10  (Figure  8)  are  used  to  obtain 
values  for  the  torsional  and  translational  spring  rates 
as  shown  in  Figure  37  .  Under  the  application  of  the  force 
Fh  ,  the  peripheral  movement  at  point  b  is  about  20%  of  the 
peripheral  movement  at  point  a  (Figure  37  above).  The  ex¬ 
pression  for  the  f^  otprint  spring  rate  from  Reference  1  is 

(4a. 25)  K:*  =  <t-Pr)^So/d 

Where  for  the  F-111  with  a  vertical  tire  load  of  25,000  lb, 

d  =  46.65  in.  =  Tire  diameter 

P  =  150  psi  =  Tire  operating  pressure 

Pr  =  150  psi  -  Tire  rated  pressure 

So  =  2.75  in.  =  Operating  (static)  deflection 

Thus, 

(4a.  26)  Kx  =  C,fc)(‘Ii..65)(5-)(;50}  ^2,75/44.65  -^Ifo /t/,„ 

The  application  of  Fh  =  8150  lb.  causes  the  footprint  to 
move  one  inch.  At  point  b,  the  movement  is  .2  inches. 
Assuming  that  the  movement  at  point  b  is  all  due  to  rota¬ 
tion,  the  apparent  torsional  spring  rate  is: 

(4a.  27)  Crt  “  -  Czs.S2)i. zaS'7X0^^o)  fiz/tpyrld 


Since  .8  inches  of  the  1.0  inch  footprint  motion  is  due  to 
tread  C.G.  fore  and  aft  translation,  the  apparent  spring 
rate  is; 


(4a.  28)  CjT  ~  ^  =  91^0  =  10,200  \y/ii 

.8  .8 


Figure  38  Tire  Damping  Models 

It  is  well  publicized  and  generally  accepted  that  the 
elastic  and  damping  characteristics  of  tires  and  other 
structural  devices  are  not  accurately  described  by  the 
mathematically  convenient  linear  spring-viscous  damper 
representation  over  a  wide  frequency  range.  The  behavior 
of  rubber-like  materials  is  particularly  (iifferent  than 
that  described  by  the  conventional  m-rdel.  To  establish 
suitable  mathematical  descriptions  of  the  various  damping 
forces  for  tires  and  other  elements  of  uUis  study,  several 
models  were  explored.  Figure  38  depicts  the  two  types  of 
elastic  systems  which  are  used.  The  Type  1  model  is  a 
conventional  system  with  viscous  damping  and  Type  2  is  a 
visco-elastic  system,  having  elasticity  and  damping  which 
varies  with  frequency.  To  compare  the  two,  consider  the 
effects  of  driving  each  with  a  variable  force  F(t)  =  Fo 
.  In  each  case,  the  resultant  deflection  is 
s  oco  oe^(ujt  -  w) 


The  loss  coefficient,  /3,  is  defined  by  /3  ,  For  a 
conventional  system  (Type  1  withe  and  k  constant),  the 
loss  coefficient  which  is  a  measure  of  the  damping  is 
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given  by: 

(4a.29)  /3  =  Cco/k 

Reference  1  assumes  that  c  is  of  the  form  c  =  '^  k/oj  where 
and  k  are  constanto  From  this  ; 

(4a.  30)  /?>  = 

Reference  8  seems  to  indicate  (p.  55)  that  the  loss  co¬ 
efficient  for  tire  tread  rubber  is  somewhere  between  the 
above  two  values. 

For  the  Type  2  system,  shown  in  Figure  38, with  c,  n,  and  k 
constant,  the  loss  coefficient  is  given  by: 

(4.. 31)  ,3  .  (=«)/(,  ,  (3rj(,,g)) 

To  represent  a  tire,  the  values  (c/k)  =  1,56  x  10”^  sec,  and 
(n/k)  =  0.520,  were  used  to  compute  values  of  /3  for  a 
range  of  frequencies.  /3  verstxs  co  is  shown  in  Figure  39 
for  both  Type  1  and  Type  2  models,  along  with  values  of  /3 


Loss 


The  tire  elastic  and  damping  coefficients  are: 


(4a.32) 


Crr  =  I5--9  lb  s«c/in 
Err  ^  (.51)^77  -  5300  Ib/in 
Dft7  =  C  Get  -  30^400  LW  ii?  5fic/»'ad 

Ecr  =  Ls^)Ca,T=  lo.tsxio*’ in 


Frequency  •  rad/sec 


Figure  39  Model  Loss  Factors 


The  equation  (4a«16)  for  the  rolling  radius  Rt  is  a  re¬ 
statement  of  equation  (76  b)  of  Reference  1*  To  allow  for 
circumferential  decay  length  other  than  those  equal  to  the 
outside  free  tire  radius,  a  coefficient  is  provided. 
For  this  study  Urr  is  set  equal  to  1.0. 


Axle  Parameters 


The  observed  torsional  natural  frequency  of  the  axle  (with 
brake  stators)  is  125  cps.  The  calculated  value  for 
its  moment  of  inertia  is  16.8  LBF  SEC^/IN.  Thus,  the 
torsional  spring  rate,Ces  >  is  established  as: 

(4a.33)Ce5=  (  2Tr  125)^  ( ir^lb/rad 

For  the  steel  axle,  a  value  for  7  (in  the  Type  1  system  in 
Figure  38)  is  probably  something  less  than  .01  (Reference 
7).  Thus,  at  resonance,  if  cco/k  =  then  the  damping  co¬ 
efficient  is  established  as: 

(4a. 34)  Dgs  =  k7?/oj  =  Cio.4x  io‘’)(.oi)/C3tti2s)  -  »32  Inlkeec/rtiJl 
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Tire  Rolling  Resistance 


From  Figure  17a  of  Reference  2,  the  rolling  resistance  co¬ 
efficient,  /V  is  given  by  pr  =  .012  +  1  x  10“^v  where  v  is 
the  axle  speed  in  INCHES/SEC.  Thus, 

(4a.  35)  ~  J^r  Rr  C.012  +  V) 

Or  alternately,  =  C-O‘2  1  )  C  FeR./S)(S  )  Rr 

Since  =  Cmt  ,  the  rolling  resistance  coefficients  are 

established  as; 

(4a. 36)  =  .012  Cmt  Rt  =  (.oi2)(sS3o)C2o.st)  =  235-o\t 

(4a.37)  Dvc  =  Xto.si)  =  40.3 

Figure  40  shows  the  friction  coefficient  for  a  tire  slid¬ 
ing  (i.e.  full  skid)  on  a  dry  concrete  runway  as  a  func¬ 
tion  of  velocity.  This  data  is  taken  from  Reference  3  and 
is  applicable  to  a  typical  runway  contaminated  with  rubber 
deposits  from  previous  airplane  operations.  Table  8 
below  lists  the  appropriate  coefficients  for  equation 
(4a. 22)  which  apply  for  dry  and  wet  runway  surface  condi¬ 
tions. 

Table  8  Runway  Friction  Characteristics 


SYMBOL 

UNITS 

WET 

CONCRETE 

DRY 

CONCRETE 

Ut, 

Ut2 

£t 

<< 

SEC/ IN 

SEC/ IN 

.050 

.180  _3 
,065  X  10  3 
1.0  X  lO’-^ 

.200 

.450  _3 
,065  X  10  o 
2.5  X  10 

Initial  Conditions 


All  initial  conditions,  except  wheel  and  tire  rotational 
speed,  will  be  set  to  zero.  From  the  airplane  system  at 
time  =  0,  Vp  =  2400  and  *  2.243.  Using  equation 
(4a. 16)  results  in: 

(4a.38)  Rr  =  -  5  5m  =  '  5  (2  245)  ^  22,fc7  in 

In  order  that  Vjj  be  zero,  equations  (4a.  18)  and  (4a.  19) 
show  that: 


.  ^ 

J 

i  . 

r 

0  eoo  4C0  MO  800  iOOO  i200 


Tire  Footprint  Velocity  (in /sec) 


Figure  40  Tire  Sliding  Friction  Coefficient 
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Wheel  and  Tire  System  (Flywheel)  Parameters 
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WHEEL  ^ND  TIRE  SYSTEM  (3  DEGREE) 


Figure  41  shows  the  components  of  the  wheel  and  tire 
system.  The  wheel  and  tire  system 


Figure  41  Components  of  the  Wheel  and  Tire  System 

for  the  3  degree  airplane  system  is  essentially  the  same  as 
for  the  flywheel  model.  The  Airplane  System  still  furn¬ 
ishes  the  tire  deflection  and  the  tire  vertical  load 

.  The  ground  speed,  however,  is  no  longer  furnished  by 
the  Airplane  System,  but  is  found  by  summing  forces  on  the 
tire,  wheel,  brake,  and  axle  mass.  The  horizontal  force 
exerted  on  the  axle  by  the  airplane  is  calculated  by  obtain¬ 
ing  the  translational  and  rotational  (0g)  gear 

positions  from  the  Airplane  System. 
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Figure  42  Tire  Horizontal  Model 

Referring  to  Figure  42  equation  (4a. 1)  in  the  flywheel 
system  changes  to 

(4b.  1)  =  C(^( 

Equations  (4b. 2)  throuj?h  {4b. 9)  are  listed  for  com¬ 
pleteness,  although  they  are  the  same  as  (4a. 2)  through 
(4a. 9)  . 

(4b  .2)  FVt  ~  EfT^Xfr  Xyj 

(4b.3)DTTCXv"X^j^)~  ETrCXrr“Xv} 

(4b. 4)  W&w  Xvi;  -  +  Efj- 

(4b.  5)  W7£  Xyr  ■  “Err  “  ^st 
(4b. 6)  E&gh  ~  Rk.0  Wr 
C4b.7)  E?rv  ”  ^ k:e  W j ^  ctrz. 0 t ^ 


Fif,ure  43  shows  the  rotational  model  of  the  wheel,  tire, 
axle,  and  lowe'^  strut  with  the  gear  rotation  0(i  added. 
Including  the  effect  of  there  follows: 

(4b.f.)Tj2]'  ■  Cj2f(^iv”  ^r)  '^y) 

(4b .  9)  "  t- CT (  " 0y  ) 

(4b.lO)Ts  =  Ci?.s(0s  DrsC^s'^g) 


Figure  43  Tire  Rotational  Model 


Equations  (4b. 11)  through  (4b. 16)  are  the  same  as  (4a. 11) 
through  (4a  .  16)  . 


11) 

Hr  =  Re 

; 

(^b. 

,12) 

Wr 

>  0 

Tre  "  j 

It 

W  T 

-  0 

1  5m  C- 

i  f 

Vl  j 

<L  C 
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(■4b.  13)  Wjy 0^  -  Ts 

(4b ,  14)  Wtw  ©IV  “  “'fizr  “Tg-p 

(4b.l5)WxT©r  '  ^rfttr  ■‘■TfeT  "  "^GR 

(4b.  16)  R-p  "  RpT  “  3  Sm  ~  ^  ^  TT  ~ ^ 

Because  the  ground  is  now  stationary  equation  (4a. 17) 
becomes 

(4b.  17)  Ves  ■=  Xrr 

The  remaining  equations  are  unchanged  except  for  noting 
that  the  outputs  Xwivj  and  for  the  Airplane 

System  ^are  obtained  by  renaming  Xrr  *  Thus  X^vm-  Xrr  and 
X^,VK,  =  Xrr  •  Continuing, 


(4b.  18)  Vr  -- 

^Rs  ■  Sr  Wr 

(4b.  19)  Wji  •= 

V\/'vv  “  Ws 

(4b. 20)  Tump 

C  *  ■“  ©  rtY  ^  ) 

(4b. 21)  Fst  ■= 

(4b. 22)  Ur  = 

Uti  ■^CUTi-ErVRje'"''^^  l:f 

Vr>0 

Vr  -  0 

-U-n  -  (U,T2~£rV^riS  )  e  id- 

Ve<  c 
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B.  Parameter  Evaluation 

The  parameter  values  for  this  system  are  essentially  the 
same  as  for  the  wheel  and  tire  system  that  corresponds  to 
the  flywheel  system.  One  difference  Is  the  value's  for 
C&  and  which  are  derived  in  the  Airplane  System  (3 
degree) .  These  values  are 

(4b. 23)  j  Cg-  =  ^CO,O00  Ib/i'n 

(  D<3.  -  7*3.6  lb  stc  / in 

Since  this  system  moves  V7ith  the  airplane  the  initial  con 
ditions  should  match  the  Airplane  model.  Thus 


(4b . 24)  1 

^  TTO 

r.  0.0 

in 

=  240C* 

in/  ^cc 

(4b. 25)  1 

1  X  vvo 

-  0.0 

in 

j 

^  X  vVO 

=  2400 

J  n  /  Stf  c 

From  equation  (4b,16) 

(-'lb. 26)  Rr  =  R^r  =  23,32  -.80  =  22. SZ  IN 

Thus  for  a  "spun  up"  tire,  we  have  from  equation  (18) 

('4b.27)  Wr  -  V(?s/Rr  =  2400/22,52  =  106,. 7  rad/sfic 


Then 


(4b. 28)  &ro  =  =  10G..7  rad/<ec 

Also  from  equation  (10),  choose  0so 

(4b. 29)  6so  =  =  .0329  racl/scc 
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Table  10  Wheel  and  Tire  System  (3  Degree)  Parameters 
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Table  10  (Contd) 
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YMBOL  TYPE  VALUE  '.KITS  DESCRIPTION 


FooCprint  Velocity 


4c.  TIRE  AND  VTflEEL  SYSTEM  (6  DEGREE) 


The  wheel  and  tire  system  for  the  six  degree  problem  is  the 
same  as  that  described  in  the  three  degree  system  except  for 
inclusion  of  the  lateral  mode.  Equations  (4c. 1)  through 
(4c. 17)  are  the  same  as  (4b. 1)  through  (4b. 17)  in  the  three 
degree  system. 


Mathematical  Description 

(4c. 1) 

F,.  = 

"^gCXax'Xw) 

Xw) 

(4c.2) 

Frr  ■ 

“  C7-r(Xrr  +  E_^^(Xy7-- 

-Xy) 

(4c. 3) 

Drr 

CXy~Xw)  -  E.-rT(XTt“Xy 

) 

(4c. 4) 

Wfiw 

X  w  "  ■*'  Ftt 

(4c. 5) 

We 

Xtt  ~  “  Frr  "  Fbt  ■*" 

(4c. 6) 

(4c. 7) 

f^RV 

<6-x  > 

(4c. 8) 

“  ^r)  F  gT  (0  vv 

(4c.  9) 

Dpr(6>'“^r)  ”  Eet^^w'^y) 

(4c. 10) 

Ts  = 

-  +  I>iis  (05  “ 

©J 

(4c. 11) 

Hr 

(4c. 12) 

“  /  5fvj  (  Dsc  -^^yicW-r} 

'li  \Vr  O 

is  Wt  =  o 

(■"  ■1'  Wt) 

if  W'r<o 

(4c, 13)  =  ^37  Is 

(4c,  14)  Wiw'^w"  ”^£7  “  ^sr 


(4c.l5) 

Wi  T  S  7  - 

Hr  f^T  ^ZT  ^EE 

(4c. 16) 

~  ■ 

-  j  Sm  “  ^ec  (  Xrr  ~ 

(4c. 17) 

Vrs  =  Xtt 

”  Xvviv)  also  XwM  “ 
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Equation  (4b. 18)  which  gives  the  relative  velocity  between 
the  footprint  and  the  ground  is  changed  to  account  for  the 
lateral  footprint  velocity  y„  . 

(4c.  18)  V,  =  V  Ve*  + 

Equations  (4b. 19)  and  (4b. 20)  are  unchanged. 

(4c.  19)  W0=Ww-Ws 

(4c.  20)  CHyCVes/Vny)  ) 

Now  is  the  relative  velocity  in  the  x  direction  so 
(4c. 21)  Vex  -  1/es  "  Rt  Vr 

Thus,  the  angle  which  defines  the  friction  force  direc¬ 
tion  as  shown  in  figure  45  is  given  by 

(4c. 22)  Ym/Vex) 


Figure  45  Footprint  Friction  Components 
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Figure  46  Wheel  and  Tire  System  (6  Degree)  Equation  Flow  Diagram 


Thus,  Fst  is  now  given  by 


(4c.  23)  FBr=  Ur  >  +  Dmy 


and 

(4c. 24)  Ur  =  I  Uti  +- (  Ut2  -  Ey  Ves)  e  ^r>^- 

(  O  »  f  V  J2  ^  o 

The  lateral  friction  force  F^y  is  given  by 

(4c.25)F;v-  Enmp  Ur  w'O.^rL'  ^/3r  ) 

The  airplane  produc^.  ’  a  lateral  position  Yax  of  the  axle. 

If  is  the  lateral  ocation  of  the  footprint,  then  a  lat¬ 
eral  force  is  proa  ced  and 

(4c.  26)  pj-p  =  C5T(yAx“Y^^)  ^  ~  Yivi) 

Finally,  forces  can  be  summed  laterally  on  the  footprint  to 
obtain 


(4c.  27)  VVr5  Y,„  =  Fjr  -  Fffv 
B,  Parameter  Evaluation 

Wheel  and  Tire  System  Parameters 

Most  of  the  parameter  values  were  derived  in  the  wheel  and 
tire  system  that  was  used  with  the  flywheel.  The  only  addi 
tion  is  the  evaluation  of  C^t  and  D^t  , 

From  reference  1  (p.lS) 

(4c. 28)  =  T;,  w  (P  t.Z4.PR)(  I  -  .7(5o/w)) 

Assuming  that  ==  2.245  in.,  then 

(4c.  29)  (2)(i8)(l.24)(l5c)(l-  ,7(2.24  5/18}) 

(4c.30)Csr-  Ib/tn 

Using  OJn  =  V  =  Va4foO/-(^725  =  298,^ 

Using  ^  =  .1  (from  page  50  of  Ref.  1)  results  in 

(4c. 31)  Dir  =  =  (.l)(c.+  4.o)/2584=  2.IU5 


Initial  Conditions 


The  only  difference  between  this  and  the  three  degree  sys¬ 
tem  is  evaluating  Ymo  and  V„o  .  It  should  be  that  Vmo  = 

Since  this  system  is  used  for  both  right  and  left  sides,  then 
Ymo  will  differ.  Assume  that  this  system  is  used  as  the 
right  wheel  and  tire.  Then,  Ymo  =  V^xc-  as  in  equation 
(4c. 56).  From  this  equation  at  Time  =  0,  (since  R  =  P  =  0), 
then 

(4Co32)  "  y©  +■  SsiAJ  ~  O  +  —  (oO  trj 

Similarly 

(4c. 33  Vmo  =  VftXR.  =  Vo  =  o.o  \n/s€c 
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SYMBOL  TYPE  VALUE  UNITS  DESCRIPTION 


5. 


VJHEEL  SPEED  SENSOR 


The  primary  input  parameter  to  an  electronic  antiskid 
control  circuit  is  an  airplane  wheel  speed  signal.  For 
conventional  control  circuitry  the  input  must  be  a  direct 
current  voltage.  The  wheel  speed  sensor  may  have  any  one 
of  several  forms  such  as  a  D.C.  tachometer  or  an  A.G. 
tachometer  with  variable  voltage  or  frequency  converted 
to  a  direct  current  voltage  by  suitable  electronic  cir¬ 
cuitry.  The  control  circuit  input  signal,  is  a 

function  of  the  wheel's  angular  velocity  relative  to  the 
ci-Ale  (tachometer  mount)  and  the  characteristics  of  any 
associated  electronic  circuitry  used  for  radio  interference 
suppression  and/or  for  conversion  of  A.C.  frequency  or 
voltage  signals  to  voltage.  To  provide  the  means  fer 

n.athematicslly  describing  the  control  circuit  input  signal 
for  a  variety  of  wheel  speed  sensors,  two  approaches  are 
taken.  The  first,  identi.^ied  as  Option  1,  is  applicable 
whenever  there  is  a  perceptible  phase  lag  between  actual 
wheel  speed  and  the  antiskid  circuit  inpu  as  is  generally 
the  case  where  A.C.  voltage  signals  are  converted  to  D.C. 
or  where  a  D.C.  tachom€^*er  is  driven  through  an  elastic 
coupling.  A  second  simpler  ntatbematical  description, 
called  Option  2,  is  provided  to  miriimlzr  computation 
dii^ficvlty  and  expense  inhere  no  significant  phase  lag 
exists . 

A.  Mathematical  Description 
Option  1 

Assume  that  a  D.C.  tachometer  generator  is  mounted  on  the 
axle  and  is  driven  by  the  wheel.  The  output  of  the  hypo¬ 
thetical  generator  is  assumed  to  be  applied  tc  a  linear 
force  motor  which  acts  upon  a  single  degree  of  freedom 
damped  spring  mass  system  as  sho^vn  on  Figure  47  The 
control  circuit  input  signal,  £<3  ,  is  proportional  to  the 
mass  displacement.  By  adjusting  the  relative  character¬ 
istics  of  the  linear  force  motor,  hypothetical  generator, 
spring,  mass  and  damper  a  mathematical  description  of  a 
wide  variety  of  wheel  speed  sensors  can  be  accommodated. 


146 


Figure  47  Wheel  Speed  Signal  System 

The  output  of  the  hypothetical  r/encrator  is  proport¬ 

ional  to  the  u-heel' s  an"ular  velocity  relative  to  the  axle, 
\Al8y  define^i  by  ccjuation  (5.1).  Angular  velocity ,  vVfl  , 
is  obtained  as  an  output  of  the  tire  and  v;heel  systen>. 


(5.1) 


Syvs  ^  G  ws  )/)! 6 


The  force  produced  by  the  linear  force  r:Otor,FvVf,  is  pro¬ 
portional  to  the  generator  output, as  defined  by 
eqviation  (5.2)  . 


(5.2) 


Fwi  -  tw/j 


The  hypothetical  mass  displacement , is  obtained  from 
equation  (5.3)  which  results  frorri  summing  forces  on  the 
hypothetical  masSjW’wj. 


(5.3) 


CtA/i' 


The  antiskid  circuit  wheel  speed  input  voltage  signal,  , 
is  proportional  to  the  hypothetical  mass  displacement,  X^s> 
as  defined  by  equation  (5.4). 


(5.4)  tS  -  IhNS  + 


In  equation  (5.^),  is  any  extraneous  "noise"  which 

might  be  present  due  to  the  operation  of  other  aircraft 
systems,  etc. 


The  equation  flow  diagram  is  shown  on  Figure  48. 
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Figure  48  Wheel  Speed  Sensor  Equation  Flow  Diagram  (Option  1) 


Option  2 


For  cases  where  the  wheel  speed  transducer  is  a  D.C.  tach¬ 
ometer,  or  equivalent,  driven  through  a  rigid  coupling 
(such  as  the  F-104  and  F-111)  there  is  usually  very  small 
difference  between  the  actual  wheel  speed  and  antiskid 
control  circuit  input  (i.e.,  very  low  phase  lag  or  atten¬ 
uation)  and  the  extra  mathematical  complication  incurred 
by  using  a  very  high  gain  second  order  equation  is  not 
justified.  For  these  cases  the  antiskid  control  circuit 
input  voltage  may  be  considered  proportional  to  the  wheel's 
angular  velocity  as  defined  by  equation  (5.5). 

-  C{4iQC  ^  Es// 

No  equation  flow  diagram  is  shown  for  Option  2. 

B.  Parameter  Evaluation 
Option  1 

Tne  objective  of  using  a  single  degree  of  freedom  damped 
spring  mass  system  to  describe  the  antiskid  control  cir¬ 
cuit  input  is  to  provide  a  mathematical  "tool"  whereby 
phase  lag  within  the  wheel  speed  sensor  device  can  be 
accounted  for.  Consequently,  the  values  for  mass,  spring 
rate  and  damping  coefficient  are  chosen  to  produce  the 
desired  effect  rather  than  to  describe  physical  devices. 

The  other  coefficients  are  chosen  to  achieve  compatibility 
with  the  control  circuit.  For  the  F-111  modulated  antiskid 
circuit  let  the  hypothetical  tachometer  coefficient  be  the 
same  as  for  the  actual  F-111  tachometer,  12  volts  per  thous¬ 
and  FPM.  Therefore: 

(5.6)  Tooo'CrJf  =  voij-s^c/^/to 

Let  the  force  motor  coefficient,  the  elastic  system  spring 
rate  and  output  voltage  coefficient  all  be  equal  to  unity 
so  that  for  steady  state  conditions  the  control  circuit 
input,  £<5,  is  the  tachometer  output.  Therefore: 

Cw5  '  i’O 

C.C.QV  -  J.Q  ^/OLT  / 
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Based  on  Information  furnished  by  the  Goodyear  Aerospace 
Corp.  the  component  characteristics  and  arrangement  which 
Is  usually  utilized  for  converting  A.C.  frequency  to  D.C. 
voltage  produces  about  30  degrees  (or  greater)  phase  lag 
at  “3  cps.  The  following  equations  from  reference  12 
describe  the  single  degree  of  freedom  system’s  behavior 
when  an  oscillatory  force  is  applied. 

X  ^ 

Xo 

(5.8) 


In  these  equations  ^  is  the  phase  angle  0^s/z\^\f>/s  UJ/i 

IS  the  damping  factor,  u)a  -  fCwi/W^vs  is  undamped 

natural  frequency,  uj  Is  the  frequency  of  applied  oscil¬ 
latory  loading,  and  X/Xa  is  the  magnification  factor. 

If  the  degree  of  attenuation  and  phase  angle  are  known  at 
a  particular  frequency,  the  undamped  natural  frequency  and 
damping  factor  are  established.  Assuming  two  percent 
attenuation  and  30  degree  phase  lag  at  5  cps,  the  equations 
above  give  an  undamped  natural  frequency  of  14.6  cps 
(91. P  and  a  damping  factor  of  0.746. 


For  an  undamped  natural  frequency  of  91. P  /s^c  and  a 

spring  rate,Cwi  ,  of  1.0  I  tot-// v  the  mass,  is  estab¬ 
lished  as  The  damping  coefficient, 

0^^  t  is  established  from  the  mass  and  damping  factor  as 

(5.9)  lh^rec//rJ 


Option  2 


For  use  with  the  F-111  modulated  antiskid  control  circuit, 
use  the  actual  F-111  tachometer  output  of  12  volts  D.C.  per 
1000  RPM.  Therefore: 


(5.10) 


\HOC 


/Zk  60 

/OOOK  2.n 


^oi~r 


^50 


For  use  with  the  on-off  antiskid  circuit  as  installed  on 
the  F-104  (and  B-58)  the  tachometer  output  is  20  volts  per 
1000  FPM.  To  make  the  on-off  circuit  compatible  with  the 
F-111  requires  that  the  difference  in  tire  size  (46.5  inch 
dia.  for  F-111  versus  22  inch  dia.  for  B-58)  also  be 
accounted  for.  Therefore,  for  the  on-off  antiskid  circuit 
use: 


(5.11) 


G\f/oc=  X 
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6a.  MODULATED  ..ANTISKID  CONTROL  CIRCUIT 


After  introduction  of  on-off  type  antiskid  systems,  it 
became  apparent  from  various  analyses  and  studies  of  test 
results  and  operational  performance  that  braking  effec¬ 
tiveness  could  be  increased  if  the  number  of  antiskid  cycles 
and  their  intensity  could  be  minimized.  To  minimize  anti¬ 
skid  cycli.ng  occurrences  and  intensity,  it  is  necessary  to 
control  the  amount  of  brake  torque  being  applied  such  that 
the  available  friction  torque  is  not  exceeded  for  as  much 
of  the  time  as  is  possible.  A  number  of  devices  utilizing 
various  principles  of  operation  have  been  used  for  this 
purpose.  These  devices  predominately  utilize  the  principle 
of  regulating  or  "modulating"  brake  pressure  to  keep  its 
value  as  near  as  possible  to  that  which  will  produce  a  skid. 
One  of  the  first  of  these  type  devices  is  a  hydraulic  pres¬ 
sure  modulator  comprised  of  an  orifice  and  accumulator 
installed  upstream  from  the  pilot's  metering  valve  and 
configured  such  that  repetitive  antiskid  cycling  causes  a 
temporary  reduction  in  pilot's  metered  pressure.  The 
Convair  Model  880  airplane's  Hytrol  MKI  antiskid  system 
with  hydraulic  modulation  is  a  typical  example  of  this  type 
installation. 

A  subsequent  development  was  the  Bendix  system  which  is 
used  on  Grumman  A6A  and  Lockheed  Cl4l  aircraft.  This 
system  combines  hydraulic  modulation  accomplished  within  the 
off-on  type  control  valve  with  two  levels  of  skid  detection, 
(i.e,,  brake  pressure  reduction  in  two  steps  controlled  by 
skid  intensity)*  Further  improvements  have  been  achieved  by 
utilizing  a  servo  type  pressure  regulating  valve  with  elec¬ 
tronic  control  to  achieve  a  wide  range  of  control  charac¬ 
teristics  and  better  accommodate  widely  varying  runway 
friction  conditions  encountered  dur5,ng  aircraft  operation. 

The  Goodyear  Adaptive  system  used  on  General  Dynamics  F-111 
aircraft  and  the  Hytrol  MK  II  system  used  on  McDonnell- 
Douglas  F4C  and  LTV  A7A  aircraft  are  examples  of  the  servo 
valve  type  systems.  Within  each  of  the  types  or  classes  of 
systems  there  are  a  number  of  variations  in  circuitry  and 
component  arrangement  depending  upon  the  aircraft  type, 
landing  gear  arrangement  and  configuration,  and  the  airplane's 
mission  requirements.  For  this  program  a  mathematical  model 
of  the  F-111  airplane's  Goodyear  Adaptive  Antiskid  Control 
Circuit  is  developed.  Models  for  other  type  circuits  can 
be  developed  using  similar  procedures. 
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Figure  49  is  a  block  diagram  showing  the  basic  functional 
elements  of  the  Goodyear  adaptive  antiskid  control  circuit 
as  used  on  the  F-lll  airplane  and  showing  the  relationship 
of  the  control  circuit  to  the  other  brake  system  components. 
During  antiskid  circuit  operation,  a  wheel  speed  signal  is 
provided  as  an  input  to  a  deceleration  detector.  Within  the 
deceleration  detector  the  wheel ‘s  deceleration  rate  is  com¬ 
puted  and  compared  to  a  threshold  value  provided  by  a  skid 
detection  threshold  circuit  element.  The  deceleration 
detector  produces  a  skid  signal  proportional  to  the  amount 
by  which  the  wheel's  deceleration  rate  exceeds  the  threshold 
value.  The  skid  signal  is  applied  to  a  valve  control 
amplifier  which  in  turn  produces  a  valve  control  signal 
proportional  to  the  input  skid  signal  plus  any  pressure  bias 
signal  which  might  exist.  The  valve  control  signal  is 
supplied  to  the  antiskid  control  valve  (a  servo  type  pres¬ 
sure  regulator)  for  brake  pressxire  control  and  to  a  modu¬ 
lation  circuit  element.  The  modulation  circuit  element 
interprets  the  valve  control  signal  and  provides  a  pressure 
bias  signal  to  the  valve  control  amplifier  and  a  threshold 
control  signal  to  the  skid  detection  threshold  circuit 
element.  The  wheel  speed  signal  is  also  supplied  to  the 
locked  wheel  prevention  circuit  elements  consisting  of  an 
airplane  speed  reference  and  a  wheel  soeed  comparison 
element.  When  the  airplane  speed  reference  indicates  that 
the  airplane's  speed  exceeds  "locked  wheel  arming  speed" 
(usually  2C  mph)  and  simultaneously  the  wheel  speed  is  less 
than  that  which  should  exist  for  a  slightly  lower  airplane 
speed  (usually  10  mph),  the  wheel  speed  comparison  circuit 
element  produces  a  skid  signal  sufficient  to  fully  release 
the  brake.  Locked  wheel  arming  speed  is  chosen  as  some 
reasonably  low  speed  below  which  a  locked  wheel  is  not 
particularly  detrimental.  The  locked  wheel  feature  is 
deactivated  below  locked  wheel  arming  speed  so  that  the 
airplane  can  be  brought  to  a  complete  stop.  The  aircraft 
circuit  also  incorporates  circuit  elements  for  failure 
detection,  automatic  cutoff  and  prevention  of  brake  appli¬ 
cation  prior  to  touchdown.  These  logic  type  functions  do 
net  affect  aircraft  stopping  performance  and  are  not 
included  in  this  analysis. 

A.  Modulated  Antiskid  Circuit  Mathematical  Description 

A  simplified  schematic  diagram  of  the  Goodyear  adaptive 
antiskid  circuit  for  one  wheel  as  used  on  F-lll  type 
aircraft  is  shown  on  Figure  50  .  This  circuit  accoin- 
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plishes  deceler •’tion  skid  control  as  previously  described 
in  the  control  circuit  functional  description  as  follows. 

An  input  voltage,  Ecj  ,  is  provided  by  a  wheel  driven  D.  C. 
tachometer  generator  (Geu).  £<5  charges  -he  deceleration 
detector,  capacitor, C/  ,  through  resistance  ^<5  and  diode  Oe 
during  wheel  spin-up.  For  normal  wheel  deceleration  rates, 
vith  no  incipient  skidding,  the  generator  voltage  will 
decrease  relatively  slowly  and  a  small  current  will  flow 
from  the  positive  side  of  Cf  through  £<;  ,  the  generator, 
and  transistor  Qt  to  the  negative  side  of  Ci, 

This  current  discharges  capacitor  O  and  causes  its 
voltage  to  closely  follow  ,  Transistor  Qi  is  the  skid 
detection  threshold  circuit  element.  Qi  is  a  current- 
limiting  de.'ice  that  offers  very  low  impedance  to  current 
below  its  chreshold  value  and  extremely  high  impedance  to 
any  current  above  that  threshold.  The  threshold  is  con¬ 
trolled  by  KZ  •  Diodes  Oz,  and  Oj  provided  bias  voltage 
for  the  operation  of  Qi .  When  an  incipient  skid  occurs , 
the  generator  voltage  decreases  rapidly  and  since  Qi  limits 
Lhe  discharge  current  into  Ci  ,  the  voltage  at  the  negative 
side  of  Cl  decreases  rnd  causes  current  to  flow  through  ^S/i^ 
Rb  ^  Qz^  and  /??  .  The  current  into  the  base  of  Ckt  is 
amplified  hyC?i,d?3  »  and  ^7,  (the  valve  control 

amplifier)  to  produce  a  voltage  across  Rv  ,  the  antiskid 
valve  coil.  Voltage  applied  to  the  antiskid  valve  causes 
brake  pressure  to  be  reduced  proportionally  and  thereby 
alleviate  the  incipient  skid.  Antiskid  valve  voltage  is 
feedback  to  the  amplifier  input  through  to  stabilize 
amplifier  gain  against  changes  due  to  temperature  and 
component  characteristic  variations. 

Antiskid  valve  voltage  pulses  are  transmitted  to  the  modu¬ 
lation  circuit  elements  through  capacitor  C4.  ,  Within  the 
modulation  circuit  element,  consisting  of  C4,  Ru  ,R'^, 

D7,  C3  ,  Cz  ,  D?  and  each  increase  in  voltage  to  the  valve 
produces  an  increase  in  the  charge  on  Cj  ,  Voltage  on  C5 
causes  CJ4  to  charge  £2  .  Since  £2  discharges  through  l?90,l?eA^ 
and  ,  the  voltage  onCi  provides  a  threshold  control 
signal  to^?,.  The  charge  on  Cj  ,  and  in  turn  on  Cz  .  is 
a  function  of  the  amplitude  and  frequency  of  valve  voltage 
pulses.  Voltage  onCz  is  also  applied  to  through 
and  Oi  to  provide  a  pressure  bias  signal  to  the  valve  con¬ 
trol  amplifier.  The  operation  of  the  modulation  circuit 
element  results  in  an  automatic  threshold  change  to  the  skid 
sensing  circuit  and  a  bias  to  the  valve  control  amplifier 
to  match  the  braking  conditions  being  encountered. 
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The  equations  listed  in  Tables  13,  14,  15,  16,  and  17 
describe  the  circuit's  operation  in  accordance  with  the 
equation  flow  diagram  shown  on  Figure  51  ,  The  assumptions 
and  procedures  used  to  develop  these  equations  are  described 
in  Appendix  I,  The  circuit  has  twelve  possible  modes  of 
operation  depending  upon  the  instantaneous  conditions  which 
exist.  The  conditions  which  define  the  circuit's  mode  of 
operation  at  a  particular  instant  are;  (1)  the  current, 

/Ic4,  into  capacitorfTf  is  either  positive,  negative  or 
zero,  (2)  the  valve  control  amplifier  is  operating  either 
in  the  cutoff  mode  or  in  the  amplification  mode,  and  (3) 
the  modulating  circuit  element  is  either  providing  a  pressure 
bias  signal  or  it  is  not.  Table  18  lists  the  twelve  cir¬ 
cuit  modes  resulting  from  combinations  of  the  above  circum¬ 
stances.  The  valve  amplifier  condition  is  indicated  by 
current  being  equal  to  or  less  than  C607  in  the  cutoff 

mode  and  being  greater  than  C607  in  the  amplification  mode. 
The  pressure  bias  signal  is  indicated  as  existing  when 
voltage  Va  is  greater  than  zero  and  as  not  existing  when  Vs 
is  equal  to  or  less  than  zero. 

The  circuit  condition  which  exists  at  a  particular  instant 
is  established  by  assuming  a  condition  and  using  the  equa¬ 
tions  applicable  to  the  assumed  condition  to  test  for  the 
validity  of  the  assumption.  For  instance,  circuit  condition 
number  1  assumes  that  current  ^^4-  is  positive;  therefore, 
equation  6a-N8-l  from  Table  16  must  indicate  a  positive 
value  of  for  the  assumption  to  be  correct*  If  so, 

then  the  equations  for  and  Va  are  similarly  tested. 

If  the  assumed  condition  is  correct,  the  applicable  equations 
are  used  to  compute  the  various  currents  and  voltages.  If 
the  assumed  condition  is  foucid  to  be  incorrect,  other 
conditions  are  successively  assumed  and  tested  until  the 
correct  condition  is  found. 

B.  Parameter  Evaluation 

Table  19  lists  the  parameters  defining  the  modulated  cir¬ 
cuit's  operation.  The  values  for  the  constants  are  computed 
from  various  circuit  element  characteristics  (resistance, 
capacitance,  etc.)  as  described  in  reference  13  and  in  the 
semiconductor  component  manufacturers  catalogs. 


Table  13  Modulated  Antiskid  Circuit  Equation  Summary 


Equation  No. 

Equation 

(6a-l) 

\4:(  =  y  Vci  cH: 

(6a-Al) 

Va  =  Ac\  C&>oB 

(6a-2) 

\Jc-ir  f  Vc2 

(6a-A2) 

\ic2r  ncz  Cbo? 

(6a-3) 

f  Vcs  dt 

(6a-A3) 

/ici  C&IO 

(6a-4) 

11:^.=  J  \/c<i-  di: 

(6a-4A) 

\lc4-  /?C4  Cbli 

(6a.5) 

~  F,s-  \/c.i 

(6a-6) 

Alms  Fofi.  \/f>  C  MD  EeC CC>iL 

•=lo  Mf-CGi£ 04  Es  F-Cbii, 

(6a-LW-l) 

hi  hi  -  hEAZ  *■  /li-ivr 

(6a-N3) 

/ic3 -iAc4  -  Va  C&I9-A8&4  Foiz/k4>o 

Fa-P  -  A0H4  Foh.  Iki  ^ o 

(6a-N4) 

£V  =  Ads  Oi-oi,  •- €407 

(6a-N5-n) 

See  Table  15 

(6a-N8-n) 

See  Table  16 

(6a-N10) 

(6a-Nll) 

At~  ~/?cal 

(6a-N14) 

/}cz  -  t-  y?c?4-l/c2  Foii  ^  0 

-  ]/cz  Cs/i  Fofi.  flcA-O 
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Table  13  Modulated  Antiskid  Circuit  Equation  Summary 
(Contd) 


Equation  No. 
(6a-Q-lC) 


(6a-Q4) 

(6a-RlO) 


(6a-VB-n) 

(6a-VQ4) 


Equation 

/4ill  -  Fofi- 

fbe-Vdi,- 

/Ie^4-  C6l^  fls(i4- 

fioS  -  tE&-'^c,)d6M-  ^6x1 

EoH  (Gs-^ci)  >  C6:ii ^/’^XQ 

■=  O  FOC  Lek-\fu)  £  C6ii/E^2^ 

See  Table  14 

Ci2z  - 

Esa.  i >  C^2iJCc,z . 
~Jo  -  l/ii)  E C6>2i/Ci,x 


A  *•  •- 


•W'^’C'-Wnv- 


Table  14  Pressure  Bias  Signal  Condition  Test  Equations 


Circuit 


Pressure  Bias  Signal  Test  Equation 


Condition,  n  (Equation  6a-VB-n) 


iS  =  -  C444  i£<s-Vci)  ~ C4i>s (\/cj i-\/c<)  ^€4^^ 

Ms-  \/cJiC/^  -C4^7 

ycxCitfi-CssaiEs'-^a] 

7  &  8  \/b-  )/cxC^  (£G-^ci) 

9  &  10  ^clCm  ~CSB2.i£<^-~)/ci] bkit^:^fCss4. 

11  &  12  ^8-  ^CxCm-Cs'SS(£s~^Ci)-CB8^l)(:z-*^c^  7CSB7 


Table  15  Siammary  of  Equations  for  Computing 
Current  Ads' 


Circuit 
Condition,  n 

Applicable  Equation  6a-N5-n  (See  Note) 

1  &  2 

AIdb  -  -  Ac4 

3  &  4 

ft  OS  -  Cm  C404  --Cscs  -/7c^ 

5  &  6 

ftps-  ftv/il 

7  &  8 

ftos-  ft//tr  C404  ^C4-os 

9  &  10 

flcs  -  ftVBt  Cd>0^  -  ftc4- 

11  &  12 

ft  OS  -  ftvAX  C6o^  C404  -  Csos  -  ftc^ 

Note:  For 

all  circuit  conditions  if  ftos<~0^  set  ftoS^O 

Table  16  Capacitor  C4  Current  Mode  Test  Equations 


Circuit 
Condition,  n 

Applicable  Equation  (6a-N8-n) 

1 

"  (^0 C4‘Ti>  -  L^C3fi/c'4)C47?  -  C47ii^ 

2 

/?C4  -  l/cz  - (^6 -i/oi  C479  - Ll/CStVc'4)  Cw  ■' Cif-yz 

3 

/?C4  -  '(^<S-/cy)C^i  -  i\/c^tV<.'4)^V7a  -  ^47/ 

4 

/Ic4  ^  Vcz  CHi  -  LVCi  tVc^C4n-C47S^ 

5-P 

5-N 

*  flc4  -  lE0-  IE/)  CaS^  r  (i/cx  r  li'i)  Cs’Ji  -  C^oi 

6-P 

-^/\c4  -  SCzCssz.-lS;~y^)Cs^3  \E4)C^S7  '^^^4  /fw 

6-N 

-  Ec4  ~~yCx^^^>  *^(E'4~k:JCs^S'4  ^LVCii^C^Cs'Vj  '^S:»7 

7-P 

•f  f\c4  '  -  -  i4’/jG7fc  -  0/c3  4  Vc^  C^3?  t  /3r^O 

7-N 

8-P 

^  ll\c4  ~  LL'z  EfiU  •(Eti-yc’/)Es'f3  -ll/d '//'44/f//fi 

8-N 

'  /Jl'4  C^i-f  •*'lEk  ^ci)  Cs'H  *'((62,*^^^)  '  //s'47 

9 

fic4  “  -C/dzrlE^JCs'^/  3'Es'7o 

10 

lh4  -  yC2.C^97.-‘(fd^-l/c/)E^s'9H  -Ll^x-*/c4)Cs'9fr^S94 

11 

flc4  -  '  -O/cxtl/c^C^iS’  7-  Csyy 

12 

flcA  '  Vc^Cm  -liL^I^)Cs'9/9'^a£>/ 

Table  17  Valve  Amplifier  Operating  Mode  Test  Equations 


Circuit 
Condition, n 

Applicable  Equation  (6a-Q2-n)  (See  Note) 

1 

/Jecsz.  =  - - /c/J  yc<t)C4^7  -t  C-^s'e 

2 

-  Vcz  C-f^/  -  ^Ycf)  C4S9  -  (i>^3  C4i>£>  -f 

3 

/lmz=  “(£<^-Vcl)C44i?~(i^^'f'^C4-)C4‘47  +0448 

4 

^£(3i=  VctC^Sb  -  lOs-]/ct)C‘449  +  +C4SZ 

5 

~  (£<5~^ci)  €£.7/  ^Cs'SZ 

6 

/}£Ql  -  Vcz  Csss  ~(£c'‘Va)CSS4  'Csss 

7 

ftEc^z  =  '  (£<s-\/a]CSzc.  +CS27 

8 

/laat  =  Vcz  CszF-fEs  'l/c/}Cszs  -Csso 

9 

Reat  -  ^(^^'Vci)Cs^s ’IVcz+Vc^Cs^c.  ^■C'sc? 

10 

/Jeqz  ~  l/ct  CSS>?-  (£s']/ci)CsZ9  'ti/cz4l/cV)CS7£>  '  CS7/ 

11 

flscit  =  'i£<S-\/ci\Cs?S~(^z+Vc4)  Cs7(,  i-cs?? 

12 

Z^EGiZ  ='  Vet  Cs?s  -  -  \Zct)CS79'‘(VczfVc4)Cs&0  -  Cset 

Note:  For  all  Circuit  Conditions  if 

Table  18  Modulated  Antiskid  Circuit  Conditions 


Circuit 

Condition 

Capacitor  C4 
Current  Mode 
(See  Note  1) 

Valve  Amplifier 
Operating  Mode 
(See  Note  2) 

Pressure  Bias 
Signal  Condition 
(See  Note  3) 

1 

/1C4.>0 

fl^QZ  ~  C^7 

Vb  -  C> 

2 

O 

fl£.eiz~C6o7 

Vs  >0 

3 

/^c^>  O 

A^QZ  >C6o7 

Vs  -<D 

4 

fic4>0 

Vs  >0 

5 

Aeqz  ^C^o7 

Vb 

6 

/lc4 

^C^d7 

Vs^o 

7 

/-ic4-  -O 

Asqz  >  C607 

Vs  -0 

8 

Ac^-o 

^  Ci:>07 

\(s>o 

9 

A^QZ  ^Ci^o7 

Vs  -0 

10 

Aesiz  ^  C^o7 

Vs^o 

11 

A  Eat  >  C^o7 

Vs  ^0 

12 

/lc-(-4  0 

AeGlX.  >  C607 

Vs  >0 

Notes : 

1.  Capacitor  is  charging  for  /Ica->0^  static 

for  /?C4-  =o  and  discharging  for  <  O, 


?  The  valve  amplifier  is  amplifying  for 

Aeqi  >€£>07  and  is  cutoff  for  /^bqz  ~C6o7. 

3.  A  pressure  bias  signal  exists  for  Vs>0 
and  does  not  exist  for 
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Modulated  Antiskid  Circuit  Equation  Flow  Diagram 


Type  Value  Units  Description  (See  Note) 
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Note;  All  equation  numbers  in  Description  are  prececieci  by  6a.  For  example, 
EQU  VB-n  means  equations  number  6a-VB-n. 
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Symbol  Type  I  Value  I  Units  Description  (See  Note  page  167) 
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Type  Value  Units  I  Description  (See  Note  page  167) 
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Symbol  Type  Value  Units  Description  (See  Note  page  167) 
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6b.  ON-OFF  ANTISKID  CONTROL  CIRCUIT 


Most  aircraft  on-off  type  antiskid  systems  operate  accord¬ 
ing  to  the  functional  block  diagram  shown  on  Figure  52 
The  various  functional  elements  may  be  electrical,  mechani¬ 
cal  or  a  combination  of  electrical  and  mechanical  devices. 
If  during  braking  the  brake  torque  applied  to  the  V7heel 
exceeds  the  amount  which  can  be  reacted  by  friction  at  the 
tire-ground  interface,  the  antiskid  system  operates  to  pre¬ 
vent  tire  skids  as  follows.  A  wheel  speed  signal  is  pro¬ 
vided  to  a  deceleration  detection  element  where  the  wheel*s 
deceleration  rate  is  computed  and  compared  to  a  threshold 
value  which  is  provided  by  a  skid  detection  threshold  ele¬ 
ment.  The  deceleration  detector  produces  a  skid  signal 
whenever  the  wheel's  deceleration  rate  exceeds  the  thres¬ 
hold  value.  The  wheel  speed  signal  is  also  supplied  to  a 
wheel  speed  reference  element  and  a  wheel  speed  ccxnparison 
element.  The  wheel  speed  reference  element  is  a  "memory” 
device  which  produces  a  "comparison  index."  The  "compari¬ 
son  index"  is  the  wheel's  initial  unbraked  speed  minus  an 
adjustment  to  account  for  the  aircraft's  deceleration.  The 
wheel  speed  comparison  element  compares  wheel  speed  to  the 
"comparison  index"  and  produces  a  skid  signal  whenever  the 
wheel  speed  is  less  than  the  "comparison  index."  The 
deceleration  detection  element  initiates  a  skid  signal  and 
the  wheel  speed  comparison  element  maintains  the  skid 
signal  until  the  wheel  has  regained  most  of  its  initial 
speed.  The  skid  signals  from  both  the  deceleration  detec¬ 
tion  element  and  the  wheel  speed  comparison  element  are 
transmitted  to  a  valve  control  element  which  acts  to  con¬ 
trol  the  antiskid  valve  such  that  the  brake  is  released 
when  a  skid  signal  exists  and  the  brake  is  applied  when  a 
skid  signal  does  not  exist. 

An  electrical  system  of  the  form  shown  on  Figure  53 
or  a  mechanical  device  as  shown  on  Figure  55  are  the  most 
couimou  means  used  for  implementing  the  on-off  antiskid 
system  function. 

EJecLricdl  On-Off  Antiskid  System 

Figure  53  is  a  schematic  diagram  of  the  Goodyear  elec¬ 
trical  on-off  antiskid  control  circuit  as  used  on  the  Lock¬ 
heed  F104  and  General  Dynamics  B-58  aircraft.  This  circuit 
accomplishes  on-ofl  an.iskid  control  according  to  the  pre¬ 
ceding  functional  description  as  follows:  The  wheel  speed 
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ff  Antiskid  Control  Functional  Block  Diagram 


(A)  Schematic  Diagram 


(B)  Mathematical  Identification  Showing  Transistor^? 
and  Diodes  in  Terms  of  their  Equivalent  Circuits 

Figure  53  Electrical  Oa-Off  Antiskid  Control  Circuit 
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indication  element,  a  wheel  driven  D.  C.  tachometer 
generator  (GEN) ,  supplies  an  input  voltage  EG  which  is 
proportional  to  wheel  speed.  EG  charges  capacitor  Cl 
through  R5,  Diode  D1  and  Resistance  RG  during  wheel  spin- 
up.  Capacitor  Cl  is  both  the  deceleration  detector  and 
the  wheel,  speed  reference  element.  For  normal  wheel  decel¬ 
eration  rates,  with  no  incipient  skidding,  the  generator 
voltage  decreases  relatively  slowly  and  a  small  current 
will  flow  from  the  positive  side  of  Cl  through  R5  and 
through  d2  and  Rl,  (GEN)  and  RG  to  the  negative  side  of 
Cl.  This  current  discharges  Cl  and  causes  its  voltage  to 
closely  follow  EG.  The  amplifier  comprised  of  R2,  Q1,  R3, 
R4  and  Q2  acts  as  the  skid  detection  deceleration  threshold 
element,  the  wheel  speed  comparison  element,  and,  in  con¬ 
junction  with  relay  Kl,  the  valve  control  element.  When 
an  incipient  skid  occurs,  the  generator  voltage  decreases 
rapidly.  R5  and  Rl  limit  the  discharge  current  flow  into 
Cl  so  that  the  voltage  at  the  positive  side  of  Cl  increases 
The  value  of  the  voltage  at  the  positive  side  of  Cl  is  pro¬ 
portional  to  wheel  deceleration  rate.  The  amplifier  charac 
teristics  are  set  so  that  when  the  voltage  at  the  positive 
side  of  Cl  is  a  value  Vsot  or  greater,  enough  current  flows 
into  the  base  of  Q1  to  cause  Q2  to  conduct  sufficiently 
for  relay  Kl  to  actuate.  When  relay  Kl  actuates  the  supply 
voltage  is  applied  to  the  antiskid  valve  coil  LV.  The 
voltage  across  the  antiskid  valve  coil,  EV,  is  equal  to  the 
supply  voltage.  Actuation  of  relay  Kl  also  causes  Rl  to 
be  disconnected  from  ground  so  that  the  resistance  in  the 
discharge  path  of  Cl  is  increased  to  aid  Its  action  as  a 
wheel  speed  reference.  Voltage  Vsor  is  the  skid  detection 
threshold.  More  modern  versions  of  this  circuit  utilize 
transistors  to  perform  the  function  of  relay  Kl;  however, 
their  operation  is  the  same.  Resistor  R5  is  the  speed 
reference  adjustment  element. 


A-1  Electrical  On-Off  Mathematical  Description 


The  mathematical  description  of  the  electrical  circuit’s 
operation  is  developed  from  Figure  53(b)  which  is  the 
schematic  from  Figure  53(a)  with  the  transistors  and 
diodes  shown  in  terms  of  their  equivalent  circuits  and 
the  appropriate  currents  and  voltages  identified. 

The  voltage  across  capacitor  Cl  is  defined  by; 


(6b-l-l) 


j  Vci 


where  (6b-l-Al)  Vci  =  AciC'?os 


'705  *  ^ /Ci 


Current  ACl  is  established  by  summing  currents  at  node 
^  as : 


(6b-l-Nl) 


/Ici  =  /lot  +  *  /toz 


Using  Ohm*s  law  and  summing  voltages  around  the  loop  of 
which  RDl  is  a  part,  current  ADI  is  established  as: 


(6b.l-Rl) 


/?oi  =  E<5  -  l^/  '  t.Q\ ^ Ro\  >0 

■  ^  (£e; ^ 


To  combine  constants,  write  equation  (6b-l-Rl)  as: 

/Jos  =  (£<5 -Vci)  C'70^  -  C'?o^ 

Noting  that  because  of  diode  Dl,  P\x>\  is  restricted  to 
positive  values  only. 

In  a  similar  manner,  using  Ohm’s  law  and  summing  voltages 
around  the  loop  containing  R4,  current  /^fl<^is  established 
as : 


(6b-l-V2) 


Slimming  currents  at  node  gives  : 
(6b-l-N2)  /!dz= /^eai  * /Iri 

By  Ohm's  law  the  voltage  across  RD2  is 
(6b-l-v3)  \/o2  =  ^oz  R.OZ 


For  the  case  where  no  skid  signal  exists  and  relay  K1  is 
not  actuated,  R1  is  connected  to  ground  and  a  current 
ARl  may  flow.  Using  Ohm's  law  and  by  summing  voltages 
around  the  loop  Rl,  D2,  Cl  and  (GEN)  is  established 
as  ; 

(6b-l-V4)  '  tZoz  -  \^z  J 


By  substituting  equation  (N2)  into  equation  (Nl)  current 
/\ci  is  established  as: 

(6b-l-Nl)  '  f\c\-  /foi  ^ 


Since  the  variables  EG  and  VCl  are  always  used  in  the 
form  of  their  difference,  define  the  difference  as: 

(6b-i-3)  (e<j^i/c,')  =  £5-14, 

By  substituting  (6b-l-V3)  and  (6b-l-N2)  into  (6b-l-V4), 
(6v-1-v4)  '  -toz)  £jsa  \ 

By  summing  currents  at  node  (^,  current  AEQl  is 
(6b-l-N3)  ftizri\=  fisaz  ^ 

By  summing  voltages  around  the  loop  containing  R3  and  the 
base  and  emitter  of  Q2, 


(6b-l-V5) 


O  =  Vk3  ' 


Note:  For  Q2  the  base-emitted  junction 

potential  has  been  omitted  to  reduce 
mathematical  complexity.  This  is  jus¬ 
tified  beca'vise  whether  or  not  Q2  is 
conducting  has  negligible  effect  on 
current  /?cj. 

By  substituting  (6b-l-N3)  along  with  the  Ohm's  law 
expressions  = 

and  the  transistor  characteristic 
into  (6b-l-V5)  and  solving  for  * 


(6b-l-V5)  ' 


z  ' 


By  substituting  (6b-l-V5)  *  and  (6b-l-N3)  into  the  Ohm's 
law  expression  Vgg  ^ 


(6b-l-2) 


By  summing  voltages  around  the  loop  R3,  REQl,  EQl,  RBQl, 
R2,  RD2,  Cl  and  (GEN) 

(6b-l-V6)  O'  V^3  *  ^ 


By  substituting  (6b-l-2)  and  (6b-l-V4)  along  with  the 
Ohm's  law  expres  sions  V^ai  ~  ^  ^si  -  Qnd 

\/j^2  -  transistor  characteristic 

into  (6b-l-V6)  and  solving  tox  ^60.1  \ 

(6b-l-V6)  ^  iVc/- 1:^)  Cvo/  '  ^?£>z. 

--  O 

For  the  case  where  relay  Kl  is  actuated  and  R1  is  dis¬ 
connected  from  ground  the  same  substitution  is  made  except 
that  ^02^  is  used  in  place  of  equation  (6b-l-v4y. 

For  the  actual  circuit  components  used  on  the  aircraft, 
the  resulting  equation  has  coefficients  that  are  negli¬ 
gibly  different  from  (6b-l-V6)  *;  therefore,  equation 


(6b-l-v6)  *  will  be  used  for  both  cases. 

The  value  of  /?^^i  which  causes  relay  K1  to  be  actuated 
is  defined  as,  C700,  the  skid  detection  threshold  current. 
From  this  definition  and  equation  (6b-l-V4) * 

/}^i  -  C709  "  ^  C^// 

/^ea\  ^  C^oa 

-  O  ^  C 

is  not  actuated  EV  ■  0,  when  relay  K1  is 
VS ;  therefore, 

£^U  -  o  /^OA  A3^\^C?£>o 

'  Vs  ^  C  7oo 

The  equation  flow  diagram  for  the  electrical  on-off  con¬ 
trol  circuit  is  shown  on  Figure  54  . 


(6b-l-V4-l) 

When  relay  K1 
actuated  EV  “ 

(6b-l-EV) 


B-1  Electrical  On-Off  Parameter  Evaluation 


Table  20  lists  the  parameters  and  their  values  as 
applicable  for  the  General  D3niainics  B-58  control  circuit. 
(The  same  circuit  is  used  on  the  Lockheed  F-i04,) 


Figure  54  Electrical  On-Off  Circuit  Equation  Flow  Diagram 
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Mechanical  On-Off  Antiskid  Device 


Figure  55  is  a  schematic  drawing  showing  the  operating 
principles  of  a  coomonly  used  mechanical  on-off  antiskid 
device  of  which  the  Hydroaire  Hytrol  Mk  i  and  Dunlop 
Maxaret  units  are  typical  examples.  The  device  operates  as 
follows.  The  rotor  (the  wheel  speed  indication  element)  is 
connected  to  the  aircraft  wheel  by  some  positive  means  such 
as  a  direct  connection  or  gear  train,  etc.,  so  that  the 
rotor's  angular  velocity  is  a  constant  ratio  of  aircraft 
wheel  speed.  During  spinup  the  motion  of  the  rotor  ic 
transmitted  through  the  clutch  to  the  drum.  The  clutch  is 
configured  such  that  it  is  self-energiaing  fo..-  rotation  in 
the  direction  of  wheel  rotation  associated  with  forward  air¬ 
plane  motion,  shown  here  as  counterclockwise.  Stop  (DSl) 
on  the  drum  engages  stop  (FSl)  on  the  flywheel,  thereby 
transmitting  torque  to  cause  the  velocity  of  the  flywheel 
to  be  the  same  as  the  drum.  The  flywheel  and  the  drum  are 
connected  by  spring  (CA)and  damper  (DA),  As  the  aircraft 
wheel  and  rotor  decelerate,  a  clockwise  torque  is  trans¬ 
mitted  through  the  clutch  to  the  drum  and  from  the  drum 
through  spring  (CA)  and  damper  (DA)  to  the  flywheel.  The 
amount  of  this  torque  is  proportional  to  the  product  of  the 
rotor's  deceleration  rate  and  the  flywheel's  inertia.  The 
torque  compresses  spring  (CA)  so  that  the  flywheel  moves 
counterclockwise  with  respect  to  the  drum.  For  steady 
airplane  wheel  deceleration  the  amount  of  relative  motion 
between  the  flywheel  and  drum  is  proportional  to  the  decel¬ 
eration  rate.  A  suitable  mef^hanism  (usually  a  set  of  elec¬ 
trical  contact  points  or  a  cam  device)  connected  between 
the  flywheel  and  drum  causes  a  valve  to  be  actuated  so  that 
brake  pressure  is  relieved  whenever  a  pre-established  amount 
of  relative  motion  occurs.  The  clutch  is  also  configured 
so  that  when  the  torque  from  the  rotor  to  the  drum  is 
clockwise,  the  torque  capacity  is  limited  to  some  slightly 
greater  amount  than  that  required  to  initiate  brake  release. 
If  the  rotor  experiences  greater  deceleration  than  that 
required  to  initiate  brake  release,  the  clutch  slips  and 
allows  the  drum  and  flywheel  to  overrun  the  rotor.  The 
flywheel's  inertia  reacted  by  the  drag  of  the  clutch  main¬ 
tains  a  torque  on  spring  (CA)  so  that  the  rel:itive  motion 
between  the  drum  and  flywheel  (skid  signal)  is  sustained 
until  the  flywheel's  kinetic  energy  is  dissipated  or 
until  the  rotor  has  regained  sufficient  speed  to  eliminate 
clutch  slippage.  For  this  device  the  flywheel's  inertia 
causing  displacement  of  spring  (CA)  and  damper  (DA)  is  the 


deceleration  detector  element,  the  clutch’s  overrunning 
drag  torque  on  the  drum  is  the  reference  adjustment 
element,  the  clutch  is  the  wheel  speed  comparison  ele¬ 
ment  and  the  rotational  kinetic  energy  of  the  flywheel  is 
the  wheel  speed  reference  element, 

A-2  Mechanical  On-Off  Mathematical  Description 

The  mathematical  description  of  the  mechanical  on-off  anti¬ 
skid  device  is  developed  by  referring  to  figure  55  f  which 
defines  the  applicable  parameters  and  shows  flywheel  stop 
(FSl)  represented  by  a  spring-damper  system.  Also,  a 
spring-damper  system  is  added  between  the  rotor  and  clutch 
carrier  to  represent  the  small  motion  which  actually  occurs 
during  clutch  operation. 

At  flywheel  stop  (FSl)  there  is  a  torque,  TS,  which  is 
exerted  on  the  flywheel  by  the  drum,  if  drum  stop  (DSl) 
is  in  contact  with  FSl.  If  the  mass  of  FSl  is  considered 
small  in  compariscn  to  the  stop  spring  (CS)  and  stop  dam¬ 
per  (DS)  then,  setting  the  sum  of  torques  on  FSl  at  zero: 

(6b-2-l)  Ts  =  Cs  (K-y)  -DsY 


Where  Csiyo-/)  is  the  scop  spring  torque,  (-Osy)  is  the 
stop  damper  torque  and  Yo  is  the  free  length  of  spring  CS. 

Since  TS  results  irom  a  contact  force,  it  cannot  be  less 

than  zero;  therefore,  if  Y ^  L&o-Of)  is  less  than 

then  TS  ■  0.  Rewritirig  equation  (1)  solving  for  ^  gives; 

(6b-2-2)  Cs  (Yo-y) -Ts/ Os 

Y  is  then  established  by: 


y  as  computed  from  (■6b-2-2)  and  (6b-2-3)  is  compared  to 
^fs-C&o-Bp)  to  establish  TS.  If  TS  is  other  than  zero,  it  is 
computed  from  (6b-2-l)  using  Y  ^  &Fs-(^o-eF)  ANfi  -(Sd-Sf). 
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Substituting  the  above  expressions  for  and  f  into 
(6b-2-l)  gives: 

(6b-2-l-l)  TS  -  Cs[/o-<^Fs  +  (^p-'^r^  ^Os 

Foa.  I  leo-^p]  -BfsJ  -o 
-  O  Fo/Z.  T-^' /■  C^D  “B'f')  -BfsJ  ^  o 

Siamning  torques  on  the  flywheel  gives: 

(6b-2-4)  ^  ^  Ts  (Bo-Bf) 

Summing  torques  on  the  drum,  gives : 

(6b-2-5)  ©0  =  C"Ts  ~Dfl(SD-SF)^~^c^  j W'^o 

Where  Ic  is  the  clutch  torque. 

Subtracting  (6b-2-4)  from  (6b-2-5)  results  in: 

(6b-2-6) 

^  7c/Wo 

By  integrating  (6b-2-6)  twice,  ('(&□- ^r) 
and  established  as 

(6b*2-7)  (&o*Op)  ~  J" i^o-  Of)  eft 

(6b-2-8)  [Oo-^f)  =  J(eD-eF)d± 

Substituting  values  for  and  I&q-Of)  computed  from 

(6b-2-7)  and  (6b-2-8)  into  equation  (6b-2-4)  and  integrating 
once  establishes  0f=  as  follows: 

(6b-2-9)  Sf 

Combining  the  results  from  (6b-2-7)  and  (6b-2-9)  establishes 
Oti  as : 

(6b-2-10)  ©0=  (©o-©f)^Op 
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The  clutch  will  now  be  examined. 

The  torque  exerted  on  the  clutch  carrier  by  the  rotor, Tc  , 
is  defined  by: 

(6b-2-ll)  Tc  =  Cc  (  +  Oc 

If,  as  for  the  flywheel  stop,  it  is  assumed  that  the  clutch 
carrier  inertia  is  negligibly  small,  the  torque  between  the 
clutch  and  the  drum  equals  the  torque  between  the  rotor 
and  the  clutch  carrier.  In  this  case  equation  (6b-2-ll)  may 
be  solved  for  (©V-©c)and  by  integrating  once  is 

obtained : 

(6b-2-12)  (S'/€-<S»c)=  J' 

Where  is  obtained  from  the  following  version  of 

(6b-2-ll) 

(6b-2-ll-l)  ?  [Tc  ^  Cc  ^  Dc, 

It  follows  that: 

(6b-2-13)  Sc 

If  the  clutch  is  configured  so  that  there  is  no  slipping 
for  counterclockwise  torque  on  the  drum,  <9c  must  equal  ©o 
and  any  difference  between  Bk  and  ©'o  must  be  relative 
velocity  between  the  clutch  carrier  and  the  rotor  (i.e, 

&■  Sc)  •  If  So  is  substituted  for  ©c  in  equation  (6b-2-ll) 
the  .  the  resulting  equation  can  be  used  to  compute  the 
UOtq\e  required  to  force  ©c  to  be  equal  to  ©d  .  There¬ 
fore,  making  this  substitution, 

(6b-2-ll-2)  Tc^  Cc  +  Dc(©«-an) 

Equation  (6b-2-ll-2)  adequately  describes  the  component 
of  clutch  torque  due  to  relative  velocity;  however,  the 
component  due  to  relative  displacement  is  not  satisfactor¬ 
ily  described  because  the  torque  direction  is  independent 
of  relative  position.  To  compute  the  clutch  torque  for 
all  conditions,  equation  (6b -2 -11 -2)  will  be  modified  and 
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a  procedure  for  establishing  the  clutch  condition  will  be 
defined.  The  clutch  condition  is  established  by  the  tor¬ 
que  direction.  The  torque  direction  is  determined  by  exam¬ 
ining  the  direction  the  drum  is  attempting  to  move  relative 
to  the  clutch.  The  direction  of  the  drum's  attempted  move¬ 
ment  relative  to  the  clutch  is  established  by  comparing  the 
drum  velocity,  Oo  *  to  the  velocity, ,  of  a  hypothetical 
or  "index”  clutch.  The  "index"  clutch  will  be  permitted  to 
have  slight  slippage  on  the  drum  for  counterclockwise 
torque  so  that  there  is  a  preceivable  circumstance  to 
indicate  torque  direction.  To  describe  the  "index"  clutch 
motion  relative  to  the  rotor,  equation  (6b-2-ll-l)  is 
modified  by  substituting  ©cv<  and  ©ch  for  Oc  and  ©c  as 
follows : 

(6b-2-ll-lrl)  (eK-ecn)  =  [To  -Ccfe^-ecJ)]  / Oc 


The  clutch  torque,  Tc  ,  is  defined  by  equation  (6b-2-ll-3). 
f©-«-©cH)is  obtained  from  equation  (6t-2-12)  and  <9c^  is 
then  established  from  equation  (6b-2-13),  noting  that  in 
each  case  Sen  and  Qch  are  used  in  place  of  Oc.  and  . 

The  clutch  condition  is  established  by  the  difference 
between  Och  and  ©V  as  follows: 


For  >0  Clutch  torque  is  positive  on  the 

drum  (clutch  attempting  to  have 
positive  velocity  with  respect 
to  drum) 

For  (OcH'S^  =0  Clutch  is  not  attempting  to 

move  relative  to  drum 


For  Clutch  torque  is  negative  on 

the  drum.  (Drum  is  attempting 
to  have  positive  velocity  with 
respect  to  clutch). 


Now  that  the  clutch  condition  is  defined,  equation  (6b-2-ll- 
2)  is  modified  so  that  the  torque  direction  is  established 
by  the  direction  of  relative  velocity  between  the  drum  and 
the  clutch  as  follows: 
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(6b"2-ll"3)  Tc  ~  Gc  C  \  Cc  ( ^R''©cn'll  ^  Dc  (Sn  “^p) 

I  >‘C.7SO 

-  Cr^o 

/^<|Gc<ecrt-(9t.>/cc  ('i%-(3:/*)|+Cc  (©«.©£,)  ^Cvs-o 

The  function  <ic <©cit-6o>  is  defined  as  follows: 

(6b-2-14)  6‘c<©ck-©o>  -  ^lo  Fo/Z.  >0 

-  O  Pp>€  (©cw~©o) 

~  ~I.C>  (©C,4  "^d)  O 

The  constant  C750  is  the  value  of  clutch  drag  torque  when 
the  drum  is  overrunning  the  clutch. 

The  amount  of  relative  motion  between  the  fl3rwheel  and  drum 
(©o-/9p)  is  the  skid  signal.  To  be  compatible  with  the 
electrical  antiskid  control  circuits,  assume  the  skid 
signal  is  produced  by  a  set  of  electrical  contact  points; 
therefore, 

(6b-2-15)  /j  paR  -  C7S/ 

-  o  L&-0  ^  C9S/ 

C751  is  the  skid  detection  threshold  value  of 
Also,  fcr  compatibility  with  the  other  parts  of  the 
analysis,  let  the  input  be  derived  from  the  wheel  speed 
sensor  output,  EG,  as  follows; 

(6b-2-16)  Cts'z  Eq 

C752  is  tlie  conversion  coefficient.  The  equation  flow 
diagram  for  the  mechanical  on-off  antiskid  device  is 
shown  on  Figure  56  . 

B-2  Mechanical  On-Off  Parameter  Evaluation 

No  parameter  evaluation  has  been  accomplished  for  the 
mechanical  on-off  device  because  it  is  not  applicable  to 
the  aircraft  being  considered. 
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Figure  56  Mechanical  On-Off  Device  Equation  Flow  Diagram 


7.  ANTISKID  CONTROL  VALVE 


Aircraft  antiskid  control  systems  typically  utilize  a  two- 
stage  electrically  operated  pressure  control  valve.  The 
first  stage  contains  an  electro-mechanical  device  such  as 
a  torque  motor,  solenoid  or  linear  force  motor  which  posi¬ 
tions  a  hydraulic  flow  regulating  element  (flapper,  nozzle 
or  spool)  such  that  a  control  pressure  is  produced.  The 
control  pressure  is  a  function  of  the  valve  input  pressure 
and  the  electrical  input  signal.  The  first  stage  control 
pressure  is  applied  to  the  second  stage  hydraulic  flow 
controlling  power  spool.  The  second  stage  spool  is  posi¬ 
tioned  by  forces  produced  by  the  control  pressure  and  valve 
output  pressure  in  a  manner  such  that  output  pressure  is 
controlled  in  proportion  to  the  first  stage  control  pressure, 

A.  Mathematical  Description 

First  Stage 

The  function  of  the  first  stage  can  be  described  mathe¬ 
matically  by  considering  the  control  pressure  producing 
element  to  be  a  single  degree  of  freedom  damped  spring 
mass  system  as  shown  in  Figure  57  acted  upon  by  a  force, 

Fcv  ,  proportional  to  the  electrical  input  signal, 

=■  Cscv  Ev 


Figure  57  First  Stage  Spring  Mass  System 
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The  first  stage  control  pressure,  Psc  ,  is  defined  as  a 
function  of  the  mass  position,  Xsc  ,  according  to  Figure  58  . 
Xsc  is  established  by  equation  (7.2)  which  results  from 
summing  forces  on  the  first  stage  mass,  l/Vsc. 


(7.2) 


-£i£ 


Xsc 


Osc  ; 
Wsc 


Figure  58  First 


(7.3)  . 

(7.4)  , 


Stage  Control  Pressure  -  Mass  Position 
Relationship 

7?5c  LPmv  -  ^  Pcif/e 

[  X.O  //r  )(sc  —  A  scm 


Xsc^  ~  A  Sc  /fi-  ^  Xsc  c.  /Xsc /I 

Asc ^  -  7C 

O  ip  Xsc  ^  Xscjc 
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Second  Stage 


The  physical  arrangement  of  the  F-111  antiskid  valve  second 
stage  is  shown  schematically  in  Figure  59  .  Most  other 
antiskid  valves  have  the  same  operating  principles. 
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Figure  59  Antiskid  Valve  Second  Stage 


As  described  in  the  hydraulic  system,  the  metering  valve 
output  pressure ,  ,  is  supplied  to  the  antiskid  control 

valve  second  stage  inlet.  When  the  second  stage  spool  is 
displaced  in  a  positive  direction  a  fluid  passage  opens 
permitting  hydraulic  flow  from  the  metering  valve  to  the 
antiskid  valve  outlet  cavity.  When  the  second  stage  spool 
is  displaced  in  a  negative  direction  a  fluid  passage  opens 
permitting  hydraulic  flow  from  the  outlet  cavity  to  return. 
Therefore,  the  second  stage  spool  position  defines  the 
hydraulic  flow  areas.  The  second  stage  spool  position, 

>  is  esLo’^lished  by  equation  (7.5)  which  results  from 
summing  forces  on  the  spool  mass,l4/<;^.  Figure  60  shows  a 
schematic  of  a  single  degree  of  freedom  damped  spring  mass 
system  representing  the  antiskid  valve  second  stage  spool. 
Springs  ,CcK^  ,  and  dampers ,  ,  are  stops  representing  the 

spool's  longitudinal  restraint  caused  by  its  contact  with 
the  valve  body.  The  forces  acting  on  the  spool  are  the 
positioning  spring  force,  damping  force,  stop  spring  and 
damper  forces  and  forces  due  to  outlet  cavity  pressure ,  A**/ , 
and  control  chamber  pressure ,Ajc^ . 
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Figure  60  Second  Stage  Spool  Forces 
Snniming  forces  on  the  second  stage  spool  give.: 

/  7  In  n  \  I  w  i,  \ 


Ccv  Oc^il 


(^*5)  =  (P:>cv  -Pc^/}  iKc\i/> 'Xcv)  " 


+  ££i^  Fc^s^ 

^cv 


V^CI/ 


(7*6)  - 


ir  Xcv/—  “'J’c 


(7.7)  r  (0  y^=  Xct"  - 

=  ^  ,(=  Ae/ 

The  control  chamber  pressure,  ,  is  established  by  the 

first  stage  control  presure  causing  flow  through  the  control 
orifice , ,  and  the  control  chamber  volume,  ]/scv  ^  as 
follows : 

(-J  f  -  n  *v  \ 


(^.P)  -  JpL 

Vsc^  ^ 

(7.9)  i  Xc\/c  'i' K  c\^ 

(7.10)  C^sc  =  Xf r^c  y  PjCi/P^ 


The  function  ^-S  defined  in  the  hydraulic  system. 

The  hydraulic  system  contains  provision  for  leakage  flow 
associated  with  first  stage  pressure  regulation  and  spool 
fit.  Since  these  small  flows  have  no  effect  on  the  valve *s 
performance  in  the  case  under  consideration,  they  have  not 
been  computed.  Therefore,  the  following  equations  apply: 

(7.11)  Qc\//--0 

(7.12) 

(7.13)  Qevs-O 

The  Control  Valve  Equation  Flow  Diagram  is  shown  on  Figure  61 


Figure  61  Antiskid  Control  Valve  Equation  Flow  Diagram 


B.  Parameter  Evaluation 

The  parameters  used  for  describing  the  antiskid  valve *s 
first  stage  behavior  are  established  from  measured  fre¬ 
quency  response  performance  characteristics  along  with  some 
features  of  its  physical  construction.  Frequency  response 
test  results  from  the  F-111  antiskid  valve  show  25  degrees 
phase  lag  at  5  cps.  From  various  experiments  it  is  known 
that  the  first  stage  accounts  for  most  of  the  phase  lag. 

The  F-111  valve’s  approximate  700  cps  undamped  natural 
frequency  is  quite  high  compared  to  a  more  usual  100  cps 
value.  Since  antiskid  cperatior.  is  generally  10  cps  or 
less  and  since  the  low  frequeno  phase  lag  can  be  accu¬ 
rately  described  with  a  lower  natural  frequency  system,  an 
undamped  natural  frequency  of  100  cps  will  be  used  to  mini¬ 
mize  computation  difficulty. 

Coefficients  C^c\/  and  Csc.  are  set  arbitrarily  so  that  static 
values  of  Xsc  will  be  compatible  with  values  of  Xse.m  and 
(which  are  also  arbitrarily  chosen)  and  the  proper 
valve  voltage  -  output  pressure  relationship  is  achieved. 

In  this  example  the  following  values  are  assigned: 

Csc^/ ^  /,Z  Ib9/Vo{t 
C.  sc.  ~  /,C 

(7.14)  'Ascm  -  ^*4  //je/ftfs 
Ksca.  ' 

For  100  cps  (628  /S£c  )  undamped  natural  frequency  and 

Csc  - 1.0  ibf////  I  ®3ss  Wse.  is  computed  from  LuJaY'^  Csc/\Alsc 

(7.15)  ^SC  -  J6(- 

Using  the  equations  relating  natural  frequency  and  phase 
angle  listed  in  the  wheel  speed  sensor  parameter  evaluation 
and  assuming  the  first  stage  has  20  degrees  phase  lag  at 
5  cps,  the  damping  factor  is  established  as  5^.  ^5  .  For 

the  values  of  and  Csc.  above  this  damping  factor 

results  in  : 


(7.16) 


Dsc.  =  /I.Cf  rClo  ^  lh(~  s^c 


The  area,y^c//»,  the  stop  clearances , and  ,  and  the 

liiass  of  the  second  stage  spool,  ^  are  computed  from  the 
spool*s  physical  dimensions  as  shown  on  the  valve  drawing. 

-  Ot  OS'  //y  ^ 

WcV  -  ^/,S K  sec^//^ 

(7 .17)  Jcv/f  =  o.oS 

SavA  ^  /ycM 

The  positioning  spring  rate,<rc'k’,  and  spool  damping  coc^rfi- 
cient  were  established  based  on  the  valve *s  transient 
response  characteristic  where  it  was  observed  that  a  50  cps 
about  .5  critically  damped  transient  pressure  oscillation 
appeared.  From  this  observation 

Ccv  -  4'-^ 

(7.18)  Dc\/ ^  /S‘0  A  /Jbf- sea J 

The  stop  spring,  ,  and  damper,  Oci/i  ,  characteristics 
are  arbitrarily  chosen  to  be  as  high  as  possible  within 
computation  capability. 

CtV5  •»  Sooo 

(7.19)  0cv6  ^  /iFse^//^ 

The  control  orifice  sreai/^cv'i;  ,  and  control  chamber  length, 
j  are  established  by  Lhe  valves  physical  dimensions  as 
shown  on  the  valve  drawing. 

O.Od?  //V^ 

(7.20)  OJ  /^cM 

The  control  chamber  fluid  bulk  modulus  is  that  of  MIL-H-5606 
hydraulic  fluid  as  used  in  the  hydraulic  system. 

The  undeflected  positioning  spring  length  was  computed 
assuming  it  produced  approximately  the  same  force  on  the 
valve  spool  as  25  PSI  pressure  differential. 

(7.21)  Xc^fi  -  O^Z  /a/cM 

Table  21  lists  the  paramete'  and  their  values  which  are 
applicable  to  thti  F-111  ar.u  ^^id  control  valve. 


Table  21  Antiskid  Control  Valve  Paramet«irs 
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8.  HORIZONTAL  TAIL  CONTROL 

In  the  3  degree  and  6  degree  airplane  models,  the  tail 
position  can  be  controlled  by  two  different  means.  The 
first  is  simply  to  require  that  the  horizontal  tail  rotation 
be  fixed  at  some  value  Sht •  second  is  to  fix  the  input 

commands  Se5x  s-nd  and  then  let  the  stability 

augmentation  system  adjust  the  tail  setting  Sht. 

A.  Mathematical  Description 

Figure  62  shows  a  control  system  representation  of  the 
stability  augmentation  system, 


Figure  62  Stability  Augmentation  System 


Using  figure  62  as  a  guide,  the  following  equations  describe 
the  stability  augmentation  system. 

(8.1)  SsT^  =  Gpp  Fpx 

Where  is  the  force  exerted  by  the  pilot  on  the  stick, 

(8 . 2)  Pfc  =  Gca  Sstk 
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Let  Uq  and  Uq<9  be  defined  by 

(8.3)  Ud  = 

(8.4)  U.<}<?  =  f  Uddt 

Also,  let  Urj  and  Uerj  be  defined  by 

(8.5)  Ue,  =  /R.dt 

(8.6)  /uR,dt 


Then 


(8.7)  “  (  U,QQ  +  B|o  U<^  +  "  UcEJ  "B|i  Upi  )/B; 


(8.8)  R,-  Rfc 

Let  U-Ei  ,  Uj3  ,  IJ-KE3  ,  U'Rf  ,  and  be  defined  by 

(8.9)  Upi’/Pzat 

(8.10)  fz^dt 

(8.11)  Uej23  - 

(8.12)  LLe4.  =  jk^dt 

(8.13)  Uer*  =  /Ue^dt 

Then 

(8.14)  R3  =  (^P  U.RJ  “  Ug3  )/ B;3 

(8.15)  R4.  —  (  lifZfZS  U,EE4  “  8)4  i^Q.4-')^  824. 

(8.16)  Rs-  ~  R4.  +  Sesp  ~  Sgp,^ 

Because  of  rate  and  position  limits,  the  equations 
that  describe  S^t  terms  of  R5  must  be  modified  to 
reflect  these  limits.  Let  Shtp  defined  by 

(8.  17)  SuTp  =  (  Rfl-  -  Sht)/  ®I6H 


Then 


(8 . 18)  Shti  ~ 


(8.19)  5„t 


Shtpmx 

ii  Shtp  ^  ^HTP/vix 

If  -Shtpmx  -Shtp  - 

"  ^HTPMX 

cf  SttTP  "  Shtomx 

min  {  0.0^  Shti  } 

if  Smt  ^  Smtmjix 

^HTI 

ii  SHTMtw  <  Sht*^  Sht 

max^o.OjSrtTiS 

?>  Sht  ^  5n7miN 

Finally  the  stick  position  positioned  as 

a  function  of  time  by  specifying  two  times  and  two  loads 


(8.20) 


Fpxi 

Fpxa 

Fpxi  *»*  ^  Fpxi 


li  Tpx^  -  *Fpxi 

1+  T  ^  Tpxi  ^  "Fpx^ 

^pxi  Tpxz  - 

-  Ppxi )CT-Tpxi  VCTpxz“Tpxi) 

ii  Tpxi  <>  T  -c  Tpxz 


B.  Parameter  Evaluation 

The  values  for  the  F-llU  system  parameters  are  listed 
in  Tuble22.  In  using  this  system  for  a  braking  problem 
the  usual  procedure  is  to  first  choose  a  steadystate 
value  for  Sht(5^ts»).  Set  Smts  -  *  >BST  and  set 

Fesro  (Tpxl  -  Ipj-Z  =  .  .  . 

Set  all  other  initial  conditions  to  zero. 
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Table  22  Horizontal  Tati  Control  Parameters 
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Table  22  (Contd) 
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9a.  WmAY  SYSTEM  (FLY\raEEL  AND  3  DEGREE) 


This  runway  system  is  essentially  the  same  as  the  runway 
system  for  the  6  degree.  In  fact,  the  relation  between 
the  t.70  is  given  by  and 

.  Even  though  this  system  is  like 
the  6  degree  system,  the  equations  are  listed  below  which 
take  advantage  of  the  fact  that  it  takes  less  computer 
time  to  calculate  ZaD<'>^>  than  . 

The  data  describing  the  runway  is  in  tabular  form  and 
consists  of  runway  elevation  values  as  shown  in  table  24 
as  described  in  discussion  of  the  6  degree  runway  system. 

The  data  is  from  the  center  strip  from  station  4574  to 
station  6574. 

A.  Mathematical  Description 

Let  ^  =  I >  2,  •  j  looi  denote  the  elevations 

at  two  foot  intervals.  As  an  example,  CS)  =■ 

If  X  is  a  distance  measured  down  the  runway  where  x 
is  in  inches,  then  z  -  and  u> « 

correspond  to  the  elevation  in  inches  and  the  slope  in 
inches  per  inch.  The  values  for  z  and  uj  are  determined 
as  outlined  below.  The  function  Z&d  will  have  the  pro¬ 
perty  that  Z(s.i><o>  o. 

Let  Xlho  3  constant  such  that  zooo. 

The  input  Xlrf  feet  is  derived  from  %  such 
that  o  ^  Xiffp  <  2000  and  for  some  integer  k 

(9a,  i)  Xlrf  “  Xlro  ^2.  “  2000  k 

Let  n  be  an  integer  such  that  rj  Xlpo  ^ 

and  define  Z^co  by 

(9a.  2)  Zsco"  -^(HpcCn  +  i)  -  ))(X  LEO  2h  +  2  )/2 

If  m  is  an  integer  such  that  2(m-j)-XLEp  <  2m 
then  2  and  w  are  given  by 

(9a, 3)  2  -  12  "  Hc.c(nn)j(XLEF -2m+ Z  K-  ^c&o 

(9a. 4)  uj  ==  2 
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Parameters  (Flywheel  and  3  Degree) 


9b.  RUNWAY  SYSTEM  (6  DEGREE) 


The  runway  system  is  not  actually  a  "system"  in  the  same 
sense  as  the  brake  system,  for  example.  The  runway  system 
is  simply  a  function  called  by  the  airplane  system  to 
supply  values  for  ground  slope  and  elevation.  The  data 
describing  the  runway  is  in  tabular  form  and  consists  of 
runway  elevation  values  as  shown  in  Table  24 ,  Except  for 
a  slight  modification,  the  data  in  Table  24  is  taken  from 
station  4574  to  station  6574  of  runway  25  from  reference  11. 
The  left  elevations  and  right  elevations  are  10  ft.  to  the 
left  and  10  ft,  to  the  right  of  center  respectively.  The 
elevations  have  been  modified  slightly  so  that  the  eleva¬ 
tions  at  station  4574  match  those  as  station  6574,  This 
is  done  to  provide  an  essentially  "endless"  runway  by 
repeated  use  of  a  ba'^ic  2000  ft.  strip, 

A,  Mathematical  Description 

Let  1,23  •••jiool  denote  the 

elevations  at  two  foot  increments  of  the  right,  center, 
and  left  runway  strips  respectively.  As  an  example, 

“  9,550  and  HccC5)  =  9.686,  If  x  is  a  distance 
measured  down  the  runway,  and  v  is  a  distance  measured 
out  from  the  center  of  the  runway  where  %  and  if  are  in 
inches,  then  z  =  u)  = 

correspond  to  the  elevation  in  inches  and  the  slope  in 
inches  per  inch.  The  values  for  t  and  wj  are  determined 
as  outlined  below.  The  function  will  be  chosen  in  such 
a  way  that  Z6d<o,^>  =  0.0  inches. 

Let  and  be  the  inputs  in  feet.  Thus, 

(9c. 1)  -  15/. Z 

Let/.LRobe  a  consta  nt  such  that  0  ^  X^eo  <  2000,  Xlec 
is  a  number  such  that  0^  Xi-ef  2000  and  which  also  satis¬ 
fies  the  following  equation  for  some  integer  K. 

(9c. 2)  Xlbf  =  X\.Ro  +  -  2000  K 
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Now  let  n  be  an  integer  such  that  2(n  -1)  i  Xlro  <  . 

Define  ^sbo  ^  ^ &co ,  ^  slo  as  follows; 

(9c. 3)  Esro  =  He«<n)*g(Hee(«<-i^-HEe<"'')^Xi.ec,-2n  +  2) 

(9c.4)  Zsco  =  HEcWi-i(HscC«4i)-HEcC'i'))(X,jio-2n  +  2) 

(9c.5)ZaL0  -  Hg|_Cn)  +  4  CHeJo+i^- Hg^_(nO(  Xleo' 2n  +  2  ) 

Now  let  vn  be  an  integer  such  that  2(v>i-l)  ^  X^ef  <  2m 
Define  ^Gcx  5  and  as  follows; 

(9c. 6)  =  Hj2(2k.nij  +  2  -  HeeC>^)}(X(_2f  -2m  +  2)  “  ^&fLo 

(9c.7)26cx  "  ■’■  2  + 

(9c.8)^ux  ”  "HeLCm))(XLZF  +2.)  "  ^slo 

If  Xlep^o  ,  then 

(9c. 10)  2  =  12  (  C 2&RX  ~  jCVuEp/lo) 

(9c.ll)  -  (  Vlef  /  20)  (  Hee  C'^t« )  -  H  ERt-no)  “  m+ll  -I-  H 

If  )  then 

(9c.  12)  H  =■  i2-(  ^Qcx  ^^scx^Z-^LX  )^'XLe  p  /lo  ) 

(9c. 13)  uj  =  (  yi.eF/2o)( 

■*■  i  (  H  EC  —  Hr^c  0 


Table  24  Three  Track  Elevation  Profiles 
C  Stations  and  elevations  are  in  feet  ) 


STATION 

ELEVATION 

LEFT 

CENTER 

1 

RIGHT 

Table  24  (Contd) 


Elevation 

Center 

RlSht 

9.667 

9.567 

9.668 

9.564 

9.668 

9.561 

9.6  76 

9.568 

9.676 

9.567 

9.564 

9.564 

9.66C 

9.557 

9.577 

9.561 

9.676 

9.561 

9.677 

9.567 

9.678 

9.574 

9.678 

9.575 

9.683 

9,574 

9.686 

9.^80 

'9.582 

9.536 

9.7Ci 

9.  587 

9.  7C  4 

9.  594 

9.7CJ 

9.597 

9.  7C2 

9.  588 

9.697 

9.590 

9.696 

9.5  86 

9.699 

9.593 

9.  737 

9.588 

9.698 

9.585 

9.694 

9.579 

9.697 

9.  581 

9.698 

9.593 

,9. 700 

9.583 

9.693 

9  579 

9.686 

9.579 

9.679 

9.574 

9.6  76 

9.574 

9.680 

9.568 

9.686 

9.571 

9.6f9 

9.68V 

9.693 

BriU;vl 

9.596 

9i594 

9.698 

9,599 

9,7'>i 

9.6f)3 

9.  TJO 

9.598 

9.704 

9.  596 

9.691 

9.597 

9.689 

9;596 

9.696 

9.595 

9.694 

9.  598 

9.6  83 

9.595 

9.693 

9,  595 

9.696 

9.593 

9,694 

9,587 

9,690 

r  9.505 

9.686 

9.595 

9.583 

9.602 

9.680 

L  9.601 

9^.692 

9.603 

6.  AAS 

9.630 

9,685 

9.  598 

9.678 

9.596 

9.677 

9.  595 

9.673 

9.594 

9.685 

9.  594 

9.698 

9.595 

9.6  93 

9.594 

9.699 

9.596 

9.704 

9.599 

9.703 

9.598 

9.719 

9.597 

9.715 

9.  596 

9.715 

9  .533 

Elevation  j 

Center 

Right 

9.  716 

9.603 

9.714 

9.582 

9.  717 

9.  585 

9.714 

9.580 

9.  711 

9.  5  84 

9.716 

9.578 

9.  713 

9.585 

9.  712 

9.  581 

9.711 

9.585 

9.710 

9.  592 

9.705 

9.585 

9.696 

9.  578 

9.685 

9.576 

9.  695 

9.  572 

9.70! 

9.571 

9.  702 

9.  573 

9.711 

9.  571 

9.  707 

9.568 

9.  699 

9.  569 

9.692 

9.567 

9.684 

9.  566 

9.687 

9.573 

9.  688 

9.  566 

9.682 

9.562 

9.6  92 

9.  5  67 

9.679 

9.  569 

9.678 

9.578 

9.  687 

9.  578 

9.682 

9.579 

9.689 

9.585 

9.686 

9.581 

9.678 

9.585 

9.695 

9.589 

9.6  99 

9.588 

9.709 

9.587 

9.696 

9.  585 

9.685 

9.  585 

9.682 

9.584 

9.690 

9.579 

9.685 

9.583 

9.590  9.694 
9.509  9,692 
9.580  9.691 
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9.917 

9.926 
9.921 
9.931 
9. 937 
9.936 

9.926 
9.939 
“9.946 
9.953 
9.95? 
9.956 
9.953 
?t.96l 
9.973 
9.972 

9.978 
9.977 

9.979 
9.  9S6 
9.993 
9.999 

10.'3n6 

10.308 

10.011 

10.317 

10.021 

lOO?3 

.10.091 

1)*)43 

10.039 

13.039 

10.096 

10.039 

10.033 

10.333 

io.')2r 
10.328 
J0,037 
13.399 
10. 099 
10, JS? 

I  4%  >9 

13.15? 

10.052 


10.028 

10.332 

10.029 

10.324 


13.3  36 
10,04? 
1 3. 1 4  4. 
10.044 


10.  ‘>55 
10.156 
10.348 
10.950 

10.354 
ir.fi49 

13.354 
10.053 
13.343 
10. 041 
13.343 


Right 

9. 

862 

_5. 

361 

9. 

862 

9. 

366 

9. 

S76 

9. 

875 

9. 

874 

9. 

376 

B A  f  An 

1  Elevation 

Left 

Center  . 

Right 

5740 

9.964 

11.058 

9.944 

5742 

9.967 

10.052 

9.943 

5749 

9.964 

10.056 

9.946 

5746 

9.962 

10.049 

9.946 

5749 

9.970: 

10.056  1 

9.959 

5750 

9.972 

10.(J64 

9.  966 

5752 

9.972 

10.089 

9.974 

5754 

9.967 

10.088 

9.977 

5756 

9.966' 

10.092 

9.981 

5758 

9.974 

10.091 

9.983 

5760 

9.982| 

10.098 

9.985 

5762 

9.980i 

10.101 

9.982 

5764 

9.9761 

10.C97 

9.985 

5  766 

9.975! 

10.P98 

9.983 

5768 

9.979 

10.100 

9.976 

5770 

9.978 

10.087 

9.975 

5772 

9.967 

10.094 

9.976 

5774 

9.970 

10.087 

9.975 

5776 

9.970 

10.092 

9.974 

5778  1 

9.967 

10.097 

9.975 

5780 

9.963 

10. 095 

9.976 

5782 

9.963 

10.089 

9.979 

5784 

9.957 

10.  C84 

9.976 

5786 

9.943 

10.377 

9.977 

5788 

9.934 

10.067 

9.9  79 

5790 

9.930 

10.067 

9.983 

5792 

9.927 

10.  064 

9.982 

5794 

9.914 

11.062 

9.980 

5796 

9.915 

10.066 

9.978 

5798 

9.912 

10.074 

9.  975 

5600 

9.910 

10.071 

9.976 

5802 

9.973 

10.073 

9.984 

5804 

9.925 

10.091 

9.987 

5806 

9.927 

10.088 

9.994 

5909 

9.933 

10.105 

in.oci 

5910 

9.942 

10.131 

10.1)02 

5812 

9.947 

10.125 

9.994 

5914 

9.963 

10.123 

9.991 

5316 

9.965 

10.108 

9.984 

5618 

9.973 

10.113 

9.987 

5820 

9.970 

10.134 

9.986 

5822 

9.975 

10.101 

9.985 

5024 

9.970 

10.095 

9.986 

5S?<) 

9.967 

10.  091 

9.987 

5828 

9.963 

10.095 

9.999 

5830 

9.961 

10.  105 

10.004 

5832 

9.963 

10.13? 

13.311 

5834 

9.960 

10.118 

10.015 

5836 

9.959 

13.121 

10.014 

5838 

9.963 

10.115 

10.013 

5843 

9.960 

13. in 

10.017 

5842 

9.967 

10.105 

10.316 

5844 

9.963 

10.103 

10.CC9 

5846 

9.957 

10.115 

!  10.304 

5848 

9.964 

,  10.112 

10.30  .1 

5850 

9  977. 

1  10.103 

9.996 

5852 

•  .985 

1  10.  II5 

9.995 

5859 

9.990 

'  10. Ill 

9.993 

5856 

:  9.991 

;  l').105 

9.993 

5858 

9.983 

1  10.105 

9.994 

5860 

9.987 

;  10.  100 

10.  103 

5862 

9.999 

10.095 

10.303 

5864 

9.983 

10.091 

9.999 

5866 

9.971 

:  10.088 

9.992 

5868 

9.962 

10.082 

9.983 

5870 

9.957 

10.081 

9.  985 

5872 

9.952 

13.088 

9.979 

5874 

9.950 

,  10.073 

9.076 

^Sn 

10.064 
'  10.061 

9./j5 
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station 

Elevation 

Left 

Center 

Right 

5880 

9,945 

13.055 

9.957 

5882 

9,943 

13.041 

9.954 

5884 

9,94  7 

13.025 

9.947 

58  8  6 

9,943 

U.026 

9.  946 

5883 

9,945 

13.028 

9.952 

589C 

9,952 

1  ).021 

9. 957 

5892 

9,957 

13.038 

9.967 

5894 

9,  961 

13.053 

0.  98-3 

5896 

9.963 

13.0  65 

9.994 

5898 

9.971 

13.071 

13.300 

590r 

9.973 

13.0  75 

10.336 

5972 

9.983 

1 3.0  72 

11.3)8 

59(  4 

1C. 032 

10.081 

10.C09 

59{  6 

13.C78 

13.089 

13.337 

5908 

1C. 017 

13. 1-32 

10.003 

5913 

n.C32 

13.112 

9.997 

5912 

10,027 

13.114 

13.032 

5914 

13.039 

10.117 

10, 3  )4 

5916 

10.038 

10.116 

13.3)8 

5018 

10.040 

13.118 

10.010 

5920 

13.041 

13.118 

13.335 

5922 

10.C38 

13. 109 

10.001 

5924 

10.C3i» 

13.123 

1).3J2 

5926 

10.013 

10.113 

9.  998 

5928 

9.994 

13.105 

10,030 

593'‘ 

9.903 

10.111 

13.334 

5932 

9,930 

13.110 

13.0)4 

5934 

9.937 

10.098 

10.333 

5936 

9.987 

10.134 

9.997 

5938 

9.93) 

13.099 

9.99) 

594'* 

9.982 

13.1  93 

9.  991 

5942 

0.983 

13.091 

9.936 

5944 

9.983 

10,099 

9.930 

5946 

9.976 

13.098 

9.975 

5948 

9.  982 

13.092 

9.980 

5953 

9.993 

13.091 

9.984 

5952 

10.033 

13.395 

9.907 

5954 

13.3)8 

n.ioi 

9.992 

5956 

10.009 

13.125 

9.997 

5958 

13.310 

13.123 

13.310 

596'' 

10.033 

1C. 135 

13.  3)7 

5962 

n.OOf! 

10.138 

10.002 

5964 

0.995 

13.128 

9.997 

5966 

9.994 

13.115 

9.988 

5963 

9.983 

13.133 

9.934 

5970 

9.987 

iC.  1-33 

9.977 

5972 

9.937 

1 3  .0  94 

9.974 

59  74 

9.984 

n.f  97 

9.963 

5976 

9.937 

13.099 

9.959 

5978 

9.991 

13.095 

9.954 

5980 

9.993 

IC.f  98 

9.95-: 

5982 

9.992 

13.095 

9.96) 

59  8  4 

9.99  7 

10.090 

9.  947 

5986 

9.995 

13.092 

9.955 

5988 

9.98  7 

10.39-V 

9.957 

5993 

9  .934 

13. C3? 

9.963 

5992 

9.  983 

13.383 

9.967 

5994 

9.976 

13.C38 

9.966 

5996 

9.  976 

13.088 

9.964 

5998 

9.975 

10.G86 

9.965 

60ro 

9.977 

13.075 

9.964 

6102 

9.  9M> 

n.C71 

8.  95  9 

6074 

9.95) 

13.373 

9.958 

6006 

9.945 

13.073 

9.965 

6308 

9.95) 

13.0  71 

9.963 

6010 

9.956 

13.068 

9.977 

6C12 

9.955 

13.0  77 

9.981 

6014 

9.  963 

13.0  79 

9.9)4 

6'‘  16 

9.967 

13.083 

9.  987 

6018 

9.0S6 

I3,0  7:i 

9.986 

XatloD 

- - - 

Elevation 

left 

Center 

Right 

63  2? 

9.962 

IO.C82 

9.979 

6(-22 

9,954 

10.083 

9.975 

6)24 

9.953 

IC.386 

9.972 

6026 

9.953 

10. 083 

9.975 

6328 

9.957 

10.076 

9.971 

6333 

9.960 

1C. 385 

9.973 

Fifi 

9.963 

10.396 

9.  972 

inn 

9.969 

10.392 

9.969 

6036 

9. 973 

10.096 

9.  964 

6333 

9,973 

13.393 

9.960 

6C43 

9.973 

10.089 

9.959 

6)42 

9.973 

13. ■'89 

9.958 

6044 

9.  9  76 

10.087 

9.955 

6)46 

9.973 

10.088 

9.955 

6043 

9.973 

13.034 

9.951 

6)53 

9.972 

10.034 

9.  946 

6)52 

9.933 

10. "73 

9.946 

6054 

9.931 

10.U82 

9.  946 

6366 

9.992 

13.189 

9.951 

6058 

9.  938 

10.093 

9.949 

6363 

9.9?a 

10.092 

9.954 

6C62 

9.930 

10.391 

9.952 

6364 

9.979 

10.091 

9.942 

6366 

9.  970 

10.393 

9.943 

6063 

9.970 

IO.OB2 

9.  942 

6373 

9.973 

10.376 

9.955 

6072 

9.976 

10.083 

9.963 

S374 

9.977 

10.086 

9.968 

6076 

9,933 

10.088 

9.973 

6373 

9.933 

13.188 

9.981 

608-) 

9.930 

10.133 

9.983 

6382 

9.933 

10. 108 

9.986 

6004 

9.937 

10.396 

9.935 

6036 

3.9  90 

10.089 

9.996 

6388 

9.933 

13.392 

9.998 

6093 

9.  931 

10.003 

9.996 

6392 

9.977 

10.396 

9.996 

6<-94 

9.976 

10.089 

9.991 

6)96 

9,964 

13.392 

9.982 

60  98 

9.  965 

10.063 

9.983 

6139 

9.953 

13.335 

9.936 

61  >32 

9.956 

10.076 

9,968 

6104 

9.941 

10.081 

9.938 

61U6 

9.957 

10.379 

9.996 

6133 

9.955 

1C. 076 

9.987 

6irt 

9.952 

10.374 

9.938 

6112 

9.946 

10.072 

9.  995 

6114 

9.944 

10.377 

9.932 

6116 

9.943 

10.075 

9.975 

6113 

9.9'.l 

13.371 

9.977 

6120 

9.941 

10.079 

9.975 

612? 

9.9)8 

10.063 

9.963 

6124 

9.9)8 

10.37? 

9.969 

6126 

9.9)5 

10.074 

9.968 

6128 

9.  940 

10.373 

9.972 

6130 

9.936 

10.072 

9,977 

6132 

9.9)3 

11.176 

9.974 

6134 

9.  937 

10.059 

9.970 

6136 

9.943 

n.377 

9.969 

6133 

9.  94f- 

10.070 

9.967 

6143 

9.9)4 

10.061 

9.970 

6142 

9.9)5 

10.373 

9.971 

6144 

9.937 

10.066 

9.  9  72 

6146 

9.93) 

10.?67 

9.977 

6148 

9.9)6 

1C. 066 

9.  981 

6153 

9.9)2 

13.374 

9.987 

tlf2 

9.930 

10.069 

9.996 

6)54 

9.9)3 

13.?66 

10,334 

6156 

9.9)? 

10. J82 

n.))s 

6153 

9. '>36 

13.071 

It'.'. ^ 

station 

Elevation 

Left 

Center 

Right 

6160 

9.937 

10.C79 

lO.Ill 

6162 

9.937 

13.C86 

10.012 

6164 

9.  938 

10.078 

10.015 

6166 

9.940 

12.396 

11. -"U 

6163 

9.940 

13.093 

10,?14 

61T0 

9.938 

10.096 

ic.r"9 

6172 

9.  943 

13.076 

9.999 

6174 

9.945 

10.085 

io.ur-6 

6176 

9.942 

13.080 

13.)35 

6178 

9.939 

lO.C  79 

i(!.ro2 

6181 

9.942 

13.170 

io.:iC2 

6182 

9.  938 

11. 062 

9.999 

6134 

9.942 

11 .059 

9.995 

6136 

9.955 

10.  063 

9.995 

6188 

9.946 

13.057 

9.986 

6190 

9.943 

10.065 

9.985 

6192 

9.934 

13.062 

9.  973 

6194 

9.936 

13.049 

9.962 

6196 

9.93C 

10.043 

9.955 

6198 

9.929 

13.L59 

9.954 

6200 

9.92  7 

10.C49 

9.954 

6202 

9.929 

13.081 

9.955 

6204 

9.922 

10.075 

9.  949 

6206 

9.914 

13. {’71 

9.946 

6203 

9.912 

1U.C61 

9.946 

6210 

9. 90S 

10,071 

9.949 

6212 

9.  897 

10.060 

9.953 

6214 

9.894 

10. 058 

9.949 

6216 

9.887 

13.362 

9.949 

6218 

9.889 

1  ).C61 

9.944 

6220 

9.887 

13.063 

9.941 

6222 

9.886 

10.068 

9.935 

6224 

9.880 

10.061 

9.929 

6226 

9.877 

10.055 

9.925 

6228 

9.873 

10.063 

9.  916 

6233 

9.867 

10.061 

9.914 

6232 

9.860 

10.055 

9.912 

6234 

9.853 

13.?43 

9.916 

6236 

9.850 

10.044 

9.922 

6238 

9.846 

13.041 

9.921 

6240 

9.841 

13.031 

9.719 

6242 

9.844 

13  .021 

9.91b 

6244 

9.854 

1.3.015 

9,909 

6246 

9.348 

13.013 

9,903 

6248 

9.843 

9.997 

9.99) 

625G 

9.844 

10.009 

9.899 

6252 

9.840 

10.018 

9.996 

6254 

9.836 

13.(  13 

9.  899 

6256 

9.833 

lO.t‘22 

9.91/3 

6258 

9.827 

10. Cl ) 

9.899 

6263 

9.323 

13.301 

9.891 

6262 

9.  321 

11.004 

9.  886 

6264 

9.321 

9.993 

9.99) 

6266 

9.826 

9.956 

9.879 

6268 

9.328 

9.958 

9.891 

62  73 

9.3  30 

9.963 

9.88? 

6272 

9.8)C 

9.964 

9,  9  91 

6274 

9.331 

9. 972 

9.979 

62  . 

9,832 

9.973 

9.  9S4 

62  78 

9.3  )3 

■».934 

9,995 

628  • 

9.827 

9.  9.34 

9,?  95 

6282 

9.82) 

9.974 

9.993 

6284 

9.824 

9.  973 

9.991 

6286 

9.824 

9,967 

9.  99u 

6288 

9.82) 

9.958 

9.99; 

6290 

9.82? 

9.949 

9.  979 

62  92 

9.315 

9.957 

9.97* 

6294 

9.  a'-b 

9.  952 

9.  36) 

6296 

9.837 

9 . 9  8 

9.955 

Hh 

9.  316 

9.  9-,  1 

- 

9.  -15  ) 
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Table  25  Runway  System  Parameters  (6  Degree) 
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'A’  This  input  allows  starting  the  airplane  on  a  different  part  of  the  runway 
profile  even  though  its  distance  down  the  runway  is  the  same. 


SECTION  IV 


TOTAL  SYSTEM  ANALYSIS 


Three  different  total  system  mathematical  models  have 
been  formulated  to  perform  antiskid  analysis.  The  first 
model  which  is  referred  to  as  the  flywheel  system  repre¬ 
sents  an  antiskid  system  installed  on  a  wheel  and  brake 
which  are  mounted  on  a  dynamometer.  The  second  system, 
referred  to  as  the  three  degree  system,  represents  an 
antiskid  system  installed  on  a  wheel  and  brake  mounted  on 
a  rigid  airplane  which  is  allowed  three  degrees  of  freedom 
(longitudinal  translation  down  the  runway,  translation 
vertically,  and  pitch  rotation).  The  third  system,  referred 
to  as  the  six  degree  system,  represents  a  rigid  airplane 
having  all  six  degrees  of  freedom  and  equipped  with  a  con¬ 
ventional  single  wheeled  main  landing  gear  incorporating 
independent  antiskid  control  of  each  brake.  All  of  these 
systems  are  created  utilizing  the  models  described  in 
Section  III.  The  basic  reason  for  utilizing  three  models 
is  economics.  The  six  degree  system  takes  at  least  twice 
as  long  to  run  as  the  flywheel  system  and  not  all  antiskid 
system  parameters  require  the  sophistication  of  the  six 
degree  system.  However,  it  might  be  necessary  to  check 
certain  effects  under  the  most  ccmprehensive  circumstances. 

The  "Basic  Control  System"  is  made  of  the  following 
models  as  described  in  Section  III: 

1.  .?rake  System 

2.  Hydraulic  System 

3.  Wheel  Speed  Sensor 

4.  Control  System 

5.  Antiskid  Control  yzlv& 

To  form  the  flywheel  system,  the  "Basic  Control  System" 
is  combined  with  the  3a.  Airplane  System  (Flywheel),  4a. 
Wheel  and  Tire  System  (Flywheel),  and  the  9a,  Runway  System, 
To  form  the  three  degree  system,  the  "Basic  Control  System" 
is  combined  with  the  3b.  Airplane  System  (3  Degree),  4b. 
Wheel  and  Tire  System  (3  Degree),  8.  Horizontal  Tail  Con¬ 
trol  System,  and  9a.  Runway  System,  The  six  degree  system 
incorporates  two  separate  "Basic  Control  Systems"  and  two 
separate  4c.  Wheel  and  Tire  Systems,  These  are  combined 
with  a  3c,  Airplane  System  (6  Degree)  which  utilizes  the 
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8.  Horizontal  Tail  Control  and  9c.  Runway  System  (6  Degree). 
The  model  flow  diagrams  are  shown  in  Figures  64  ,  65  ,  and 

66  .  When  utilizing  the  "Basic  Control  Systems"  with  the 

six  degree  system,  the  variables  communicating  with  the 
airplane  model  are  reidentified  to  correspond  to  the  right 
or  left  side  of  tho  airplane.  Thus,  Xaxr  is  Xax  in  the 
right  side  and  is  X^x  in  the  left  side. 

The  high  degree  of  modularity  used  in  this  analysis  is 
desirable  for  three  reasons.  The  first  reason  is  that  it 
is  easy  to  combine  the  component  models  together  to  form 
different  types  of  overall  systems.  This  is  true  not  only 
from  modeling  considerations  but  especially  from  program¬ 
ming  aspects.  As  an  example,  the  only  basic  change  re¬ 
quired  to  acconrmodate  a  twin  or  tandem  gear  would  be  to 
remodel  the  strut  in  the  airplane  system.  The  second 
reason  for  modularity  is  the  difficulty  in  being  completely 
general.  Should  a  component  arise  which  is  not  described 
by  the  existing  models,  it  is  easy  to  create  a  new  program 
for  the  new  model  without  having  to  modify  the  operation 
of  other  systems.  Thus,  from  the  programming  point  of  view, 
to  incorporate  a  new  wheel  speed  sensor  for  example,  the 
new  model  program  can  fall  back  on  the  existing  read,  write, 
and  logic  statements  of  the  existing  wheel  speed  model.  The 
input  and  output  variables  of  the  new  ca^iponent  model  are 
automatically  incorporated  properly  into  the  overall  compu¬ 
tational  procedure,  unless  some  new  variables  are  defined. 
The  third  reason  for  using  a  modular  approach  is  to  take 
advantage  of  the  different  response  characteristics  of  the 
different  component  models.  In  the  digital  procedures  uti¬ 
lized  for  computation,  essentially  a  "fixed  step"  integra¬ 
tion  technique  is  employed.  The  step  size  used  in  a  "slow" 
responding  model  does  not  have  to  be  as  small  as  one  used 
in  a  "fast"  model.  Thus,  different  component  models  uti¬ 
lize  different  size  integration  steps  which  minimizes  the 
overall  time  of  computation. 


_ _  9c ,  Runway 

3c.  Airplane  - - - \  System 


Figure  66  Six  Degree  System 


SECTION  V 


SAMPLE  CASE  ANALYSIS 


A  computer  solution  of  the  composite  total 
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solution  of  the  same  equations. 
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APPENDIX  I 


DERIVATION  OF  EQUATIONS  DESCRIBING 
THE  OPERATION  OF  THE  GOODYEAR  ADAPTIVE 
ELECTRONIC  ANTISKID  CONTROL  CIRCUIT 

The  mathematical  description  of  the  operation  of  the 
Goodyear  adaptive  electronic  antiskid  control  circuit  as 
shown  on  Figure  50  is  developed  with  conventional  cir¬ 
cuit  analysis  techniques  using  Kirchhoff’s  Laws.  Figure 
68  is  the  schematic  diagram  from  Figure  50  with  the 
transistors  and  diodes  shown  in  terms  of  their  equivalent 
circuits  and  the  various  currents  and  voltages  identified. 
The  transistor  and  diode  equivalent  circuits  are  adapta¬ 
tions  of  equivalent  circuits  developed  and  described  in 
references  13,  14  and  15.  Some  of  the  diode  forward  resis¬ 
tances  are  combined  with  other  resistance  in  series  with 
the  diodes  and  are  not  shown  separately.  Also,  since  the 
current  through  Rq  (the  output  resistance  of  the  wheel 
speed  signal  source)  has  three  non-mutually  influencing 
components,  Rci  is  included  in  Ra  ,  Ro&  and  Ru  to  simplify 
equations.  Other  simplifications  will  be  described  and 
discussed  during  the  development  of  equations. 

Referring  to  Figure  68  ,  the  circuit  equations  are 
developed  as  follows: 

The  voltage  across  capacitor  Ct  is  defined  as: 

VC|  -  y* VCi  dt 

(Al)  Vc\  '=  on  \/ci  - /}c\  CCoOS 


/Jci  is  the  current  through  G  and  Ct  is  the  c^acitance. 
/Icj  is  established  by  summing  currents  at  node  w/h  as: 


(Nil)  Aoe- A>Z. 


Using  Ohm's  law  and  summing  voltages  around  the  loop 
of  which  Rd0  is  a  part,  /?09  is  established  as: 

(RIO)  WpS  -  (^G-Vci'-Eo9)/lio9  FofiiE6-^ci-toi)>0 

-  o  Fon  iEiS - iiii -A?/ o 
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68  Modulated  Antiskid  Schematic  with  Mathematical  Identification 

and  Incorporating  Equivalent  Circuits  for  Transistors  and  Diodes 


Figure  68  (Contd) 


Noting  that  because  of  diode  Og,  is  restricted  to 
positive  values  only. 

To  combine  constants,  write  equation  (RIO)  as; 

(RIO) 

'0  /^/L  Cc2d^CbZ/ 


By  summing  currents  at  Node 
as : 


current  is  established 


(N12)  /)u  -  /Ifl3  t  /)cai 


Substituting  equation  (N12)  into  equation  (Nil)  gives: 
(Nil) '  /4:/  =  /}oe '/?/??  -  ^cci I 


To  compute  ftcbli  it  is  desirable  to  first  obtain  equations 
for  the  voltages  at  the  base  and  emitter  of  6i\  ii:  terms  of 
the  base  and  emitter  currents  and  the  appropriate  voltage 
sources.  The  vol^ge  at  the  base  of  is  ijfZ\  .  Summing 
currents  at  Node  establishes  current  as: 

(HI)  /iR\= 


Summing  voltages  around  the  loop  to  l/c’z 

gives;  ^  \Jcz  •  If  the  current  through 

diode  W  is  assumed  negligible,  then  - 

(Because  of  the  relative  resistances,  the  current  through 
O4.  is  a  very  small  fraction  of  the  current  through  flse') . 
By  Ohm's  law,  ^  LUSB  and 

yR.1  ~  /^Ri  /l\  >  Substitution  into  equation  (Nl)  and 

solving  for  gives:  j 


(Nl)' 


Vc, 


I  +- 


/3d/i  r^QH 


)/(- 
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The  voltage  at  Node  (^)  is  ^OZf-}koz  .  (Here  two  diodes 
Oz  and  O3  as  shown  on  Figure  50  are  combined  in  an 
equivalent  single  diode).  Summing  currents  at  Node  my 
establishes  current  fioz  as.' 

(N2) :  /?02  -  j9m.4  A/iS  -  /J^ai 


By  Ohm*s  law,  =  ,  and 

i.  kiiAi-  (Note:  Because  of  the 

relative  resistances,  is  a  very  small  component  of 

the  current  in  Rsb  and  is  assiimed  to  be  zero  when  computing 
\Iaox  Summing  voltages  around  the  appropriate  loops 
establishes  that  -\/boz.  - toz.  and  V/Q^A  - 

Eoz.  •  By  substitution  into  equation  (N2)  and 
solving  for  Vfioi  gives: 

my 

(,  60Z  ^4-  j 

Summing  voltages  around  the  loop  thrc  gh  which 
(the  base  of  current £$i)  flows  results  in: 

(Q-1)  ~ 


By  substituting  ,  )/SCl\  *  H.Ba\ 

amd  \l3Ci\  '  A6u\  (From  Ohra*s  Law)  along  with  equations 

(Nl) * ,  (N2) ' ,  and  the  basic  transistor  relationship 

(.hFFifO  ^6(31  equation  (Q-1)  and  solving  for /fe?j 


Cooz~ 


(Q~l)  *  Cq>0Z  -  Cu,o3  i^CZ. 

fe)  ^  '/ 


fi8Q\  ^ 

(hFr, +0  (hfe 


C&oB-. 


Bt  ^R.e 


r,)/0*  '/  i  Lt.*  Hi 
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And  where  in  the  above  (Zs'~  and  =  t  Use  * 

For  (2i  to  operate  as  a  transistor,  l/c«i  must  be  positive. 
Using  equation  ,  the  basic  transistor  characteristic 

(Q2)  /?cci-h«i/fe3i/teH}and  by  writing  the  voltage  loop 
equation  through  Dz  andt^^/i/)  it  can  be  shown 

that  is  positive  for  negative;  therefore, 

equation  (Q-1)  *  is  applicable  only  for  ( ire?;  -  l/dj ) O  • 
Substituting  (Q-l)  into  (Q-2)  gives  the  following  equation: 

(Q-lC)  /Ic^  I  "  C<o04  -  CCiOi  l4‘t  ^  O 

-o  f=vjt  ^i>6^-/C6o:r 


hfFx  r\ 


C 


By  summing  currents  at  node  (^0,  current  is  established 
as : 


(Nio)  /)/?s  =  t  /itn 


To  compute  the  components  of  /?/C3  (i.e, , and )  it 
is  desirable  to  have  equation  (NIO)  in  a  form  where  is 
expressed  in  terms  of  the  appropriate  voltages  and  resis¬ 
tances,  By  summing  voltages  around  the  loop  ,  K.?,  C 

and  GE’M,Vl27is  established  as: 


(V7) 


Substituting  equation  (V7)  along  with  V/C"?  ~ Hs 
and  }ii  into  (NIO)  gives: 

(NIO)  *  r  Al - 

To  combine  constants  write  as: 

(NIO)  '  Aa  =  /fsdx  Cd//.  +  (■  i/cj]  C;,f 


(NIO)  ' 


234 


To  compute  current  ,  Sum  voltages  around  the  loop 

through  which flows: 

(Q-2)  V/ -  -o 


Here  voltage  Vy  is  the  base  voltage  on  Q2  and  is  either 

or  whichever  is  the  largest; 

therefore,  there  will  be  a  version  of  equation  (Q-2)  for 
each  of  these  conditions.  To  establish  which  condition 
exists,  it  is  necessary  to  compute  and 


During  derivation  of  Equation  (N2) *  it  was  observed  that 
-  ykoz.  f  ^£>z  and  it  was  assumed  that  was 

small  when  compared  to  .  Using  the  same  assumption 
is  established  as  follows  by  Ohm’s  Law: 

(R-5)  ^  /?/^S 

By  substituting 

into  equation  (R-5)  above  gives : 

(R-5)*  tl/fiuz  bto  ■2.) 


By  substituting  equation  (N2)  '  for  and  investigating 

the  influence  of  /kai  within  its  allowable  range,  it  can 
be  seen  that  for  practical  purposes  [^5:^  is  a  constant. 

To  compute  VAL  ,  it  can  be  seen  that  equals 

when  yy  is  ;  therefore,  ko  by  Ohm*s 

law. 

By  Ohm's  Law  ~  /?A3 6  .  Since  the  current  through 

D4  is  very  small  and  may  be  assiomed  zero,  and  since 
is  such  a  small  component  of  the  current  through  R1  that 
it  can  be  assumed  zero,  by  Ohm’s  Lav/: 

/^fiQQ  =  Aaq/V  ■-  l/d-'z J  L  l\K  ^ 
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.  “/"'•w'.N,*  "T*  '•*1  \  S\.  ■ 


Therefore: 


i  Vci  -\//ise  -B04)  (  /f/  - 


^9/^rfl8A 

To  establish  whether  Vy  ~  l//Cs'jQ  -  '^6  )  6>£ 


V/ 


=  L  L 


^  i  -f-  _  \  / 


A  voltage  VB  will  be  defined  as  follows: 


(VB) 


(  -  Kg /hciz 


if  V<S  >  <9,  (VY-1) 
if  1(3^0,  (VY-2) 


V  y  -  1  V(:'z  ^  1^4.^ 

Vy  -  L  V^-s''3  " 


where  above  is  defined  as  ; 

-,  / \ 


Before  proceeding  with  the  computation  of  /h^  ,  the  valve 
control  €irDplifier  and  modulation  circuit  elements  will  be 
examined  to  develop  equations  for  and  Vci  ■ 

By  summing  currents  at  Node current  ^0^  is  established 
as : 


(N5)  AdB  -  Ava-Ac^ 

The  voltage  at  Node  is  defined  as  Va  .  Summing  voltages 
around  the  loop  f(^  ,  and  V'x.  gives: 

(vx)  Ey  rtoQy^  \/kos Eds' 


(EV)  -  Aoo  C40h  ^  C40y 
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By  summing  currents  at  Node(^,  ctirrent /5i2/7  is  established 
as ; 


(N4)  /i^ii  -  hv  - 


By  subst;ituting  equations  (N4),  (N5)  ,  and  (EV)  into 
equation  (VX)  and  by  using  Ohm*s  Law  to  establish  that 
^1/  -  /iyJlv ^  /f/?  and  fios' ^  VX  is 

established  as: 

(VX) '  V'x  =  Aos-  C4ti0  I-  C'^OI 


By  substituting  equation  (EV)  into  (n4)  and  using  the 
relationships  i=v^/^r/cy  and  \//in  ^/kn  li\i  EV  is  established 
as : 


(N4)  '  try  =  Oo6  1-  C40I 


The  operation  of  transistors  Q2,  Q3,  Q5,  Q6,  and  Q7  will 
now  be  considerred  to  develop  an  equation  for  current  Ava 
(Valve  Control  Amplifier  Output  Current). 


By  summing  currents  at  Node^^ 
as : 


current  is  established 


(N7)  r  AtL 


By  summing  currents  at  Node  (ll3;,  current  is  established 
as : 


(N13)  t  Am 


By  summing  currents  at  Node  (n^  >  current  is 

established  as: 


(NO) 


l^^io  -  -  /)B6n 


By  sunn'.ing  currents  at  Node^J^  current  /9/?,9  is  established 


as : 


(Ny) 


Accis 
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Summing  voltages  around  the  loop,  and  £.(o 

gives : 

(VlO)  y^dil  f  £(^11  t- l/ffcll 

By  using  the  relationships  \/r/o  ^c€c\ 

and  VdQl  -  ASQ.-y  ^80.'!  established  by  Ohm*s  law 

along  '  ith  the  transistor  characteristic  -  thfe-?  *~ ^)  ^6^1 
and  substitution  equation  (N6)  into  (VlO)  and 

solving  for  As^]  ; 

(VlO)  *  ^  H£ac>  liio  -  _ 

[_  I?tfa7  t  /I I o'J 

By  substituting  (VlO)*  and  (N6)  into  the  relationship 
VA/o  -  A8.to  /Jio 

(AlO)  ^  >4^6  C^oz.  r 


Summing  voltages  around  the  loop  RlO,  Htdi,  ^  ^g<a6. 

AcQg^j  and  /?^  gives: 

(V9)  O  “  14/£)  ^ 

By  substituting  equations  (AlO)  and  (N9)  into  (V9)  along 
with  transistor  characteristics  ~LkFe^r>)  Asdc, 

and  Ae^s  -  /fsas^  and  the  Ohm*s  Law  relation¬ 

ships  V/t^rzA/c^  )  1^6 G  )  1/^(2 

-  Ae^s-  f  and  ]/e^^  -  A  ear  and  solving  for 

AeojU  ' 


(V9)* 


ficcm 


A  0  /cgiit +/?e-43r  ^ 

Li^x  f  AieQ^,  f  77 - +- 


(Upii  (kf^efi)ih!V6’h/) 
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Substituting  equations  (N13),  (N6) ,  (N9),  and  (V9) *  into 
equation  (N7)  along  with  transistor  characteristics 
IhFt'i-'r/)  and  =  LkFiT^  ti) 

and  solving  for  gives: 

(N7)  ’ 

i^a3  C404- 


It  should  be  noted  that  these  operations  involving  Q5,  Q6, 
and  Q7  assume  thatAca?  is  large  enough  such  that  Ai/zh  is 
not  negative  and  that  the  applicable  supply  voltages,  VSl 
and  VS2,  are  large  enough  to  keep  ^  and  )Jci3S' 

positive  at  all  times.  The  latter  assumption  can  be 
proven  to  be  true  for  the  range  of  currents  experienced 
during  circuit  operation.  If  /hik3  is  not  greater  than 
C405/C404,  insufficient  voltage  is  developed  across  R9 
to  cause  Q5,  Q6  and  Q7  to  operate.  For  less  than 

C405/C404  all  of goes  through  R9  and  ~ 

therefore^ equation  (N7) *  has  two  forms  depending  on  the 
value  of  •  Write  these  two  forms  as  follows: 


'  FOt  Ac£i3 

Supply  voltage  VS2  is  large  enough  so  that  voltages  Vci^s 
and  l4di3  are  always  positive  and  a  small  leakage  current 
/)c(^3  0  flows.  All  the  equations  developed  here  are  for 
the  increment  of  At??  above  the  leakage  value. 

By  using  the  transistor  characteristics  Aca3  -  As613 

ActXl.  -  AeQj.  MPci/ IkFtrXt 

And  if  at  Node  N15  Asqs  -  A?Ciiz^ 

then: 

(Q3)  /icQ^s  =  C6?0(c 

L^oh  -  } - - - v — 

Ck.FF'Z.  ^ 

The  operation  of  the  modulating  circuit  element  will  now 
be  examined.  To  compute  valve  voltage  EV  from  equation 
(N4)  *  the  value  of  current  which  is  established  by 
Equation  (N5)  is  required.  Equation  (N5)  shows  that  a 
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ccmponent  of  Ad^  is  Ac<f.  .  Before  developing  equations  for 
computing  /?c4  some  observations  relative  to  the  operation 
of  C4  and  Q4  are  helpful. 

By  summing  currents  at  Node  current  Acf  is  established 
as : 


(N8) 


However,  because  of  diodes  D7  and  D9  currents  and 

have  limitations  depending  upon  the  direction  of  voltage 
across  the  diodes.  Summing  voltages  around  the  loop  C3, 
Rll,  D9,  C4  to  VX  gives: 

(Vll)  r  ~ /x  =  0 


Substituting  equation  (Vll)  into  the  expression 
VAlt-/l/iiilih3LS  established  by  chm's  law  gives: 


Because  of  D9,  fl/in  -  O 


fiJ/L  [/k  -^tOf  j-6 


Summing  voltages  around  the  loop  C3,  R12,  U7,  C4,  to  VX 
gives : 

(V12)  ^Cz  “  -d 


Substituting  equation  (Vl2)  into  the  expression 
established  by  Ohm's  law  gives: 

(V12)  '  /III, I  --  ( Z.  -ton 

X’/^ 

Because  of  D7  //hi-O  poU  i)/cx/l>7tZc<f.- ZA)  —  0 

Summing  voltages  around  the  loop  C3  through  Q4  to  C2  gives: 

(V-Q4)  /(fj -  -0 
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Since  the  currents and  in  transistor  Q4  are 
restricted  to  positive  values  only,  voltages  VBQ4  and  VEQ4  ar 
are  always  positive;  therefore,  equation  (V-Q4)  shows  that 
VC2  is  always  less  than  VC^  by  an  ai^ount  at  least  equal  to 
EQ4.  Also,  because  of  diodes  D7  and  D9,  no  current  can 
flow  from  C3  through  Rll,  D9,  D7  and  R12  to  C2,  For  these 
circumstances,  it  is  observed  (1)  that  positive, 

all  of  Ac4  passes  through  Rll  and  all  of  A/Lil  is 
and  (2)  that  for  negative,  all  of  /?c4  passes  through 
Rl2  and  all  of  is  ~ 

Since  there  cannot  be  positive  and  positive 

simultaneously,  equation  (N8)  evolves  to: 

(N8)'  Ac‘i-  - Ak\i>Ci 

^  'Fon  Ax-ii^a  =  o 

-/tail.  FbH  Alill>0 

By  substituting  equations  (N5) ,  (VX)  *  and  (N7)"  into 
equations  (Vll) *  and  (V12) * ,  equations  for  are 
developed  for  each  case. 

The  remaining  equations  for  the  modulation  circuit  element 
will  now  be  developed.  Substituting  the  expressions 
V/3a4  =  R.8(ii4  and  as  established  by 

Ohra*s  law  along  with  the  transistor  characteristic 

-  {.h  1^4  into  equation  (V-Q4)  and  solving  for 

/)/M4  Sives: 

(V-Q4)  '  ^ 

Foil  /  0 


-  0 


To  combine  constants,  write  equation  (V-Q4)  as: 

(V-Q4)  -  il/c}  -/4'i)  C6  2i  -Co23 

fo/K.  ^ 

-  O  P^A.  (i^j-^2,)  -  ^C^XlJC6  2.^. 

Also,  since  current  AEQ4  is  needed,  define  the  transistor 
characteristic  as  equation  Q4: 

(Q-4) 

where  C,C>/^  ~  (Ap^4-  ^  / ) 

By  summing  currents  at  Node  (N14)  current 
established  as: 


(N14) 


/Icz  = 


Using  the  same  assumption  relative  to  as  was  made 

for  equations  (Nl)  '  and  (V8)  and  by  Ohm's  Law  fiilitAi- 
is  established  by  equation  (V8)  as: 

(V8)  - _ _ (Repeated) 

f /irs 


By  summing  currents  at  Node  (n^ current  fic2  established 
as : 

(N3)  J^thrAn/-)- -  Asm 

Current  is  computed  from  l/c^^  =  /JdiA  £l4- 

established  by  Ohm's  law  and 4  is  computed  from 
equation  (V-Q4) '  and  Equation  (Q-4). 


242 


Equation  (N8) *  establishes  that: 

/tlilZ  -  fOJt 

-  ^  /fc^  i  O 

fo/l  Ac4  >0 

'  O  FO/i  /}c^  ±  0 

Substituting  the  above  and  equation  (V8)  into  equation 
(N14)  gives: 

(NIA)  /let,  -  /leo,  4- 1  Ac4-  -  C 0>lil  Pofi.  Ac4-  ^  6 

-  /'as4  -  C(,/S  Foe.  Aca-  ^  6 


Similarly,  by  substituting  the  above  Afl.n  to 
relationship  and  /€  / f  into  equation  (N3) 

f\oz  i-s  established  as: 

(N3)  '  /Ic3  =  /}c4  -  Ves  A^>o 

-  -  ]/cs  Co/'i  -  Fiifl.  Ac4-  ^  o 

The  voltages  across  capacitors  C2,  C3  and  C4  are  established 
by: 


(2) 

Vc^u.  = 

/  viz  cli 

(A2) 

V/c  z  - 

^6o9  Acz. 

(3) 

Vfs  ■ 

/  M 

(A3) 

\Jc3  - 

C^>/ci  Acs 

(4) 

Vc-f  = 

I  v'cleii 
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(A4)  y^^^ 

All  of  the  equations  describing  the  antiskid  circuit *s 
operation  have  now  been  developed;  however,  to  obtain  a 
computer  solution  of  these  equations,  they  have  to  be 
converted  to  a  suitable  form  so  that  there  are  no  "closed 
loops."  Also,  since  the  equations  tor  ^  c2rkj 

have  different  forms  depending  upon  which  circumstances 
exist,  a  procedure  must  be  established  to  define  which 
form  of  equations  (Q2) ,  (N7)"  and  (N8) *  applies  for  each 
instance.  There  are  twelve  (12)  possible  combinations  of 
circumstances  as  shown  on  Table  18  .  The  procedure  for 

defining  which  condition  exists  will  be  to  assume  a  condi¬ 
tion  and  develop  a  set  of  equations  based  on  the  assixnp- 
tion.  Using  these  equations,  the  assumption  will  be 
tested.  If  the  test  is  affirmative,  the  asstimed  condition 
exists.  If  the  test  is  negative,  the  assumption  is 
incorrect  and  other  assumed  conditions  are  tested  until 
an  affirmative  test  result  is  obtained.  To  illustrate 
this  procedure,  the  equations  for  circuit  condition  4  will 
be  developed: 

For  circuit  condition  4  ^<^4  is  positive greater  than 
VB  greater  than  zero.  Substitute  equation  (N5) 
and  the  applicable  version  of  equation  (N7)"  into  equation 
(VX)  * . 

(VX)  *  -4  V<  =  (  C^<i-  -t  C¥<J/ 


From  equations  (N8)  *  and  (Vll)  *  is  established  as: 
(N8)'-P  ficA- 

/?// 

Substitute  (VX)*-4  into  (N8)*-P  and  solve  for 
(N8) *-P4 

/)c4-  ~  ^ _  (  i/^4-  f  ^  i  C^ijo'Ct^o! J 
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Now  substitute  equations  (N5) ,  (N7)"  and  (N8)*-P^  into 
equation  (n4) *  and  solve  for  EV 


(N4) '-4 


Ev  ^ 


X?  i  I  ^C^oo 


~tp<i 

-f-  C^7  i  ^t\r 


Substituting  equations  (N4)*-4,  (NlO) * ,  (Vy-1), 

and  (Q3)  into  equation  (Q-2)  and  solving  for  gives: 


Substituting  equations  (Q-2) -4  and  (Q3)  into  equation 
(N8)*-P4  gives: 

(N8)  '  -P4  j  ^  (_  -  Ca-  7X' 

Substituting  equations  (Q-2),  (NlO)*,  (N4)*-4,  and  (Vy-2) 
into  equation  (VB)  results  in: 

(VB-3,4)  ]lcxC^  '  (!i47  C4'h^  i)/cs  t'C'icC 

Now  using  equations  (VB)-4,  (Q-2)-4,  (Q3)  and  (N8)*-P4 
the  assumption  that  ^C4.>0t  ^Ci33  >  ^  tiAd  \/b>0 

can  be  tested.  If  the  test  is  affirmative,  then  values 
for  and  EV  can  be  computed.  If  the  test  is 

negative  another  condition  must  be  tested  for. 

Tables  14,  16  and  17  are  a  suanmary  of  test  equations 
developed  in  the  same  manner  as  above.  Since  equation 
(Q3)  establishes  a  linear  relationship  between  /9e'az  and 
Acci.^  and  since  /9sci2,  needs  to  be  computed  as  a  step  in 
the  computation  of/h^y  the  test  equations  for^ca^  will 
be  performed  implicitly  by  computing  Ae^x,  and  comparing 
its  computed  value  to  C607  where  C607  is  defined  as: 

C607  =  -  _ 

(C40iy:6o^) 

Currents  and  are  computed  using  the  applicable 

test  condition  equations. 
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As  shown  on  Figure  49  the  locked  wheel  prevention  cir¬ 
cuit  elements  also  have  an  input  to  the  valve  control 
amplifier.  In  the  equations  thus  far  it  has  beenassumed 
that  the  locked  wheel  skid  signal,  at  node  ^1^,  is 
zero.  When  computing  the  valve  voltage,  it  is  necessary 
that  the  non-zero  value  of  be  accounted  for.  If 

is  not  zero  then  equation  (Q3)  is : 

(Q3-1)  =■  /)lWs 

Since/LiVi  is  a  two  valued  variable  (i.e.  either  zero  or 
the  value  required  to  drive  the  amplifier  as  necessary  to 
achieve  full  brake  release)  insofar  as  valve 
voltage  computation  is  concerned,  it  can  be  considered  as 
a  current  which  can  be  added  to  in  Equation  (Q3) . 

If  we  define  a  current  ,  valve  amplifier  input 

current,  as: 

(LW-i)  Asst  + 

and  tf eat  this  current  like  in  equation  (Q3)  and  if 

we  substitute  equations  (Q3)  and  (N7)"  into  equation  (N5) 
an  equation  for  computing  /ter  is  formulated  for  each  cir¬ 
cuit  condition.  Current  is  then  used  in  equation  (N4) 
to  compute  EV-  Table  15  lists  the  version  of  equation 
(N5)  which  is  to  be  used  for  computing  current  for 
each  circuit  condition. 

Since  the  variables  EG  and  VCl  are  always  used  in  the  form 
of  their  difference,  we  will  define  the  difference  as 
equation  (5) 

(5) 

For  the  cases  where  the  antiskid  control  circuit  mathe¬ 
matical  model  is  used  with  the  flywheel  system  or  three 
dimensional  airplane  system,  the  flywheel  velocity,  VF, 
or  the  airplane  velocity,  X,  as  applicable,  will  be  used 
as  the  airplane  speed  reference  circuit  element.  The 
wheel  speed  comparison  element  will  be  described  as  follows 

(6)  /llWS-  Con  Co/S' 

~  O  Fu/L  0/1  ~  Cipli, 
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